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ABSTRACT 


An experimental and analytical study has been carried out to 
determine the limiting heat transfer mechanisms of dropwise 
condensation. 


Extensive heat transfer measurements were made for drop con- 
densation of saturated "gas free steam" on a vertically oriented 
copper surface. Dioctadecyl disulphide and Teflon were used as 
promoter agents. The parametric effects of AT, pressure, surface 
roughness, vapor velocity, and time were investigated. 


Visual experiments were carried out to determine the drop size 
distribution and nucleation site density. Microphotographs taken of 
the condensing process resulted in a measured distribution of drop 
Sizes from 10 to 3000 microns. 


A model of the dropwise condensation process has been proposed 
which correctly predicts the observed parametric trends. The problem 
was divided into analyzing the heat transfer through a single drop 
and then utilizing the drop distribution to predict the average heat 
transfer through the surface. 


The results of the experimentally determined heat transfer 
measurements have been compared to those predicted by the model. It 
was concluded that the drop conduction resistance represents the major 
limitation to drop condensation heat transfer. Interfacial effects 
account for 25 percent of the total resistance during low pressure 
(T = 88 °F) and 10 percent during atmospheric pressure condensation. 
For dropwise condensation on a vertically oriented, chemically promoted 
copper surface the sum of all other resistances is less than 15 percent. 


Thesis Supervisor: Peter Griffith 
Title: Professor of Mechanical Engineering 
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tT RODUCTION 


Compared to the boiler section of modern steam power plants, 
the condenser has received very little attention and consequently 
advanced ever so slightly in terms of performance and size reduction. 
Significant improvement in condenser design will come, however, when 
the drop-wise mode of condensation is utilized. It has been shown that 
heat Betaeer coefficients for dropwise condensation can be 30 to 40 
times greater than for filmwise condensation. Converting this to 
hardware could mean a reduction in condenser size by a factor of one 
third. 

There are three areas that need thorough investigation before 
dropwise condensation will appear attractive to condenser designers: 
1. The Non-Condensable Gas Problem. The slightest amount of non- 
condensable gas in the vapor causes a drastic drop in the condensing 
rate for dropwise condensation. Drop condensation is not significantly 
better than film condensation for contaminated steam. The location and 
design of vents to reduce the build up of non-condensables to a tolerable 
limit is an area which has not received enough attention to date. 
2. The Promoter Problem. In order to obtain dropwise condensation, a 
non-wetting surface must be produced by "promoting" the surface. At 
present, a long-lasting promoter which causes good quality drop con- 
densation is lacking. The ideal promoter must be relatively inexpensive 
and must not introduce a significant thermal resistance to heat transfer. 
Without such a promoter, dropwise condensation will never be reliable 
On) prachweal. 


3. The Limit Problem. The heat transfer coefficient for dropwise 





condensation decreases as the pressure and the vapor to surface 
temperature difference decrease. Designers must know the magnitude 
of this drop in heat transfer performance to determine whether drop 
condensation is beneficial under certain extreme operating conditions. 
Information on the mechanisms which limit heat transfer for dropwise 
condensation under all conditions would be very valuable. 

The purpose of this work is to investigate the limit problem and 
determine the factors which affect dropwise condensation. In this 
thesis, the author will first present a description of the dropwise 
condensation process. The purpose of this section will be to specify 
how heat is transfered in dropwise condensation and what factors affect 
this process. Following this section will be a documentation of the 
visual observations and heat transfer measurements which have been per- 
formed by past researchers. At this time a brief description of the 
cause of parametric effects will be supplied. 

The author carried out careful heat transfer measurements to 
supply further information on some of the more important parametric 
trends. The results of these measurements serve as a basis for compari- 
son for the later analytical work. 

One of the more important parts of this work involves the modeling 
of the heat transfer process of dropwise condensation. The proposed 
model correctly predicts all of the parametric trends observed to date. 

The primary experimental work involved the measurement of the drop 
distribution. This distribution is used in conjunction with the model 
to predict the heat transfer performance for dropwise condensation. 

From the experimental data and the heat transfer model, conclusions 


have been drawn concerning the limiting mechanisms of dropwise condensation. 
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ifs” DESCRIPTION OF THE DROPWISE CONDENSATION PROCESS 


Neaily forty years have elapsed since ene and his co-workers 
[1] first examined the heat transfer process for dropwise condensation. 
During this period, conflicting views have arisen concerning the 
mechanisms of dropwise condensation. The controversy:focuses on two 
important aspects of the process. The first question involves the 
origin of the droplets observed on the surface. The second point 
needing clarification is whether the heat is being transfered through 


the drops or through the "bare" area between the drops. 


2.3 The "Film Theory" 

The early investigators in the field of drop condensation believed 
that a thin film of liquid existed between the visible drops. As a 
group, they contended that heat is transfered both through the drops 
anid through this liquid film. Furthermore, they explained that this 
film was the source of the larger visible drops. 

In 1936, Jakob [2] proposed that vapor molecules condense on the 
bare surface between drops forming a film of liquid. This film grows 
until it reaches a critical thickness when it "rolls itself together 
and fractures" to form droplets. Eucken [3] also assumed that conden- 
sation takes place on a continuous film between the drops. However, 


"surface dif- 


he stated that by means of a mechanism which he called 
fusion", the condensed liquid migrates into the base of drops. Emmons 
[4] proposed a slightly different mechanism for heat transfer. He con- 
tended that vapor molecules condense in the region between the drops, 


lose some thermal energy, re-evaporate, and then condense again on the 


surface of the drop. 





ie 


These three viewpoints have two ideas in common: that vapor con- 
denses in the region between the drops and that the resulting liquid 
film is the source of the drops. A number of investigators [5-11] have 
utilized this general model in developing equations to predict the 


heat transfer in dropwise condensation. 


2.2 The "Nucleation Theory” 


More recent investigations tend to disprove the "film theory". 
Adherents of the "nucleation theory” believe that drops originate from 
discrete nucleation sites randomly distributed on the surface. They 
contend that there is no liquid film between drops and that the only 
path for heat transfer is through the drops. 

Tammann and Boehme [12] first suggested in 1935 that droplet for- 
mation is a result of a nucleation process. He observed that drops 
form at the same sites during successive condensation cycles. Several 
recent investigators [13 to 20] have substantiated his observations. 

Umur and Griffith [21] supplied a piece of evidence which strongly 
proves that dropwise condensation is a nucleation phenomenon. Based 
on thermodynamic considerations and an optical study designed to deter- 
mine film thickness, they showed that a film no greater than a monolayer 
thick can exist in the area between drops. This work proved that there 
is no liquid film between drops. 

Several investigators, McCormick and Baer [22], Le Fevre and Rose 
foi, Rose [24], Gose et al. [25], Hurst [26], Mikic [27] and Glicksman 
[28] have utilized the nucleation theory in their modeling of the drop- 


wise condensation process. 


2.3 Description of Process 





Dropwise condensation is one of the most complex heat transfer 
processes which exist. Strictly speaking, the process is neither 
steady nor uniform over the condensing surface. Furthermore, many 
different heat transfer mechanisms are involved in drop condensation. 

When a vapor comes in contact with a cold non-wetting surface, 
dropwise condensation results. Because of the non-wetting or hydro- 
phobic nature of the surface, the liquid condensate forms as droplets 
as opposed to a continuous film as in the case of film condensation. 

When the vapor first encounters the condensing surface, the vapor 
makes the phase transformation at discrete nucleation sites. These 
tiny microscopic drops grow by direct condensation onto their surface. 
As the drops grow and occupy more surface area, they bump into each 
other and coalesce. At this stage, drops grow both by condensation and 
coalescence. As a drop increases in size, its rate of growth by conden- 
Sation decreases due to conduction limitations. These drops now grow 
primarily due to smaller drops feeding them by means of the coalescence 
process. When a drop reaches a critical size, the surface tension is 
eecone by gravity and the drop departs. This is the end of one com- 
plete cycle. 

As drops coalesce and depart, bare surface area is exposed to the 
vapor. Primary drops form at the uncovered nucleation sites thus sup- 
plying a fresh batch of growing drops. In the case of a vertically 
oriented surface, the departing drops slide down the condenser surface 
sweeping the area in its path of all drops. This exposure of surface 
area by coalescence and departure is essential to the continuation of 


the process. 
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On a condensing surface one encounters a range of drop sizes from 
a fraction of a micron up to a few millimeters in diameter. Drops can 
be said to be in one of four stages of development: 1) nucleation, 

2) condensation and coalescence, 3) coalescence, and 4) departure. 
The rate of growth and the number of drops of a particular size depend 
on Which stage of development it belongs to. 

As a drop grows in size and passes from one stage to the next, its 
rate of growth decreases. Furthermore, due to the mechanism of coales- 
cence, the number of drops of a particular size decreases with increasing 
diameter. It is through the numerous small drops that the majority of 
the heat is transfered. 

Assuming the drops are hemispherical in shape, there are two fac- 
tors which affect heat transfer: 

1) the condensation rate on the "active" drops, and 


" drops 


2) the number and size distribution of the “active 
The "active" drops are those which are growing by condensation. These 
drops are in the nucleation and condensation-coalescence stage of develop- 
ment and transfer the majority of heat to the surface. In contrast, 
the "dead" drops, belonging to the coalescence and departure stages, 
transfer little heat. They essentially act like insulators occupying 
a certain amount of the condenser surface area. 

hor a given fluid-surface Best E on: the condensing or growth 
rate of the "active" drops of a given size is governed by two factors: 

1) the system pressure, and 


2) the "true" surface to vapor temperature difference. 


Altering the system pressure changes the thermodynamic properties of 
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the steam and liquid phases. These property changes affect the drop 
growth rate. The "true’' temperature difference may be defined as the 
difference between the flat interface saturation temperature of the 
vapor and the actual surface temperature under an active drop. The 
measured temperature difference is often much greater than the true 
temperature driving potential. 

The number of "active drops" is determined by the nucleating ability 
of the surface and factors which affect the coalescence and departure 
mechanisms. Increasing the nucleation site density increases the number 
of primary drops which are the source of the active drops. Mechanisms 
which cause "good quality" coalescence and early departure of the large 
drops increase the amount of area available for active drops. Thus, 
nucleation, coalescence, and drop departure affect the number of active 
a . on the surface. 

It is quite easy to explain all parametric trends for dropwise 
condensation heat transfer by considering just two factors: the size 
distribution and the growth rate of active drops. Theoretical equations 
have been derived by previous investigators and modified by this author 
to predict the drop growth rate. From these expressions one can predict 
what factors affect the condensation rate on the active drops. All 
these expressions for the growth rate include a temperature difference. 
It should be pointed out that the driving potential used in these equa- 
tions should be the "true" temperature difference. The drop growth 
equations will be derived in Chapter 8. At that time, further discussion 

The actual surface temperature under an active drop differs from 
the mean surface temperature because of the non-uniform heat flux 
through the condensing surface. This "constriction" effect will be 


described later in this thesis. Any temperature drop through the pro- 
moter layer will also lower the temperature driving force. 
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of the true temperature driving potential will be given. 

The second factor which affects dropwise condensation heat trans- 
ietecie size distribution of active drops, is more difficult to pre- 
dict. No theory exists to predict fully how parameters such as tempera- 
ture difference and surface chemistry determine the nucleation site 
density. Trends can be predicted, and one must rely on these to explain 
heat transfer performance. The amount of area vacated by the dead drops 
(and thus made available to active drops) is governed by coalescence 
and drop departure mechanisms. Again, these effects can only be pre- 
dicted qualitatively. 

The present understanding we have of the subject of dropwise con- 
densation is a result of three main lines of investigation by past 
researchers: 1) modeling of the drop condensation process, 2) visual 
observations, and 3) heat transfer measurements. Since so many of the 
previous models relied on the erroneous film theory of drop condensation, 
it is not worthwhile to make an extensive literature search in this 
area. The model proposed by this author in Chapter 8 will be presented 
with reference to the more recent model of Rose [23,24]. 

Visual observations of the dropwise condensation have been bene- 
ficial in describing many of the parametric trends. Areas of investi- 
gation have included determination of nucleation site density, drop 
population, and size and shape of coalescing and departing drops. 

These works will be reviewed in the next chapter. 

A wealth of heat transfer data exists illustrating the parametric 
effects on dropwise condensation. This data is useful in testing the 
accuracy of our heat transfer models and knowledge of the drop conden- 


sation process. 
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Upon completion of the review of visual observations and parametric 
effects on heat transfer, the work carried out in this investigation 


will be fully described. 
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III. VISUAL OBSERVATIONS OF THE DROPWISE CONDENSATION PROCESS 


eae, every experimental paper to be found on the subject of 
dropwise condensation describes some visual observations of the process. 
Frequently, the visual portion of the work is merely a verification 
that good quality drop condensation is taking place. However, some 
researchers have taken great pains to measure parameters such as nu- 
cleation site density, area covered by drops, drop population, departure 
sizes, and sweeping cycles which are very basic to the full understanding 
of the dropwise condensation process. 

It is unfortunate that a good deal of the visual data obtained by 
some of the earlier investigators cannot be directly applied to heat 
transfer measurements. The reason for this is that temperature and heat 
flux measurements carried out in conjunction with the visual observations 
are in error for most of these experiments. The effects of these heat 
measurement errors are discussed in detail in the next two chapters. 

As will be shown in Chapters 8 and 9, the correct determination 
of the drop distribution is essential for the accurate prediction of 
the heat transfer. Several investigators have counted and measured the 
number of drops on the condenser surface and presented their results in 
many different ways. Unfortunately, because of omittance of important 
pieces of information on their methods of taking and processing their 
raw data, the results of these investigations are of little use for 


the purpose of predicting drop condensation heat transfer. For this 


reason, an experimental program was carried out by this author to supply 


the necessary measurement of the drop distribution. 


The results of past visual investigations of the dropwise conden- 
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sation process will now be reviewed. The principal purpose in doing 
this is to supply information to help explain the basic mechanisms of 
the heat transfer process. In addition, the results of several of the 
works serve as a basis for comparison for the experimental work carried 


out by this investigator. 


3.1 Drop Population and Area Coverage 


Dropwise condensation heat transfer cannot be predicted unless 
the drop size distribution is known. To perform this calculation, an 
expression for the number of drops of diameter D to D + AD per unit 
area as a function of D is needed. Unfortunately, this expression is 
very difficult to predict analytically and equally difficult to deter- 
mine experimentally. 

Fatica and Katz [5] were the first to attempt such a measurement. 
Te results, reproduced in Fig. 1, give the number of droplets per 
unit area as a function of time after the beginning of a cycle. From 
the distribution data, Fatica and Katz computed that 45 percent of the 
area was covered by drops of diameter 100 microns or greater. Because 
of the inconvenient form, no further information can be drawn from their 
measurements. 

Hampson and Ozisik [30] performed a similar experiment and their 
results are reproduced as Fig. 2. Their data is also not presented in 
a form useful for computing heat transfer. Figure 2 shows that at the 
end of a cycle about 54 percent of the surface was covered by drops 
greater than 125 microns in diameter. 

The drop distribution measured by Sugawara and Katsuta [11] and 
Sugawara and Michiyoshi [29] is presented in Fig. 3. From this figure, 


the drop distribution N = 2 x 10/ p can be computed. However, the 
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authors do not mention the drop measurement band width, AD, so that the 
use of this data is somewhat restricted. The area covered by visible 
drops (greater than 300u in diameter) was calculated by the authors to 
be 54 percent. In addition, the authors figure that another 19 percent 
of the surface is covered by drops greater than 10 microns but less than 
300 microns. They thus concluded that 27 percent of the surface area’ 
was available for active drops. This figure agrees reasonably well with 
the results of this research. 

The distributions measured by Sugawara and Katsuta were obtained 
using surfaces with six different roughnesses. It is interesting to 
note that within this range of roughness, there is little variation in 
the drop distribution. However, not enough information was presented 
in their paper concerning roughness measurements or surface preparation 
techniques to draw any further conclusions. 

The most recent attempt at measuring the drop distribution was 
carried out by McCormick and Westwater [15]. Their results are repro- 
duced as Fig. 4. The most interesting fact that can be drawn from this 
work is that the drop distribution for drops smaller than 160 microns is 
dependent on the temperature difference. There is too little data for 
diameters greater than 160 microns to make such a judgement. The frac- 
tional area covered by drops greater than 2 microns was calculated to 
be .64 for the lowest drop population (lowest AT) and .72 for the highest 
population. 

The temperature differences measured by McCormick and Westwater 
are not the true temperature potential discussed in the previous chapter. 
Because they operated with a stagnant system one can rest assured that 


non-condensable. gases affected their experiment. The very low measured 
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heat transfer coefficient of less than 500 Btu/hr £t7°R substantiates 
this view. This author feels that the true temperature difference was 
so low in these experiments that the surface was not saturated with 
nucleation sites. This accounts for the dependence of the drop distri- 
bution on temperature difference as reported in this work. 

The results of these four investigations along with the data ob- 
tained in this work are compared in Fig. 5 by plotting the area covered 
by drops versus a nondimensionalized drop diameter. The experimental 


conditions for these works are summarized in Table l. 


3.2 Nucleation 

The nucleation experiments carried out by past researchers were 
involved in three principal areas: 

1) Proof that dropwise condensation is a nucleation phenomenon and 
identification of the nucleation sites; 

2) Measurement of the number of nucleation sites on the surface as a 
function of temperature difference; and 

3) Measurement of the critical temperature difference needed for nu- 
Seedeion and determination of the factors which affect the process. 

All of the ten papers which will be mentioned in this section 
prove that dropwise condensation is a nucleation phenomenon. Tammann 
and Boehme [12] were the first investigators to determine this fact. 
The series of papers by McCormick and Westwater et al. [13 to 16] drive 
home the nucleation concept in a very dramatic fashion. 

A great many facts about nucleation have been brought out in these 
papers. To describe each fact and observation in detail would take a 
great amount of space. In order to bring together in one paper the 


experimental observations in this area, the following list is supplied. 
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List of Experimental Observations Concerning the Nucleation Phenomenon 


SeEacement 


Associated with Dropwise Condensation 


Reference 


OP ee 
15,16,17, 


3,19 


ike) al 
iG. 17, 
18,19 


14 ,16,19 


13,14 


eee 


Observation 


Primary drops form at the same sites 
during successive condensation cycles 
(whenever surface area is exposed to the 
vapor due to coalescence or drop depar- 
ture). 

Imperfection in the surface such as 
scratches and pits can be identified as 
nucleation sites. 

Liquid debris left in scratches by large 
sweeping drops became nucleation sites 
for the next generation of drops. 
Artificial sites can be produced by 
making cavities in the surface with a 
small needle and by spark erosion. 
Scratches made by a scalpel also served 
as nucleation sites. 

Foreign particles placed on a condensing 
surface served as nucleation sites. The 
relative nucleating ability of particles 
can be determined by comparing their net 
heat of adsorption. Particles with the 
highest net heat of adsorption nucleate 


best. 





Statement 


10 


ita 


12 


13 


Reference 


14,15,19 


1k: 


105) 


ie 


eS 


Woed4 15419 


ives rs 
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yor 


Observation 


The number of nucleation sites on a 
surface increases as the surface to 
vapor temperature difference increases. 
The nucleation site density is greater 
for rough surfaces than for smooth sur- 
faces. 

Stainless steel surfaces exhibit greater 
nucleating ability than copper surfaces. 
The wettability of the surface is an im- 
portant parameter for nucleation. As 
the contact angle decreases, the critical 
temperature difference decreases. 

The nucleation site density increases 
with pressure for the same surface at 
the same temperature difference. 
Nucleation sites can be deactivated by 
drying the surface. The shape of the 
cavity is important in determining its 
Stability. 

The critical temperature difference 
needed for nucleation depends on “Secs 
surface energy of the substrate. 

The critical temperature difference 

for nucleation varies among sites. 
Pit-like sites are prefered to slip 


lines and Scratches. The size and 
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Statement Reference Observation 


shape of pits and scratches are also 
important in determining the critical 
temperature difference. 

14 ie Natural nucleation sites are randomly 
distributed on the surface. 

5 | 19 The presence of nearby large drops 
tends to suppress formation of primary 
drops at known nucleation sites. 

16 ie Oxide residue in surface blemishes 


enhances nucleation. 


The number of nucleation sites for a particular surface and tempera- 
ture difference have been measured by several investigators. Table 2 
Summarizes the experimental conditions and results for these experiments. 
The temperature differenceSreported with these densities are probably in 
error. None of the investigators making these microscopic studies of 
nucleation took the necessary care in removing non-condensable gases. 
The AT's so reported are much higher than the true temperature difference. 
The high nucleation site densities found in this work at the low AT of 
.5 °F (non-condensable gas free system) tend to prove this point. The 
nucleation results obtained during this investigation will be discussed 


in Chapter 7. 


3.3 Drop Growth Rate 


McCormick and Baer [13], McCormick and Westwater [15], and Dolloff 
and Metzger [31] measured the drop growth rate for dropwise condensation. 


In all three experiments, the diameter squared was found to vary linearly 
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with time for all but the smallest drops (see Fig. 6). This relation- 
ship for the drop growth agrees with theoretical equations derived by 
McCormick and Baer [13], Umur and Griffith [21] and Scriven [33]. 

In Chapter 8, the Bree raezeal prediction of the growth rate of 
drops will be discussed in detail. It will be shown that for large 
drops where conduction limits growth, the diameter squared is linear 
with time. For small drops other resistances dominate and this rela- 
tionship breaks down. 

McCormick and Westwater [15] measured the drop growth rate at two 
mmoesuress They estimated that increasing the pressure from 19 mm to 
39 mm of mercury increased the drop growth rate by a factor of about 2. 
Dolloff and Metzger [31] found a similar increase when the pressure was 
increased from 44.7 psia to 89.7 psia. The pressure effect on growth 
rate will be discussed from a theoretical standpoint in Chapters 8 and 9. 

The above authors also observed that the crowding together of drops 
on the surface slowed down their growth rate. From their movies they 
determined that a drop grows fastest when its neighbors were furthest 
away. McCormick explained this slow down in growth rate to be due to 
competition for vapor. This author believes that it is the "constric- 
tion'' resistance caused by non-uniformities in surface temperature 
that is the real cause. This effect will be discussed later in this 
paper. 

McCormick and Baer [13] noticed a sharp decrease in drop growth 
rate when non-condensable gases were introduced with the steam. The 
added diffusion resistance offered by the non-condensables accounts 


for this slower growth. 


3.4 Maximum Drop Size 
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The maximum size drop which can exist in equilibrium on a surface 
has been observed in detail in works by Fatica and Katz [5] and Sugawara 
and Michiyoshi [29]. In these papers a relationship between departure 
drop size and contact angle has been derived and then confirmed by 
experimental evidence. 

Sugawara and Michiyoshi measured the maximum drop size and the ad- 
vancing and receding contact angles for a brass surface promoted by 
oleic acid at several different surface inclinations. Their results 
are reproduced as Fig. 7. They found that the maximum drop diameter 
was independent of inclination but that the volume of the drop increased 
as the inclination approached horizontal down. Their predictions of 
maximum drop diameters from contact angle measurements, indicated by | 
the x points in Fig. 7/7, agree well with measured values. 

It should be noted that these measurements were made while the 
condensation was stopped. Therefore no dynamic effects such as drop 
vibration caused by coalescence are included. Sugawara obtained this 
data at room temperature (T = 20 °C). 

Fatica and Katz performed similar experiments in obtaining contact 
angle and maximum drop size information for six metal-promoter systems 
during atmospheric pressure condensation on a vertical surface. Ad- 
vancing contact angles ranged from 90 to 110 degrees and receding con- 
tact angles from 50 to 85 degrees. The diameter of departing drops varied 
from 1.6 mm to 3.8 mm. Their measured and calculated results agreed in 
a similar fashion as did those of Sugawara and Michiyoshi. 

A discussion of factors which affect departure drop size is pro- 
vided in Appendix B. The importance of departure mechanisms will be 


emphasized in that appendix. 
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3.5 Conditions Required for Dropwise Condensation 

In order to maintain dropwise condensation, nucleation sites must 
continually be exposed to the vapor. When two drops coalesce, surface 
area is exposed provided the liquid can recede over the surface. 
McCormick and Baer [13,34] have determined that the receding contact 
angle must be greater than zero for drop condensation to continue. 
According to MacDougall and Ockrent [35], the receding contact angle will 
be greater than zero provided that the normal contact angle measured on 
a horizontal surface is larger than 50 eecocne In terms of surface 
tension, McCormick and Baer [13] demonstrated that for surface tensions 
less than 46 dyn/cm, no area was vacated after coalescence and mixed 
condensation resulted. 

McCormick and Westwater [15] analyzed the coalescence problem and 
determined that the maximum area is vacated when coalescing drops are 
the same size. For equal size drops with 90 degree contact angles, a 
maximum of 21 percent of the previously wetted area will be vacated 
during coalescence. 

One of the interesting pieces of information obtained from the 
high-speed movies of Peterson and Westwater [16] is that a drop experi- 
ences as many as 400,000 coalescences going from nucleation to departure 
size. The large number of coalescences emphasizes the tremendous 
"activity'’ of the drop condensation process. The authors estimated this 
coalescence number for their ethylene glycol system (contact angle 72 


degrees and maximum drop diameter 1.6 mm). 


The author feels that the receding contact angle must be a good deal 
greater than zero and the normal contact angle greater than 50 degrees 
in order to have "good quality" drop condensation. Further research 
is needed to verify these numbers. 





3.6 Sweeping Drops 


For condensing surfaces in the vertical orientation, departing 
drops leave the surface by sweeping down the surface. Sugawara and 
Michiyoshi [29] determined that the sweeping cycle is reduced as the 
heat flux is increased and is shorter at the lower part of the con- 
densing surface. Increasing the heat flux results in a higher conden- 
sation rate and more departing drops. The reason why the lower portion 
of a condenser surface gets swept more often than the upper area is that 
a drop increases in size as it proceeds down the surface sweeping a 
larger Mpeceeetere it. Fatica and Katz {5] noticed these Same eitects 
but added the fact that the track width of sweeping drops becomes con-_ 
stant within 3 to 4 inches of the top of the surface. The increased 
condensate in the drop from then on elongates rather than widens the 
drop. 

Le Fevre and Rose [36] determined by heat measurements that the 
drop condensation heat transfer coefficient is independent of plate 
height in the range of 1 to 4 inches. This seems contradictory to the 
fact that the sweeping cycle decreases with height in this same range. 
The probable explanation is that area exposure due to Sonnet eee is 
more important than area exposure by sweeping at the heat fluxes in 
Ene above experiment. | 

Rose [24] observed that the fractional area covered by moving 
drops increases as the heat flux is increased from 60,000 to 600,000 
Btu/hr re His results are reproduced as Fig. 8. He concluded from 
this observation that the "blanketing" by falling drops should cause a 
reduction in condenser performance at very high heat fluxes. 


The visual observations documented in this chapter help clarify 
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the mechanisms of dropwise condensation. These results will be used 
in the next chapter to help explain the parametric effects on drop 


condensation heat transfer. 





, E362 
IV. PARAMETRIC EFFECTS ON DROPWISE CONDENSATION HEAT TRANSFER 


Experimenters have determined the effects of ten parameters on 
Etiemieat transter for drop condensation. In this chapter the author 
will carefully tabulate the results of experiments which best illustrate 
ene parametric trends. In addition, the cause of these effects will be 
discussed in the light of the previous description of the dropwise con- 
densation process. 


It is important to point out that wide discrepancies exist between 


the results of early investigators. As will be explained in detail in 


the next chapter on experimental techniques, the principal cause of 
these discrepancies is the presence of non-condensable gases in the 
steam. The slightest trace of gas will drastically reduce the heat 
transfer rate in dropwise condensation. When low coefficients are re- 
ported, the probable reason for the error is htentenesb ie gas effects. 
The reader should realize that the graphs presented in Figs. 9 
through 16 are intended primarily to illustrate trends. Frequently 
the levels 2s the heat transfer coefficient reported for a certain 
parametric study by different experimenters are not in agreement. This 
is because other effects (such as non-condensables) are also influencing 
the particular experiment. However, the parametric trends are usually 


quite consistent. 


4.1 Heat Flux 

For pure steam void of non-condensable gases condensing on a 
vertical copper plate, the heat transfer coefficient increases witi 
heat flux as shown in Fig. 9. At very low heat fluxes (low AT's), the 


coefficient rises sharply. For higher heat fluxes the coefficient has 
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been shown to remain nearly constant (curves 1, 3, and 6 of Fig. 9) 

or increase gradually with heat flux (curves 2, 4, and 5). It will 

be in that heat flux does affect the departure sizes of drops and 
that this could cause slightly increased steam side coefficients. 

One of the causes for the increasing coefficient with heat flux 
is the drop distribution. As the AT is increased from zero, the number 
of nucleation sites increases rapidly. The larger number of sites mean 
more ''active'' drops and increased heat transfer. The number of nuclea- 
tion sites soon reaches a limit which is then little affected by in- 
creasing the AT, and at this point nucleation ceases to be a controlling 
factor in the heat transfer process. 

As the heat flux is increased, however, the number and size of 
the large departure size drops are affected. Higher heat fluxes cause 
early departure of drops resulting in more area available for ele 
"active'’ drops. This in turn causes a higher heat transfer coefficient. 
This effect is expected to be small and barely measurable. Other parame- 
ters such as vapor velocity and surface roughness can also affect the 
departure size and the coefficient in a similar way. 

It should also be pointed out that the drop growth rate of the very 
small drops decreases with AT to a power greater than one. Since the 
majority of the heat is transfered through these small drops, then lowering 
the AT should decrease the heat transfer coefficient. The analysis de- 
veloped in Chapters 8 and 9 verifies this hypothesis. 

Figure 10 shows the results of all the experiments attempting to 
measure the variation of heat transfer coefficient with heat flux for 
a vertical copper surface. Curves 7 through 14 of Fig. 10 are for 


experiments where non-condensables have been a strong factor. The re- 
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maining curves (duplicated in Fig. 9) are, in the opinion of this 
author, 1Ost representative of drop condensation with “gas free' steam. 
Curves 4a, 4b and 10 of Fig. 10 are much higher than the others due to 
vapor velocity effects. Table 3 lists the experimental conditions re- 
ported in the references [5,7,8,30 and 36 to 46]. 

The earlier heat transfer data [5,7,8 and 40 to 45] for dropwise 
condensation is strongly affected by non-condensables. The low coef- 
ficients (below 30,000 Btu/hr £t?°R for heat fluxes greater than 50,000 
Btu/hr ae for atmospheric condensation of steam) are primarily due to 
non-condensable gas effects. Furthermore, many of these investigations 
have shown that the heat transfer coefficient decreases with heat flux 
after a certain heat flux is exceeded. Since the non-condensable gas 
effect increases with heat flux, this also indicates that this data has 
been governed by non-condensables. 

ati conclusion, one can say that the heat transfer coefficient in- 
creases rapidly at low heat fluxes due to nucleation and drop growth 
limitations. After this initial rise, the coefficient remains constant 
or gradually increases with the heat flux depending on how the system 
affects the departure size of drops. For steam contaminated by non- 


condensabies, low heat transfer coefficients can be expected. 


4.2 Pressure 

When describing the effects of pressure on the heat transfer 
coefficient, it is convenient to divide the data into below and above 
atmospheric pressure. As the sub-atmospheric pressure is increased, 
the heat transfer coefficient increases. This trend has been substan- 
tiated by the works of Gnam [41], Ma [47], Brown and Thomas [48], Tanner 


et al. [49] and the data obtained in this work (see Fig. 11). The in- 
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erease in coefficient with pressure is due partly to higher growth rates 
of the "active" drops caused by variation in thermodynamic properties 
and partly to the alteration in the size and numbers of "active" drops 
on the surface. The relative importance of these two factors will be 
discussed in Chapter 9. 

There is rather poor agreement between the two experiments carried 
out to measure the heat transfer coefficient for pressures greater than 
Menospheric (Fig. 12). Wenzel [38] reported a gradually decreasing 
coefficient as the pressure increased from 1 to 4 atmospheres. Con- 
trarily, in this region Dolloff and Metzger [31] found the coefficient 
to increase sharply. However, above a certain critical pressure, he 
observed that the coefficient decreases and then gradually increases 
again as the pressure is increased. 

It is difficult to say which of these pieces of data should be 
trusted. As will be pointed out in Chapter 9, increasing the pressure 
beyond atmospheric causes property changes resulting in slower drop. 
growth rates and a gradual decrease in coefficient. However, if the 
number of "active" drops continues to increase with pressure (the trend 
noticed for sub-atmospheric pressures), then this could cause an increase 
in coefficient. Depending on which effect is stronger, decrease in 
growth rate or increase in number of drops, the coefficient could in- 
crease or decrease with pressure above atmospheric. 

The experimental conditions described in the seven investigations 
into pressure effects are tabulated in Table 4. 

It can safely be said that experiments have shown that the heat 
transfer coefficient for dropwise condensation increases with pressure 


below atmospheric. This increase is due both to favorable drop growth 
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and drop distribution effects. Above atmospheric pressure there is 
discrepancy in the heat transfer data. Drop growth rates are expected 
to decrease with higher pressures but population of active drops to 
increase. The result could be a steady increase, steady decrease, or 
increase and then decrease in coefficient with increasing pressure 


above atmospheric. 


4.3 Non-Condensable Gases 

Undoubtedly, non-condensable gases can affect the performance of 
dropwise condensation more than all the other parametric changes com- 
bined. Even the smallest amounts of non-condensables allowed to accumu- 
late near the condensing surface offer such a high thermal resistance 
that it dominates over the small thermal resistance of drop condensation. 
Wenzel [38], Furman and Hampson [50], Hampson [51], Tanner et al. [39, 
49], Le Fevre and Rose [36] and others have shown that the steam side 
heat transfer drops rapidly as non-condensable gas content is increased. 

The velocity of the steam-gas mixture is actually more important 
than the concentration of gas. For a stagnant system the gases in the 
steam are carried towards the condensing surface by the steam and re- 
main in the vicinity of the surface causing a high mass transfer re- 
sistance. A large enough steam velocity past the surface, however, will 
sweep away the non-condensables reducing this resistance greatly. 

When comparing heat transfer data in the drop condensation litera- 
ture, one frequently finds low values for the heat transfer coefficient 
for experiments with supposedly "gas free steam'’. Le Fevre and Rose [36] 
have shown that no matter how long one deairates the steam, sufficient 


non-condensables remain to affect drop condensation. The six investiga- 
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tions in which the correct heat transfer coefficient for gas free steam 
have Deen measured have demonstrated the need for close venting [36,46] 
or reasonably high vapor velocities past the condensing surface [30, 37 
to 39] to minimize the resistance due to non-condensable gases. 
Non-condensables can also affect drop condensation by causing early 
breakdown of the promoter. The result is a "sloppy" drop condensation 
and a lower condensing rate. 
The conclusion that can be drawn about non-condensing gases is 
that they drastically reduce the dropwise condensation heat transfer 
coefficient. The reduction is due to the addition of a large diffusional 
resistance effectively lowering the temperature difference which seees 
the condensation. The effect of non-condensables increases with gas 


concentration but can be minimized by providing sufficient vapor velocity 


past the condensing surface. 


4.4 Vapor Velocity 


The most important role that vapor velocity plays in dropwise 
condensation is in the removal of non-condensable gases near the con- 
densing surface. However, it has been shown by Tanner et al. [39], 
Gaara et al. [52], and this author (see Fig. 13 and Table 5) that in- 
creasing the velocity in excess of the critical velocity needed to eliminate 
non-condensables also increases the heat transfer coefficient. This in- 
crease is due to early removal of departing drops. 

O'Bara, et al. [52] have shown that too high a velocity can decrease 
the heat transfer. They explained that high velocities (greater than 
5 ft/sec) cause a large shear stress which tends to flatten the drops. 
These flat drops will then occupy more surface area leaving less area 


available to the small "active" drops on the surface. 
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In conclusion, one can say that velocity increases the heat trans- 
Heme cOQerticrent for dropwise condensation. Starting from a stagnant 
system, the initial increase in velocity causes a rapid increase in 
condensation rate by removing non-condensing gases. Further increase 
in velocity above the value needed to remove non-condensables, increases 
the coefficient by early removal of departing drops. At very high vapor 
velocities, there could be an adverse effect caused by a flattening of 


the drops by the shearing action of the vapor. 


4.5 Conductivity of Condensing Surface 
Increasing the thermal conductivity of the condensing surface in- 
creases the heat transfer coefficient. Both Tanner et al. [53] and 
Griffith and Man Suk Lee [54] have shown that the steam side heat transfer 
coefficient is nearly 5 times greater for copper than for stainless steel. 
This trend is illustrated in Fig. 14 and Table oi Tanner and his col- 
leagues demonstrated that this effect is most pronounced at larger heat 
fluxes. The paper of Hampson and Ozisik [30] showed that the coefficient 
decreased in the order: Copper > Brass > Cr Ni > Monel > Staybrite. 
Mikic [28] has explained the effect of conductivity by describing 
a "constriction" resistance caused by the non-uniformity in surface 
temperature. He shows that lowering the conductivity of the surface 
increases this resistance and thus lowers the measured steam side coef- 
ficient. Lowering the thermal conductivity of the substrate decreases 
the "true" vapor to surface temperature difference. This in turn de- 
creases the growth of the "active" drops and lowers the heat transfer 
“The data of Griffith and Man Suk Lee is for a horizontal surface facing 
down. This is one reason why their data is considerably lower than 
fhat of Tanner. Also the author feels that there were errors in the 


measurement of the wall temperature in the former experiment. Neither 
of these effects alter the reported trend for conductivity and roughness. 
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Seemrictent. A further description of the constriction resistance 
will be presented in Chapter 8 and in Appendix G. 

It is interesting to note that the increased nucleating ability 
of stainless steel as compared to copper (reference [13]) is offset by 
the constriction resistance. That is, the large 'dead" drops affect 
heat transfer more than the increased number of "active'' drops for low 


conductivity surfaces. 


4.6 Surface Finish 

feanner eteal. [53]), Criffith and Man Suk Lee [54] and this author 
have found that heat transfer is higher for mirror smooth surfaces than 
for rough surfaces (Fig. 15 and Table 6). Surface finish can affect 
dropwise condensation by (1) altering the nucleation characteristics 
of the surface, and (2) changing the shape and size of the large drops. 

Mout surfaces have more nucleation sites than smooth surfaces. 

Thus rough surfaces should have a greater number of "active" drops and 
exhibit higher heat transfer coefficients. Since it has been demonstrated 
that the coefficient decreases with roughness, it can be concluded that 
nucleation effects are again offset by other considerations. 

Due to contact angle hysteresis caused by rough surfaces, the size 
and shape of drops are altered. It can be expected that more irregularly 
shaped drops will occur on rough surfaces decreasing the area exposed 
when two drops coalesce. In addition, the departing drops will be held 
back longer on rough surfaces. Both of these effects tend to decrease 


the amount of "active'’ area and consequently the heat transfer performance. 


4.7 Surface Inclination 





Surface inclination affects the dropwise condensation process by 
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varying the shape, volume, and sliding speed of falling drops. It will 
be shown in this paper that 90 percent of heat is transfered through 
the small drops (less than 150 microns). Therefore the inclinations 
that cause early departure of drops and rapid sliding speeds should 
give the best heat transfer. 

The investigators [30,32, and 38] who measured the effect of sur- 
face inclination on eee neces all found that the coefficient is 
greatest ate the plate is in the vertical position and decreases with 
departure from the vertical in ee: direction (see Fig. 16 and Table 
7). Rose [32] performed the most exacting experiments varying the in- 
clination in intervals of 10 degrees for five different heat fluxes. 

He observed that the coefficient has two maximums: one at 90° (vertical) 
and the other a 135° (45° facing down). These peaks in the coefficient 
versus inclination curves are most pronounced at the higher heat fluxes. 

In Appendix C, the topic of drop departure mechanisms and their 
effect on heat transfer will be discussed. The experimental results, 


that the coefficient is greatest for a vertical surface, is consistent 


with the fact that departure sizes are smallest for this orientation. 


4.8 Location on Condensing Plate 


Le Fevre and Rose [36] measured local steam-side coefficients at 
three locations on a vertical copper surface. They found that the local 
coefficient does not significantly depend on plate height from 1 to 4 
inches from the top. Shea and Krase [43] also measured coefficients 
at several height locations and found the coefficient to decrease with 
distance from the top surface. However, this decrease can be attributed 
to velocity and non-condensable effects. 


It is surprising that distance from the top of a vertical plate 
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aie sO eect et tect on the heat transfer coefficient. The frequency 
at which a ,. int on the surface is swept increases with distance down 
the plate. ‘ais means that the maximum drop size at points toward 
the bottom of the surface should be less than at the top of the surface. 
However, it is not certain whether the amount of area cleared by the 
sweeping action of falling drops is significant compared to that cleared 
by the coalescence mechanism at certain heat fluxes. 

Since Rose has shown that the heat transfer coefficient is not 
height dependent, then one can conclude that the number of "active" 


drops does not vary significantly with the location on the surface. 


4.9 Condensing Vapor 

Dropwise condensation has been observed for many vapors besides 
steam. Topper and Baer [55], using a Teflon-coated surface, maintained 
drop condensation with aniline, ethylene glycol, benzene, and nitroben- 
zene. Bobco and Gosman [56] observed dropwise condensation of iso- 
octane, decalin, ethanol, carbon disulphide, cyclohexane, isoheptene, 
and methyl alcohol on a surface promoted by fluorinated acids. Without 
any promoter, mercury has been observed to condense in the dropwise 
mode by Misra and Bonilla [57] and Ivanovski [58]. 

Due to the wide variation in apparatus used by the experimenters, 
it is fruitless to compare the heat transfer from vapor to vapor. Since 
the largest resistance for 'gas free'’ drop condensation is the conduction 
resistance through the drop, one can surmise that liquids with the 
largest thermal conductivity will exhibit the highest heat transfer. 
Of course, this statement is made assuming that the drop distribution 
and the quality of the drop condensation for the materials are comparable. 


Because of the large range of wetting ability of various fluids, this 
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assumption is probably not valid. 


4.10 * Promoter 

im order to obtain dropwise condensation, steam must come in 
contact with a cold hydrophobic (mon-wetting) surface. The non-wetting 
characteristic of the surface is brought about by "promoting" the sur- 
face. The best promoter is one which adds an insignificant thermal re- 
sistance due to its thickness, has a long life, and causes "good quality" 
dropwise condensation. Promoters can be categroized as: (1) chemically 
adhering, (2) physically adhering, (3) non-adhering, and (4) noble metals. 
A wealth of literature [59 to 68] exists describing studies of promotion 
techniques. These deal with the structure of the promoters, the length 
of their lives, and their applicability for industrial condensers. 

It is not advisable to measure the performance of promoters by 
comparing data of one experimenter with another. Other considerations 
such as non-condensable gases, velocity, orientation, etc. outweigh the 
promoter effect. However, Tanner et al. [53] and Le Fevre and Rose [36] 
each took heat transfer data with several promoters and thus the rela- 
tive merit of these promoters can be compared. Tanner reported that 
the heat transfer coefficient increased in the promoter series: Dibenzyl 
disulphide < Dioctadecyl disulphide < Montan wax and Montanic acid. Rose 
reports higher coefficients in the series: Dioctadecyl disulphide < 
"No. 1 Amine" (octadecylamine) < Di-S-octadecyl, 10 decanedixanthate < 
Dodecanetris (ethanethio) silane. In each investigation there was an 
increase in coefficient by about 50 percent in going from the first to 
the last promoter in the series. 

The promoter can affect the condensation by altering the contact 


angle of the drops and introducing a significant thermal resistance due 
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to its thickness. Lowering the contact angle decreases the amount of 
area exposed by coalescence and adversely affects the size and shape 

of drops. These in turn decrease the number of "active" drops on the 
surface. Introduction of an added promoter resistance lowers the "true" 
temperature difference and thus slows down the condensation rate. The 


results of both effects would be a lowering of the measured heat transfer 


coefficient. 
4.11 Summary 


The parametric effects on dropwise condensation heat transfer are 
summarized in Table 8. The primary cause for the performance of the 
heat transfer (whether due to growth rate or number of "active" drops) 


memeindicated for each case. 
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V. EXPERIMENTAL APPARATUS AND PROCEDURES 


The previous chapter has demonstrated that extreme care is needed 
in order to obtain accurate heat transfer measurements for dropwise con- 
dénsation. The series of papers by Rose [36 and 46] and Tanner [39,49, 
53] and their colleagues have carefully documented the important factors 
which affect accuracy in drop condensation experiments. The papers 
emphasize that non-condensable gases and errors in measuring the sur- 
face temperature account for the large discrepancies in heat transfer 
‘data encountered in the literature. 

Several months were wasted in this research because some of the 
recommended procedures were not closely observed. In order that future 
researchers do not encounter the same problems which this author did, a 
detailed description of the steps the author took in obtaining accurate 


heat transfer and photographic measurements will be provided. 


2.1 Description of Apparatus 
5.1.1 Design of Vapor System _and Chamber 


There were two principal considerations governing the design of 
tthe vapor system: elimination of the non-condensable gas problem and 
convenience for microscopic investigation. Luckily, the best designs 
for each consideration were compatible with each GieNcise 

The non-condensable gas problem was overcome by building a small, 
reasonably vacuum tight, "straight through" system. Vapor was generated 
in the boiler section in two 4 liter glass Erlenmeyer flasks at a maxi- 
mum: rate of one cubic foot per minute. The vapor was led to the con- 
densing chamber by By eee section of 3/4 inch ID rubber vacuum tubing 
and passed straight through" into the exit end of the system. At the 


highest condensation rates only one quarter to one third of the steam 
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was condensed. In all tests the vapor flow was great enough to sweep 
away any non-condensables from the condensing surface. 

The system could be operated at omoneureci Lone eToeaneme down to 
30 mm of mercury. For atmospheric pressure operation, heat was supplied 
to the boilers by two 1.5 KW hot plates and the vapor was allowed to 
pass through the exit valve into the atmosphere. Lower pressures were 
maintained by lowering the power of the hot plates and utilizing an 
aspirator pump. The pressure desired could be regulated by controlling 
the bleed valve in conjunction with the power input to the boiler. The 
pressure in the chamber was measured to .1 mm Hg with an absolute mer- 
cury manometer. 

A vapor flow rate of .25 cubic feet per minute (corresponding to 
an entrance vapor velocity into the chamber of 3 feet per second) was 
found to be sufficient to minimize non-condensable gas effects for both 
atmospheric and low pressure operations. The vapor flow rates were 
determined by passing the vapor produced in the boiler into a large 
condenser and measuring the condensate flow with time. A flow meter 
was also used to check these results. Agreement to within +10 percent 
was obtained. Figure 17 represents a schematic of the vapor system. 

The design of the condensing chamber is particularly important in 
insuring "gas free'’ results. The chamber should be small and free of 
potential dead air spaces. The chamber used, illustrated in Figs. 18 
and 22, was the second one designed, built, and tested. The original 
chamber was much larger (6 inches ID and 4 inches deep) than the second 
(2.65 inches ID and .5“inch deep). In the first chamber, because of 
its size, non-Gondensables were allowed to collect in stagnant pockets. 


At the time the large chamber was used, the boiler output was also 
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lower resulting in a significantly slower vapor velocity past the con- 
densing surface. Table 9 dramatically shows the increase in measured 
ee reient when the small "gas free'' system was used as opposed to the 
larger system. This again emphasizes the importance of taking great 
care in eliminating non-condensable gases in dropwise condensation 
experiments. 

While testing out the condensing chambers, the author noticed 
that when non-condensables were present, low frequency temperature 
fluctuations of up to 1 °F existed at the surface. Le Fevre and Rose 
[36] have explained that these fluctuations are due to:a mixing of the 
gasS-vapor mixture by the action of the falling drops. When the small 
chamber was used with high vapor flow rates, these temperature fluctu- 


famens at the surface ceased to exist. 


5.1.2 Design of Test Section 

Because of the high heat transfer coefficients (up to 50,000 
Btu/hr ft? °F) encountered in dropwise condensation, precise temperature 
measurement techniques are required. The papers by Le Fevre and Rose 
[36] and Citalogu and Rose [46] carefully outline some of the more im- 
portant considerations. 

Great care was carried out in this experiment to insure accurate 
temperature readings. Only high precision thermocouple wire, lead wire, 
thermocouple switches, and potentiometer were used. Seven thermocouples 
were placed along the axis of the .8 inch diameter, 2 inch long copper 
test section. In order to minimize errors due to conduction along the 
thermocouple wire, .005 inch diameter wire (copper-constantan) folded 


over twice was inserted in the .030 inch diameter holes (see Figs. 19 





ema 20 for details). In this way, the immersion depth was increased 
from -4 inch to over 1.5 inches. Appendix A describes the thermocouple 
mounting techniques and immersion tests in greater detail. In Table 9, 
the magnitude of the error due to incorrect temperature measurement 

is estimated. 

The copper test section was soldered to the stainless steel adaptor 
Sleeve which provided a tight O-ring seal in the back of the condensing 
chamber. Cork and fiber glass insulation wrapped around the copper 
section insured uniform axial conduction along its length. The con- 
densing surface protruded about .1 inch beyond the inside chamber wall 
into the vapor (see Fig. 18). 

The rear of the test section butted up against the cooling water 
chamber. Figure 21 represents a schematic of the cooling water loop. 
Ice water from a 60 gallon drum, warm water from a 30 gallon drum, or 


city water could be pumped through the cooling chamber at rates from l 
to 10 gallons per minute. The heat flux was controlled by varying the 
flow rate and temperature of the water flow. In addition, minor changes 


in heat flux could be controlled by varying the pressure of the cooling 
er oe an = amen emanaiamait aarti icra . 7 . ——— VT 7 
Ciambenwacainst the rear of the test section (thus varying the contact 
resistance between the two surfaces). 


2.1.3 Microscopic and Photographic Equipment 


One of the primary goals in the experimental program was to obtain 
detailed microphotographs of the dropwise condensation process. The 


purpose was to take picttres at several magnifications in order to 
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measure the drop size distribution. 


A specially constructed microscope mount was designed to enable 
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movement of the scope in the vertical and horizontal directions. This 
permitted exact location of the microscope's field of view on the surface. 
Long working distance objective lenses (made by Vickers) permitted high 
magnification (up to 400X) observation of the condensation process 
through the glass window. A Polaroid film pack adapted to a Leitz 
shutter mechanism was used to take the microphotographs. 

Two types of lighting were utilized in order to get the large range 
of pictures needed for complete photographic coverage of the drop size 
spectrum. Most of the pictures were taken using a standard vertical 
illuminator mounted behind the objective lens. When the 1/125 second 
shutter speed of the camera was sufficient to stop the action, an in- 
candescent bulb was used as the light source for the illuminator. For 
high magnification (200X and 400X) pictures, it was necessary to use a 
strobe lamp of duration 10 microseconds placed inside the vertical il- 
luminator to stop the action of the small "active" drops. When this 
technique was used, two vertical illuminators were mounted in series on 
the microscope. The illuminator with an incandescent bulb was first used 
to focus the microscope on the surface and then shut off. The strobe 
illuminator could then be blinked while the shutter of the camera was 


opened to take the picture. Figure 23 shows the complete test setup. 


5.2 Heat Measurement Procedures 


Peewee ceparation of the Surfaces 


Most of the heat transfer and photographic data were taken using a 
ee 
mirror smooth, chemically promoted, copper surface. The steps in pre- 
paring the surface for a typical run were as follows: 


1. The surface was polished to a mirror finish by rubbing it on 
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a metalurgical polishing cloth soaked with first .3 micron and then .05 
micron aluminum oxide polishing compound. All traces of oxide and promo- 
ter debris were removed by this process. | 

2. The polished surface was carefully washed with a cotton swab 
and distilled water. 

3. The specimen was further washed by rubbing the surface with 
a cotton swab soaked with carbon tetrachloride. At this time the surface 
was examined to insure that it was free of lint or other impurities. 
Following this washing, the surface was immersed face down in a basin of 
carbon tetrachloride for 10 minutes. 

4. The surface was then immersed face down in a basin of promoter 
solutions consisting of 1 percent by weight of dioctadecyl disulphide 
dissolved in carbon tetrachloride and allowed to stand 15 minutes. 

5. From the promoter bath, the test section was quickly placed in 
position in the condensing chamber, the cooling water turned on, vapor 
introduced, and condensation begun. 

The rough surface used in the tests was prepared by rubbing the 
mirror smooth surface in random directions on 500 grit aluminum oxide 
emery paper. The same cleaning and promoter Sequence used for the mirror 
surface followed after this roughening. 

It should be noted that on several occasions when a new test section 
was promoted for the first time, poor quality drop condensation resulted. 
After several polishing, cleaning, and promotion cycles were carried out, 
good quality drop peodens ation was obtained. The same event occurred 
when first Petecing to promote the newly roughened surface. The explana- 
tion must be that the chemical bond between the promoter and copper is 


not strong for first-time promoted surfaces. 
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The procedure used to produce the Teflon surface was recommended 
by engineers from the E. I. Du Pont De Nemours Company. Teflon 30 TFE 
fluorocarbon resin was diluted to a 5 to 10 percent solid dispersion. 

A .1 mil brushed chrome surface was first electrodeposited on the mirror 
smooth, copper pariiee,. The surface was then dipped in the TFE resin, 
withdrawn, and the excess liquid drained. A hot air gun was used to fully 
dry the surface. After drying, the test section was then sintered in a 
700 °F oven for 30 minutes. This completed the coating procedure. 

Before each condensing run, the Teflon surface was merely wiped 
clean with a cotton swab soaked with carbon tetrachloride. The same 
teflon surface was used for three condensing runs spread over a period 
of one month. During this period the quality of the condensation was as 


good as the chemically promoted surface. 


Meee neat Transfer Data Taking and Processing Procedures 


Citakoglu and Rose [46] have shown that promoting a copper surface 
Pith dioctadecyl disulphide in the above manner leaves an overabundance 
of promoter on the surface. This extra promoter thickness constitutes a 
measurable thermal resistance. As the washing action of the condensation 
process removes the excess promoter from the surface, the measured heat 
transfer coefficient increases. It is not until all but a monolayer of 
promoter is left that the coefficient levels out at its steady value. 

The procedure used in this experiment was to permit rapid conden- 
sation to take place for at least 3 hours before reportable heat transfer 
feta was taken. In the next chapter, this transient effect is further 


described. Later in the project, it was found that the transient time 





% 
The chrome prevented oxidation of the copper during the high 
temperature sintering step. 
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(time before excess promoter was washed off and the coefficient remained 


steady), could be shortened by spraying the surface with carbon tetrachlo- 
ride just after removing it from the promoter bath and before placing it 
into the chamber. 

A typical "data point" was taken by measuring the vapor tempera- 
ture with two thermocouples located inside the condenser chamber and the 
temperature at seven locations along the axis of the test section. The 
average of the two vapor thermocouple readings were taken as the vapor 
Semperature. Excellent linear gradients (see Fig. A-4) were found to 
exist along the test section and these were converted to heat flux by 
the Fourier rate equation. The surface temperature was found by extrapo~ 
lating the temperature gradient to the surface. Knowing the heat flux 


'- 


and the vapor and surface temperatures, the heat transfer coefficient can 
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be computed. 


The selection of the "best" straight line temperature gradient 
through the temperature data was done using least square techniques with 
a computer. This procedure removes any bias the experimenter might be 
inclined to add to the data. Appendix I lists the computer program 
written to process the data, and several examples of the linear tempera- 


ture gradient measured through the test section. 


a.3 Drop Distribution Measurement Procedures 


So sleeeionegraphic Procedure 


The mirror smooth chemically promoted surface was used for the 
runs in which microphotographs were taken of the dropwise condensation 


process. The surface was prepared in the same way as for a heat transfer 


* 
The two vapor thermocouples consistently agreed with each other to 


within .05 °F. 
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run. Photographs were not taken until after three to four hours of 


condensation had taken place to insure that data was taken while on the 


flat portion of the transient curve. While photographs were taken, heat 
flux-AT measurements were made to determine the operating conditions. 


Once steady condensation was achieved, the photographic work was 
begun. In order to obtain data over the complete range of measurable 
drop sizes (diameters from 10 to 3000 microns), photographs were taken 
at six meeHieicattons from 5X to 400X. Photographs were taken randomly 
PY respect to time so that there would be no agreement between the 
frequency of the pictures and the frequency of the coalescence and sweeping 
cycles. 

The ideal way to measure the drop distribution is to take a picture 
of a large section of condenser surface. By measuring and counting every 
drop from nucleation to departure size on the surface and dividing by 
the surface area, an accurate spatial average could be obtained. Un- 
fortunately, due to limitations of optics, this is not poe oben Instead, 
pictures must be taken at several different magnifications. Because the 
field of view at the high magnifications is so small, a large number of 
pictures must be taken to insure that a representative sample is obtained. 

By counting the number of drops of a particular diameter on each 
picture and then taking the average over the number of pictures, a repre- 
sentative count can be obtained for that drop size. This was the proce- 
dure used in this experiment. The microscope was focused at a particular 


kk 
point on the condenser ‘surface and a sequence of pictures taken in a 


* 
To get sufficient optical resolution to see the smallest drops (order of 
.l micron) would require magnifications of over 400X. The field of view 
at this magnification is so small that the large drops could not be seen. 


k 
The point of focus was chosen arbitrarily near the center of the surface. 
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random fashion with respect to time. In this way a time average of the 
number of drops of a particular size appearing within the field of view 
of the microscope was obtained. 

The time average can be shown to be equal to the space average by 
relying on the ergotic hypothesis. The ergotic hypothesis states that 
for random processes, the two averages are equal. It has been shown by 
McCormick and Baer [13] that nucleation sites are randomly distributed 
on the surface. Thus, it can be assumed that any size drop in the drop 
eeerum can exist on any point of the surface at any time. 

A question arises concerning how many pictures should be taken at 
each magnification. In order to obtain approximately the same confidence 
level for the drop counting data at each magnification, it is necessary 
to increase the sample size with magnification. This is because at high 
magnifications the field of view is small with respect to the large size 
drops. For this reason, 150 pictures were taken at 400X, 100 at 200X, 
55 at 120X, 40 at 60X, 24 at 10X and 16 at 5X. Appendix J discusses the 


picture taking procedures in greater detail. 


Beee Drop Counting Procedure 


After the sample of photographs were taken, the next job was to 
count and measure the drops. At first an attempt was made to count and 
measure every drop which could be resolved in the pictures. It soon be- 
came evident, after counting several hundred drops on each photograph, 
that this procedure was far too time consuming. A drop was measured by 
comparing its size with thé spacing between two parallel lines scribed 
on a clear plastic sheet. 

The procedure decided upon was to choose ten to twenty drop sizes 


ranging from 10 microns to departure size. Next, by examining the photo- 
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graphs, it was decided at which magnifications a particular size drop 
could be most accurately counted and measured. All drops within a cer- 
tain band width on either side of the mean diameter were counted and con- 
sidered to be of that size. Thus, on any one photograph all drops of 
size (D-AD) to (DAD) were counted and recorded to be of size D. Fre- 
quently the same size drop was measured and counted at two different 
magnifications to insure agreement from one set of photos to a oenCan 

The entire sample of photographs were first processed using a drop 
measurement band width of +10 percent. Because of the lack of clarity of 
the periphery of the drops at higher magnifications, it was decided that 
this band width was too small. It was just too difficult to decide for 
sure whether or not a drop should be considered to be in that small band 
width. For this reason the pictures were again analyzed using a band 
width of +20 percent. The author feels that the results obtained from 


this counting are more accurate. 


Seo imitations of Microscopic Equipment 


The smallest drops which could be accurately measured and counted 
on the surface were 10 microns. Smaller drops down to about 1 micron in 
diameter could be seen on the microphotographs but not measured and counted 
with sufficient accuracy. Greater microscopic resolution will be required 
in order to obtain data for drops less than 10 microns. 

It is difficult to get improved resolution because of two experi- 
mental problems. First, higher power objective lenses needed to give 
greater resolution have working distances so small that the large drops 
bump into the lense. Second, it is extremely difficult to stop the 


motion of the small "active" drops at high magnifications. These short- 





comings must be overcome to obtain drop distribution data down to 


nucleation size drops. 
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VI. RESULTS OF HEAT TRANSFER MEASUREMENTS 


fe General 

The parametric effects of heat flux, pressure, surface roughness, 
vapor velocity, Teflon promotion, and time have been investigated in this 
work. These parameters were chosen because they are applicable in answer- 
ing the basic question of what are the limiting mechanisms for dropwise 
condensation. The data presented here is for steady state, '"non- 
condensable ee good quality dropwise condensation of steam on a 
vertical copper surface promoted by dioctadecyl disulphide and Teflon. 

In the paper of Le Fevre and Rose [36], the authors discuss the ad- 
vantages and disadvantages of certain ways of presenting drop condensation 
heat transfer data. They pointed out that the traditional plot of heat 
transfer coefficient versus heat flux can be misleading when errors in 
temperature difference are appreciable. For this reason graphs of surface 
to vapor temperature difference versus heat flux are used to present the 
basic heat transfer data. Where appropriate, the heat transfer coefficient 
is plotted against the parameter being varied (such as pressure or velocity). 
In addition, plots of heat transfer coefficient versus temperature dif- 
ference are provided to illustrate trends used for comparison purposes in 
later chapters. 

For all cases in which heat flux-AT data was taken, the results 


were well represented by a linear relation: 


AT = mQ/At+b . (6.1) 


Ke 
The term "non-condensable free" is used to signify the absence of err 
to non-condensin ases. As discussed in the previous chapter 


extensive precautions were taken to insure that non-condensing gases were 
not allowed to accumulate near the surface. 
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The best straight line was chosen through the data by least squares 
technique with a computer (see Appendix I for listing of program). The 


heat transfer coefficient-AT relationship can be determined from equation 


G:. Lae: 


uo eon 
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meeseequation predicts that the coefficient is zero for AT <b. Thus, 


the critical temperature difference, AT . 


, needed to drive the drop 
E1Ie 


condensation process is b. The coefficient increases with AT until, for 
large AT's, the heat transfer coefficient is nearly constant and equal 
to 1/m. 

The author feels that the principal contribution made by the heat 
transfer data obtained during this experiment is the accurate coefficients 
measured at very low surface to vapor temperature differences. Up to now 
there has been little reliable data taken of AT's below 1 °F. It is hoped 
that these experiments will clarify the temperature difference effect on 


the heat transfer coefficient for dropwise condensation. 


6.2 Results for a Mirror Smooth, Chemically Promoted Surface 


Figures 24 and 27 present "gas free", dropwise condensation heat 
transfer data for a vertically oriented, mirror smooth, chemically pro- 
moted surface at atmospheric and low pressure (34 mm Hg). The data pre- 
sented in these figures was taken over a period of several months. 
Greater scatter pees in the data taken from day to day as compared to 
the data taken during one run. Over 90 percent of all the data was well 
within the +10 percent scatter band on AT-q/A plots. The straight line 
plotted through the data is the one chosen by least squares technique 
by computer. 


Figure 27 illustrates the performance of heat transfer coefficient 
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as a function of temperature diprérence.” The constant value of the 
coefficient for atmospheric pressure (38,200 Btu/hr £t°°F) is nearly 
Paree times greater than that for low pressure (13,500 Btu/hr £t-°F). 
In both cases, the pectetedent starts to drop sharply at about .5 °F. 
At a AT of .15 °F, the coefficients are 50 percent of their constant 
eimies. Ihe decrease in coefficient at the low AT's is due partly to a 
slowing down of the growth rates and partly to a decrease in the number 
of "active'’ drops. The model derived in Chapter 8 and discussed in 


Chapter 9 predicts this trend. 


6.3 Results for a Rough, Chemically Promoted Surface 


The rough surface (prepared by rubbing on 500 grit emery paper) 
showed a similar trend with AT as the mirror smooth surface (see Figs. 

Moe 26 and 29). The constant value coefficient (for AT > 1°F) for the 
rough surface was 30 percent below that for the mirror surface at atmos- 
pheric pressure and 20 percent lower for low pressure condensation. This 
drop in heat transfer performance with roughness is attributed to the 
adverse effects on the shape and size of large drops. 

An important observation can be made about the rough surface operating 
at low AT's. For a AT below .2°F, the rough surface exhibits higher heat 
transfer rates than the mirror smooth surface. The relatively high coef- 
ficient at low AT for the rough surface is due to its greater nucleation 
characteristics. At very low AT's, nucleation theory predicts that a 
rough surface will have a greater number of nucleation sites. 

From these results it can be concluded that a mirror smooth surface 
gives better performance than a rough surface at normal operating conditions. 


* 
The curve plotted through the h versus AT data was chosen by computer 
moon Eqs. (6.1) and (6.2). 
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However, when operating at extremely low AT's, a rough surface may be 
desirable due to its greater nucleating ability. 


eee Results for a Teflon Coated Surface 


Figures 26, 28, and 29 show the rather large decrease in heat 


transfer performance for the Teflon surface compared to the chemically 
= ——— LST. esas 





promoted surfaces. The reason for this is that the Teflon coating intro- 
duces a significant thermal resistance due to the conduction loss across 
feemickuess., The resistance of the Teflon layer can be calculated 
knowing its thickness and thermal conductivity and added in series with 
the condensation resistance. 

The conduction heat transfer coefficient for Teflon 30 TFE-Fluoro- 
carbon resin of thickness 1.5 microns (5 x 10> inches) is 20,000 Btu/ | 
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°F. The overall coefficient (promoter in series with condensation 
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Besistance) of the Teflon surface is about one third that of the mirror 
smooth, chemically promoted surface at atmospheric pressure and a little 


over half at low pressure. In Appendix H, the overall heat transfer coef- 


ficient is predicted as a function of the thickness of the Teflon layer. 


6.5 Results of Pressure Tests 

Heat transfer data was taken with the vertically oriented, mirror 
smooth, chemically promoted surface at 30 different pressures ranging 
from 24 mm Hg up to 760 mm Hg. A constant AT of about 3 °F was maintained 
while the pressure was varied. As ieTeR ee eaees in Fig. 30, the heat 
transfer coefficient decreases linearly with saturation temperature over 
this range. Figure 11 shows the same data plotted against saturation 
pressure and compared to the data of previous investigators. 


x 
Du Pont, the manufacturer of Teflon 30, predicts a layer thickness of 
1.5 microns when following their suggested application techniques. 





i ~64- 


The following linear relation was found to exist between the heat 
transfer coefficient and the saturation temperature in the range of T 
Ss 


ae 
metween // F to 212 °F for a AT of 3 °F: 


h = 220 (T. =) : Cor) 


at 


In this relation, h is measured in units of Btu/hr £t-°F and ee in E, 

Equation (6.3) probably specifies the relation between heat transfer 
coefficient and saturation temperature over a sizeable range of AT's. 
However, at low AT's where nucleation becomes important, this relation 
would change. Also, at very high heat fluxes other mechanisms would domi- 
nate and the equation would not be valid. 

Chapter 9 will explain that the decrease in heat transfer coefficient 
with pressure is due to two factors: the decrease in drop growth rate due 


to property changes, and the alteration in the drop distribution resulting 


in fewer active drops. 


6.6 Results of Velocity Tests 


During early tests, a low capacity boiler was used resulting in an 
insufficient vapor velocity past the condensing surface. The heat transfer 
coefficient (see Table 9) was found to be low for these low velocity tests 
due to non-condensable gas effects. 

With the higher capacity boiler, a vapor flow rate of 1 cubic foot 
per minute could be achieved resulting in a maximum vapor velocity at 
the entrance of the chamber of about 12 feet per second. It was found 
that a velocity of 3 feet per second at the chamber entrance was suf-~ 


+ 
ficient to eliminate all non-condensable effects. WEthmechis, velocity 


x 
In all tests the vapor blew across the condensing surface. 
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the downward track of the falling drops was nearly vertical. 

Heat measurements were made at four velocities to determine the 
Perect Of excess velocity on drop condensation on a vertical surface. 

The sweeping angle which the falling drops made with the vertical was 
also measured. The variation of heat transfer coefficient and sweeping 
angle with vapor velocity is shown in Fig. 31. For velocities less than 
3 feet per second, surface temperature fluctuations indicative of non- 
condensable gas effects became noticeable. Because of these temperature 
fluctuations, heat transfer-velocity measurements were not made in this 
region. The dotted line represents only the author's guess at the shape 
of the oe low vapor velocities. 

The initial rapid increase in coefficient with velocity (illustrated 
by the dotted line of Fig. 31) is due to the removal of non-condensing gases 
from the vicinity of the condensing surface. The gradual rise in heat 
transfer for velocities greater than 3 feet per second is a result of early 


departure of the large drops due to the shearing action of the vapor. 


6.7 Transient Effects 

Surfaces promoted with dioctadecyl-disulphide (and probably other 
chemical promoters of the same type) are strongly time dependent in the 
early stages of their use. During the first hour or two of condensation, 
the coefficient first rises and then drops sharply (see Fig. 32). 

The rise in coefficient is due to the removal of excess promoter by 
the washing action of the drop condensation process. As the thickness of 
the promoter layer is decreased, the conduction loss through the promoter 
is lowered and the overall steam side coefficient (promoter plus conden- 
sation resistance) increases. It is thought that all but a monolayer of 


promoter remains chemically bonded to the surface. 
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The sharp drop-off in coefficient is surmised to be due to an altera- 
tion in the wetting characteristics of the surface. This could affect 
both the quality of the condensation and the nucleating ability of the 
Surirace. 

After about three to four hours of rapid condensation, the coef- 
ficient levels off and remains nearly constant with time. The results 
of runs made on five different days are illustrated in Fig. 32. One 
notices that there is a slight difference in the steady state value of 
coefficient from day to day. The scatter of +10 percent about the average 
of 38,200 Btu/hr to °F is the maximum encountered throughout the experi- 
ments. The data indicated by dots is for the rough surface. 

Thirteen hours was the longest duration of any run made during this 
project. One should note,however, that the coefficient will begin to de- 
crease gradually after an extended period of condensation. This is due 
to a breakdown in promoter effectiveness in maintaining a non-wetting 
surface. As the quality of the condensation becomes poorer, the coef- 
ficient will decrease. 

Citakoglu and Rose [46] also reported a transient effect during their 
experiments with a surface promoted with dioctadecyl-disulphide. They 
noticed that the coefficient increased sharply during the initial stages 
of condensation similar to this experiment. However, they noticed no peak 
and subsequent decrease in coefficient. Their coefficient merely evened 
out after the initial increase. The disagreement between Rose's and this 
author's transient curves is probably due to different surface-promoter 


characteristics. 
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* 
6.8 tie, at eG) Vibieation of Surface 


Som, preliminary heat transfer measurements were made while vibrating 
the test ection in and out of the condensing acres The vibration was 
dene by scupling the test section and cooling water chamber assembly to a 
"Shak or" manufactured by Goodman Industries. The purpose of the test was 
to show that vibrating the surface would shake off the departure size 
drops prematurely and result in an increase in heat transfer. 

Unfortunately, only a limited amount of data was taken during these 
vibration tests. The shaker burned out after part of the data was taken 
for the first heat flux to be investigated. The results are summarized 
toedable li. 

The slight increase in coefficient with vibration is due to early 
remeval of departing drops. This effect is discussed in Appendix B. The 
data indicates that high vibrational frequencies are more effective to 
this end. 

It should be pointed out that vibration should show a greater effect 
for horizontally oriented and low conductivity surfaces. For these sur- 
feces, the large departing drops play a significant role in lowering the 
coefficient. Verification of this trend would be beneficial in understanding 


the importance of the departure region of the drop distribution spectrum 


omeneat transfer. 


6.9 Conclusions 


The author feels confident that the heat transfer measurements reported 





“The design and construction of the "shaker" used in the vibration tests 
were done by Mr. Howard Charney as a term project for a Heat Transfer 
Course at M.I.T. This author acted as one of the supervisors for the 
project. 
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gn this chapter correctly demonstrate the parametric trends of heat flux, 
pressure, surface roughness, vapor velocity, and Teflon promotion. The 
two experimental errors, non-condensables and surface temperature measure- 
ment, which have plagued past researchers, have been minimized in this 


work. The results are summarized in Table 10. 
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VII. PHOTOGRAPHIC MEASUREMENTS 


wo, Drop Distribution 


The principal goal of the visual experiments was to measure the 
drop distribution for dropwise condensation on a vertical copper surface. 
A secondary goal was to determine the nucleation site density. 

The shortcomings of the previous attempts to obtain these measure- 
ments have been mentioned in previous chapters. Accurate heat transfer 
data was not taken in conjunction with visual data. In addition, insuf- 
ficient background information limited the usefulness of previous drop 
distribution data. 

The resulting drop distributions obtained in this experiment for 
system pressures of 34 mm Hg and 760 mm lg and a AT of .5 °F are plotted 
in Fig. 33 and listed in Tables 12 and 13. These two distributions were 
measured using a band width, AD, of +20 percent. 

Due to the resolution limitations of the optics, the smallest drop 
which could be accurately measured and counted was 10 microns for low 

x 
pressure and 20 microns for atmospheric pressure condensation. Drops 
down to size of one to two microns in diameter could be seen on the 
microphotographs but could not be measured and counted with sufficient ac- 
curacy to include in the distribution data. Thus, the drop distribution 
below 10 microns could not be obtained. 

Since the distribution at the low end of the drop spectrum could 


not be measured, "educated" guesses were required to continue the distri- 


x 
The more rapid growth rate at atmospheric as opposed to low pressure 


condensation at the same AT results in less clarity in the microphoto- 
graphs. In addition, the much larger number of small drops at high 
pressure made counting more difficult. 
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bution down to the smallest drop size. Six distributions, illustrated 
in Figs. 34 and 35, were chosen spanning the possible choices between 
two extreme cases. The lowest distributions (L-0 for low pressure and 
A-O for atmospheric pressure condensation) approximate the case where no 
coalescence takes place between nucleation size drops and drops of 10 
microns. These distributions also assume that the drop'‘growth rate is 
constant in this size range. The highest distributions (L-6 and A-6) 
assume that the coalescence and growth rates for drops within this range 
are the same as for a drop of 10 microns. The real drop distribution 
must lie between these two extreme distributions. 

From the measured and predicted distributions, the area covered by 
drops could be computed. The results for three distributions at each 
pressure are plotted in Figs. 37 and 38. 

Because drop distribution data can be interpreted in many ways, the 
raw counting data is supplied in Appendix J. In addition, a step by step 
description of the data taking processing procedures are listed. The 
taking of the photographs and counting of the drops involved over six 
months of painstaking labor. The purpose for presenting all the raw data 
is to make it available for other possible interpretations and processing 
schemes. 

In Chapter 2, it was stated that there are four possible stages of 
development for a drop: the nucleation stage, the condensation-coalescence 
Stage, the coalescence stage, and the departure stage. The growth rate 
and coalescence rate varies for a drop as it progresses from one stage to 
Breenext. This would imply that there should be four separate distributions 
over the complete range of drops from nucleation to departure size. 


When a "primary" drop nucleates from a site, its growth is purely due 
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to condensation. Because of the very small size of the nucleating drop 
(order of .1 to 1 micron in diameter), neighboring drops appear to be 
distant from each other. Thus there is little probability that these 
primary drops will bump into each other and coalesce. Since there is virtu- 
ally no coalescence among the primary drops, their number remains constant. 
The result of this reasoning will be a drop distribution which remains 
constant in this nucleation region as shown in Figs. 34 and 35. 

As the drops continue to grow in size, the probability that coales- 
cence will occur increases. Thus, the number of drops of a particular 
size will decrease. The point at which coalescence becomes important is 
unknown at this time. The author arbitrarily chose one micron as the 
separation point of the nucleation region with the condensation-coalescence 
region. Figures 34 and 35 illustrate the change in distribution from A 
constant value in the nucleation region to a decreasing value in the con- 
densation-coalescence region. 

The distributions within the nucleation and condensation-coalescence 
regions are dependent on the nucleation site density and the growth rate 
of drops. In addition, the amount of surface area which the drops in the 
coalescence and departure regions occupy has a strong effect on the number 
of "active" drops which can exist on the surface. Thus it can be concluded 
that the distribution at the low end of the drop spectrum will vary with 
experimental conditions. 

When a drop reaches a size somewhere between 10 and 150 microns in 
diameter (depending on the pressure), the conduction resistance becomes 


so large that very little condensation takes place on its surface. At 


“The drop distribution will be constant in this region only if the growth 
rate is reasonably independent of size. Also, it is assumed that the 
sweeping actionSof the large drops do not affect these primary drops. 
For first approximations, these assumptions are valid. 
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this stage the growth of the drop is governed primarily by coalescence. 
Glicksman and Rose [27] have analyzed the coalescence and packing process 
in an attempt to predict what the distribution should be in this region. 
In their analysis they ano that the distribution would be independent 
of all variables which affect the condensation rate. They conclude that 
the distribution for good quality drop condensation within this region is 
universal for all condensing surfaces operating at normal operating con- 
ditions (where nucleation is not a limitation). 

Glicksman's predicted results in the form of area coverage as a 
function of nondimensionalized drop size are compared to measured values 
in Fig. 5. The agreement with this author's measured distributions for 
large drop sizes is quite close. 

Departure mechanisms affect the distribution of the very large drops 
on the surface. The number and size of the departing drops will vary with 
the orientation of the surface. Any of the parameters which affect contact 
angle such as surface roughness and the surface tension of the liquid also 
alter the size of departing drops. 

The author observed that the condensation rate affects the drop dis- 
tribution in the departure region. The effect of heat flux is illustrated 
in the photographs of Figs. 39 and 40. For the same surface at the same 
operating pressure, it is evident that departing drops are larger for low 
heat fluxes. As the heat flux increases, the rate at which the "active" 
drops bump into the large drops increases. This coalescing action causes 
the large drops to vibrate. Due to the greater vibration rate at higher 
fluxes, departure-size drops are more unstable and leave prematurely. A 
further discussion of the departure mechanism, including the heat flux 


effect, is included in Appendix B. 
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The photographs also show heat flux has an effect on the velocity 
of falliing drops. At high heat fluxes, the departing drops accelerate 
rapidly and never hesitate during their fall. (Note that the shutter speed 
of 1/125 sec could not stop the action of the falling drops.) But at low 
heat fluxes, the falling drops are sluggish and frequently stop momentarily 
during their downward journey. Again the explanation is that the rapid 
coalescence rate lowers the restraining force on the large drops. 

Figure 36 compares the drop distribution for atmospheric and low 
pressure condensation at the same AT and for the same surface. It is 
evident that there are fewer large drops on the surface during atmospheric 
as compared to low pressure condensation. The author believes that the 
reason for this is a combination of the higher heat flux and lower surface 
tension for the condensation at the higher pressure. Appendix B outlines 
some of the parameters which can affect the distribution in the departure 
region. The analysis is by no means complete and serves only to introduce 
an area of study of dropwise condensation which has been chiefly neglected 
up to now. 

The point separating the coalescence and departure regions seems to 
be somewhere between 500 and 1000 microns. This conclusion is reached by 
observing where the low pressure and atmospheric pressure distributions 
cross in Fig. 36. 

The principal purpose in measuring the drop distribution was to 
Emaoie one to predict the heat esoeHehES: for dropwise condensation. Un- 
fortunately, the complete distribution could not be obtained due to 
experimental difficulties in obtaining data for drops smaller than 10 
microns in diameter. Figures 34 and 35 illustrate seven distributions 


for each operating pressure which span the possible limits. The problem 
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now was to choose one of these distributions for each pressure which 
approximates the true distribution. 


"correct" distribution was to 


The procedure used to choose the 
determine which one made the predicted heat flux equal the measured 
heat flux. In the next two chapters this procedure will be described in 
detail. The distributions so chosen for each of the operating pressures 
are compared in Fig. 36 and Table 14. These are the distributions which 


the author feels correctly approximate the true distributions for the 


experimental conditions used in this investigation. 


7.2 Nucleation 

Originally, the primary goal of this thesis was the determination 
of the nucleation site density for dropwise condensation as a function 
of surface chemistry and surface to vapor temperature difference. Al- 
though the emphasis was shifted to the measurement of the drop distribu- 
tion, some observations can be made concerning nucleation. 

Figure 41 compares the density of small drops condensing at the 
same pressure Se = 88 °F) at a AT of .12 °F and .37 °F. The photo- 
graphs certainly demonstrate that the nucleation site density has increased 
in this AT range. The effect of pressure on nucleation is illustrated by 
the photographs of Figs. 42 and 43. The nucleation site density is con- 
siderably greater at the same AT for the atmospheric pressure compared 
to the low pressure conditions. 

The theory of nucleation related to dropwise condensation is out- 
lined in Appendix D. The effects of AT and pressure are discussed in 
some detail in this appendix. Briefly, it can be said that the surface 
to vapor temperature difference is the potential which drives the phase 


transformation over the nucleation energy barrier. More nucleation sites 


Base 


are activated at higher AT's. The primary effect of pressure is the 
lowering of this energy barrier due to property changes. The lower 
contact angles at the higher saturation temperatures also add to this 
decrease in the energy barrier. 

The main body of photographs taken for the drop distribution study 
were for a temperature difference of .5 °F. This AT was chosen because 
the relatively low resulting condensation rates enabled easier examina- 
tion of the process. The large sample of high magnification photos ob- 
tained at this AT for low and atmospheric pressures enabled the author 
to estimate the nucleation site density. At a pressure of 34 mm Hg and a 
fieot -> F, a site density of 200 million (2 x 10°) per square centimeter 
was estimated. For atmospheric pressure at the same AT of .5 °F, the 
density was approximately two to three times greater than at the low 
pressure. 

The nucleation site density was estimated by counting the smallest 
drops (about 1 to 2 microns) which could be resolved on the microphoto- 
graphs. When the location of primary drops coincided on several photos 
then that location was considered a nucleation site. 

The reason why only approximate values are stated is that it was 
extremely difficult to identify the great number of primary drops on 
the surface at this AT. The primary drops are so small and so close 
together at the magnifications used (200X and 400X), that resolution 
limitations prevented any greater accuracy. As the AT is lowered, nu- 
cleation studies become eee. because the density decreases (see Fig. 41). 
However, detailed counting of the nucleation site density was not carried 


out for temperature differences lower than .5 °F. 


7.3. Discussion of Photographic Study 
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The photographic data obtained in this work extends our knowledge 
of the drop distribution and nucleation site density. Accurate heat 
measurements made in conjunction with the photographic work enable one 


to use the data in predicting heat transfer performance for dropwise 


condensation. 
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VIIL. MODEL OF THE HEAT TRANSFER PROCESS FOR DROPWISE CONDENSATION 


The model of the heat transfer process for dropwise condensation 
proposed in this section is similar to that developed by Rose [23,24]. 
The problem is analyzed by first determining the heat transfer through 
a single drop and then calculating the heat flux through the surface 
from knowledge of the drop distribution. Rose's analysis successfully 
predicts heat transfer performance when compared to experimental results. 
However, four constants had to be "fit" into his equations to achieve 
this agreement. 

The basic concept of Rose's model has been continued in this analysis. 
However, because of more precise knowledge of the drop distribution, a 
bit "cleaner" presentation can be afforded. Also, two additional re- 
sistances to drop condensation heat transfer are considered in the pre- 
sent model to take into account the promoter and surface conductivity 


effects. 


8.1 Heat Transfer Through a Single Drop 


8.1.1 Resistance Concept 


The basic assumption for this model for the heat transfer through a 
single drop is that the resistances due to the various heat transfer 
mechanisms are independent and additive. Although this is not completely 
true, the originsof the resistances are sufficiently remote from one 
another that the assumption is valid for first approximations. 

Throughout this analysis drops will be assumed to be hemispherical 
in shape. Contact angle measurements very close to 90° for most of the 
promoted surfaces substantiate this assumption. Internal circulation 


within a drop due to thermocapillarity will be neglected. Lorenz and 
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Mikic [69] have shown that circulation increases the heat transfer through 
large drops by less than 5 percent. In addition, it will be assumed that 
the vibratory motion of the drop due to coalescence (discussed in 
Appendix B) does not Meer Cie heat transfer through a drop. 

Heat transfer will be assumed to take place only through drops. 
The present school of thought,that no liquid film exists between the 
drops and that no heat is transfered through the "bare" surface area,will 
be accepted. Overwhelming experimental evidence supports this view. 

The resistances to heat transfer for dropwise condensation are: 
1. Diffusion Resistance, Re? through non-condensable gases concentrated 
near the surface. The temperature drop, AT? due to this resistance can 
be estimated by considering the relationship among the partial pressures 
of the gas and vapor and the saturation pressure. However, little de- 
tailed analysis has been carried out to predict this diffusion resistance. 
The importance of non-condensables has been dramatically demonstrated in 
the discussion of past heat transfer measurements (see Fig. 10). 
2. Curvature Resistance, Ro. Because of the curved interface, the 
saturation temperature of the vapor in equilibrium with a drop is less 
than the saturation temperature of vapor in equilibrium with a liquid 
at a flat interface at the same pressure. Thus for a drop, there will 


be a temperature difference (T. for a curved interface - T. for a 


at at 


flat interface) which must be exceeded in order for it to exist on 

a condensing surface. This AT. can be thought of in terms of an equiva- 
Tent resistance in order A cote it with other effects. However, 

the "curvature" resistance is not a resistance (AT/Q/A) in the true sense 
of the word since the temperature difference is independent of heat flux. 


3. Interfacial Mass Transfer Resistance, Ri. A pressure difference is 
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required to drive the mass transfer across the vapor-liquid interface. 
This pressure difference can be converted to a temperature difference 
and a resistance calculated. 

4, Drop Conduction Resistance, Race The resistance to conduction heat 
transfer through a drop is the governing resistance for the large drops. 
This resistance is computed considering no internal circulation or 

drop motion. 

5. Promoter Resistance, SK There is a temperature drop through the 
promoter layer due to conduction loses. Since the thickness of most 
chemical promoters are extremely thin, this resistance can often be 
neglected. However, for Teflon it has been shown that the promoter re- 
sistance must be Pena i dexed. 

6. Constriction Resistance, Rom Because the majority of heat is 
transfered through the small "active" drops, the heat flux is not uni- 
form over the condenser surface. The crowding of the heat flow lines 
results in a temperature drop in the substrate which is greater than 
that due to straight conduction. This "constriction" resistance, analyzed 
by Mikic [28], has been shown to be significant for low conductivity 


condensing surfaces. 


The six resistances and the resulting temperature drops are il- 
lustrated in Fig. 44. These resistances are discussed in greater length 
in Appendices C and E through H. . 

For most drop condensation heat transfer experiments, the average 
heat transfer coefficient is measured. The usual technique is to obtain 
the temperature gradient through the bulk material on which the condensa- 
tion is occuring and convert this to heat flux using the known thermal 


conductivity. The average surface temperature is determined by extrapo- 
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lation of the straight line temperature gradient to the surface. The 
average heat transfer coefficient is then defined as 


= O/A 
h (| -T) CS) 

S W 
where T. is the saturation temperature at a flat interface. The local 
heat transfer-eeefficient—on_the surface is the coefficient for a single 
drop 


Q/A 


T -T CoyZ) 
o W 


In this equation Q is the heat transfered through a single drop and A 


1s the base area of a drop. Defined in this way, the local coefficient is 


= | 


Sail datas vcs yell 
Cc 1 


+R +R «= Ses) 
ne Pp cm 


dc 


The next step in the analysis entails evaluation of the six re- 
Sistances. For this work, the diffusion resistance due to non-condensable 
gases will be neglected. In addition, the vapor will be assumed to be 
saturated. The justification for this is that the experiments were 
carried out with "gas free", saturated steam. 

For the computation of the heat transfer through a single drop, the 
"true" temperature driving potential is (T. - ae where T. is the satura- 
tion temperature at a flat interface and i is the average temperature at 
Euembase Of the drop. The temperature drop due to constriction resistance 
(T. “ T Om and due to the promoter layer (T on - a can be neglected 
when discussing the heat transfer through a single avoon Then their 
effects can be added in series with the average resistance of all the 
ee [28] has shown that the constriction effect due to the crowding 

of heat flow through the periphery of a drop where the conduction path 


is shortest is of little importance compared to the crowding due to 
EiiemdistGibution of "dead" drops on the surface. 


en 
—- a 
a7 i —=_> 
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drops when computing the overall heat transfer coefficient (see Appendix 


G and H). Thus for this analysis, the total temperature difference to 
be considered is 


AT. = Fx +. AT. a AT. 7 AT 46 a i te Pe : (8.4) 


AT. = AT . = 5 AT. + AT Ac : (8.5) 


That is 


Sel.2 Lempéerature Drop Due to Curvature Resistance 


In order for a drop to be in mechanical and thermodynamic equilib- 
rium on a surface, the temperature of the drop must be colder than the 
plane interface saturation temperature corresponding to the system 
pressure. The derivation for the temperature drop through the interface 
due to curvature is outlined in detail in Appendix C. The resulting 
equation ie 


2 eo 


ae (8.6) 
. eee 


When the AT due to curvature is equal to the total AT. available 
fenat is T. - Ey then the minimum size drop which can exist on the 
surface is 9 ie 5 f 

r = ——— a ; (8.7) 


wee Pe aoe 


min 


Substituting (8.7) in (8.6), a simple expression for the tempera- 


ture difference to curvature can be found to be 


e.g 
AT. = a Ate (8.8) 


8.1.3 Temperature Drop Due to Interfacial Resistance 


The temperature drop due to the interfacial resistance is 


hceee eee (8.9) 
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where Q is the heat transfered through a single drop of radius r, and 

h, is the interfacial heat transfer coefficient. The coef fiietent, hh. 
a 

can be evaluated from the equation derived by Nabavian and Bromley [70] 


2 
L2H 
20 M 
eee) ee 
1 2 a” ‘2TRT T Vv 
S S 2g 








h (8.10) 
The term, &, is called the condensation coefficient and can be thought 
of as the fraction of vapor molecules striking the liquid interface which 
enter the liquid phase. 
Appendix E outlines the derivation of equation (8.10). Also in- 
cluded in this section is a discussion of the condensation coefficient. 
Throughout this paper it will be assumed that a = 1. The justifications 


for making this assumption are discussed in the appendix. 


8.1.4 Temperature Drop Due to Conduction Resistance 


The conduction resistance for a hemispherical drop of radius, r, 
can be approximated by comparing it to the resistance for a right- 
circular cylinder of radius, r, and height, r. The resistance for the 


cylinder is 








R = (8.11) 
cyl TE. 
and for the drop 
te 
Rae 7 (8.12) 
ae are 


where C is a shape factor which must be less than one. Mikic [28] has 
shown that if C = 1/4, this approximation for R,, agrees well with an 
exact solution of the problem derived by Umur and Griffith [21]. Other 
investigators [13,24] agree that the shape factor is near le 4 

The derivation of the drop conduction resistance is given in 


greater detail in Appendix F. The resulting temperature drop due to 
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conduction through a drop assuming a shape factor of 1/4 is 


ee a : (8.13) 
Tr 


oe oeeeeat. [ransfer Through a Drop 


The transfer of sensible heat from the vapor to the drop is negli- 
gible compared to the latent heat released as a result of condensation. 
For this case 


7 av 
OM a (8.14) 


where the volume, V, for a hemispherical drop is equal to 2/3 me Thus 


_ 2 dr 


8.1.6 Total Temperature Drop 


The total temperature drop due to curvature, interfacial, and drop 
conduction resistances can be found by substituting equations (8.8), 


foe) and (8.13) into (8.5): 


Te. 33 
AT. = ae ae to - BE (8.16) 
E h, 20 GkTtr 


Now if the expression for the heat transfered through a drop (8.15) 


is substituted in (8.16), then 


is pH OH 
_ min fg dr Le. weedus 
AT. =. AT. =f hy aca 7a Tar C3017) 


8.1.7 Drop Growth Rate 


Solving for <= and substituting D = 2r, one gets an expression 


for the drop growth rate 


fe Dep 
dD 4 min 
ee een Lites , 8.18 
dt oH ay D/2k + 2/h, ( ) 


fg 





a 


Equation (8.18) can be integrated to find the relationship between 


diemeter, D, and time, t. Assuming that at t = 0, D=D_., then 
| min 


2 2 


2k ! ( aay D in’? D in? a aE 1n(D De. 
eons p>) + D. In (@-pD.)) 

hy min min min 

ie 2 2 4 


If the curvature term is neglected, then (8.19) simplifies to 


ie? 2 2 4 
aoe OL) + he (D-D,) = oil, AT. tt (8.20) 


When the interfacial term is neglected, the relationship becomes 


1 Z Zz 
a b-5 CD i. D in? Ba 2D Vin (D ~ Dia) Bs ie In (D - ee 
1 Z is 4 
ioe D vin = oH AT. f : (3220) 


When only the conduction term is kept 


one? ae A 
rao, =o 7) = ou, AT. ts. (8.22) 


Throughout this derivation, the drop growth rate has been derived 
for an "isolated" drop. In reality, drops are growing in the presence 
of their neighboring drops. The constriction resistance takes into ac- 
count the nonuniform surface temperature due to the crowding of the heat 
flow through the active drops. This resistance has the effect of lowering 
the "true" AT as compared to the measured AT (see Appendix H). Instead 


of substituting the measured AT into (8.16), the "true' AT driving the 
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condensation rate should be used. However, this will not be the proce- 
dure carried out in this derivation. Instead, the computed value of the 
constriction resistance will be added in series with the resistance of 
all the drops averaged over the surface area. In this way the total 


temperature difference (AT = Te = T can be used in the above expressions. 


8.2 Heat Transfer Through the Surface 


If there is one drop of diameter, D, on a unit area of condenser 


surface, then the heat transfer per unit area, A, would be 


dV 1 2 dD 
(Q/A), = Pe, aE = PH ee (8.23) 
Sm@estituting (8.18) into (8.23) 
2 Lae nae 
(Q/A), = 1D AT. D/2k + 27h, : (8.24) 


Now if there are N dD drops of diameter D to D + dD per square centimeter 


* 
of condenser surface , then the heat transfer through that area would be 


d(Q/A), = NdD (Q/A), 
or 
— 2. 1 Din’? 
d(Q/A) . = NdD TD AT. D/ak + 2/h. , (8.25) 
if 


On a condenser surface there are drops ranging in size from D in 
(the minimum drop size) to De (the maximum drop size). If the drop. 


distribution, NdD, is continuous over this range, then the average 


heat transfer per square centimeter of condenser surface is 





D 

% ° = ° 

The total number of drops of size D, to D, is N tp N dD. According 
to this definition N = AN/AD where AN is the 1 number of drops 
counted per square centimeter in the size range D to D = AD. 
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D 
max 
(a), =] aaa), 
Die 
min 
or pmax 
— 2 i Din’? 
(Q/A) | = TAT. N D D/2k + 27h, dD (8.26) 
ere 


The heat flux for dropwise condensation can be predicted using 
(8.26) provided that the minimum and maximum drop sizes and the drep 
distribution are known. The maximum drop size was determined experi- 
mentally from the low magnification photographs. For low pressure 
condensation, D ax was found to be 3 mm whereas for atmospheric pressure 
it was 2.5 mm. These numbers were used as the upper limits in (8.26). 

Due to microscopic limitations, the drop distribution could not be 
determined experimentally over the complete range of drop sizes. There- 
fore it was necessary to use a trial and error approach to determine 
the most "correct" distribution at the low end of the spectrum. The 
details of this procedure will be discussed in the next chapter. 

The minimum size drop which exists on a condenser surface is the 
smallest drop which has just nucleated. For heterogeneous nucleation 
from a uniform temperature vapor on a plane substrate, the size of a 
nucleating drop is the same as the minimum drop size specified by Conv). 
No smaller drop than this can exist on the surface. 

Because there is a temperature gradient in the vapor near the wall, 
one can argue that the nucleating drop will be as much as twice as large 
as the minimum drop size. Therefore, the actual size of a nucleating 
drop is somewhere between D vin and 2 D vin’ As will be shown in the next 
chapter, a negligible amount of heat is transfered through drops smaller 


than one micron. Thus it is important whether one uses Das or 2D en as 
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the lower limit for the integration of (8.7). For the sake of simplicity 
the smallest drop which can exist on the surface will be considered to 


be D ce 
min 


To summarize the results of the analysis, the following equations 


are restated: 


—_ Din’? 


D/2k + 27h, dD 


i 
— 


C3527) 


20 M 


6 - e) “aR 


measured maximum drop size 


measured and proposed drop distribution 
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IX. PREDICTION OF HEAT TRANSFER FOR DROPWISE CONDENSATION 


9.1 Determination of Drop Growth Rate 


In the previous chapter, the growth rate due to condensation for 
an isolated hemispherical droplet was derived by considering the curva- 
ture, interfacial mass transfer, and drop conduction resistances. Equa- 
tion (8.18) states the relationship between growth rate, and diameter, 
D. 

Figure 45 compares the growth rate distributions (2 vs. D) for low 
and atmospheric pressure conditions. From this plot, one sees that the 
drop growth rate increases with pressure for all size drops. The maximum 
growth rate, occuring at .2 microns for atmospheric pressure, is nearly 
eight times greater than the maximum rate, occuring at .4 microns, for 
low pressure condensation. It is the interfacial resistance which pri- 
marily accounts for this decrease in the growth rate with pressure. 

The individual effects of the three resistances on drop growth rate 
can be seen by dropping off terms in Eq. (8.18). The curvature effect 
can be neglected by setting De = 0. In reality, this is saying that 
all size drops are large compared to Din’ Letting the interfacial heat 
transfer coefficient, hj» go to infinity results in the neglecting of 
interfacial effects. When both of these simplifications are made 
simultaneously, then only the drop conduction term is left. 

This procedure of dropping off terms was carried out and the results 
plotted in Figs. 46 and 47. At both low and atmospheric pressures, 
neglecting the curvature term doubles the growth rate for the smallest 
size drops. However, by the time the drops have grown to a diameter of 


5 microns, curvature effects have disappeared completely. 
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In the case of low pressure condensation, the maximum drop growth 
rate jncreases from 280 u/sec to 4000 u/sec when the interfacial re- 
sistance is neglected. Whereas for atmospheric conditions, the maximum 
growth rate increases rot 2000 u/sec to 4400 u/sec. This indicates 
that if it were not for the mass transfer effects, the growth rate for 
small size drops would be nearly independent of pressure. Only drops 
smaller than 40 microns in diameter are affected by the interfacial re- 
sistance. 

For drops much larger than the minimum size, the rate of growth is 
very nearly linear with AT. However, since the D of Eq. (8.18) varies 
inversely with AT (see Eq. (8.7)),the growth rate does not exhibit this 
linear relation for small drops. Figures 48 and 49 illustrate the 
effect of varying AT on the drop growth rate distribution. Since the 
expression for the average heat flux through the surface (Eq. (8.26)) 
involves the integral of the growth rate distribution, these figures 
show that the heat transfer rate will decrease sharply at low AT's. 

The effect of AT on the maximum drop growth rate is described in 
Fig. 50. This plot shows that the growth rate decreases rapidly for 
AT's less than 1 °F. This trend forecasts a decreasing heat transfer 
coefficient in this low AT region. This observation will be discussed 
again in Section yep or this chapter. 

The relationships between diameter and time, as derived in equations 
(8.19) through (8.22), are presented in Figs. 51 to 54. The first graph 
demonstrates how the interfacial resistance dominates at low pressures 
for small drop sizes. The difference between curves 1 and 3 of this 
figure illustrates this point. Figure 52 shows that the interface re- 


sistance loses its dominating role as the pressure is increased. (The 





=9= 


@erference between curves 1 and 3 of Fig. 52 is slight as compared to 
Mee, 51.) 

Figure 51 indicates that when the interfacial mass transfer re- 
sistance plays a dominant role, the diameter varies linearly with time. 
It appears that this will be the case only for small size drops growing 
in a low pressure system. 

Figure 53 continues the display of diameter versus time for larger 
diameters. In this range, the mass transfer effects are beginning to 
become unimportant. 

The drop continues to grow rapidly until conduction limitations 
take over. Figure 54 shows that p° is linear with time for drops greater 
than about 40 microns. This growth rate law has been experimentally 
verified by McCormick and Baer [13] (see Chapter 3). 

The growth rate process is but half the picture in dropwise conden- 
sation. The number and distribution of drops on the surface are equally 
important. The growth rate effects are "weighted" because there are 
more small drops on the surface. This concept will be explained in the 


next section. 


9.2 Selection of "Correct" Drop Distribution 


Since the drop distribution was not measured down to the minimum 
meep size, the heat transfer cannot fully be predicted utilizing the 
model outlined in the previous chapter. Figures 34 and 35 graphically 
illustrate seven distributions which span the range of possible distri- 
butions for the low end of the drop spectrum. Each of these distributions 
was substituted in Eq. (8.27) and the resulting heat flux obtained by 


integrating over the complete range of drop sizes. Numerical integration 
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was required for this job due to the complexity of the integrand. A 
list of the thermodynamic properties and the computer program written 
to perform this calculation are listed in eee in 

The "correct" distribution is the one which results in a predicted 
heat flux closest to the measured heat flux for that pressure. Figures 
55 and 56 illustrate the results of these calculations. From these cal- 
culations, distribution A-3 was chosen to be the "correct" one for atmos- 
pheric distribution and L-3 for low pressure distribution. 

Chapter 7 included a discussion of the expected shape of the drop 
distribution curve. It was noted at that time that the distribution is 
constant in the nucleation region and begins to decrease in the condensa- 
tion-coalescence region because of the mechanics of the drop growth and 
coalescence process. The nucleation drop size, the level of the constant 
feeribution, and the demarcation point between the two regions have aii 
been guessed at and then determined by "forcing" the predicted heat flux 
to equal the measured value. Other shaped distributions could have been 
chosen between the minimum drop size and the smallest measurable drop 
which would have predicted the heat flux. The author feels that the 
shape of the distribution chosen in this work is in agreement with the 
basic mechanisms governing the process in this region. The two correct 
distributions, L-3 for low pressure and A-3 for atmospheric pressure con- 


ditions, were used for all subsequent calculations. 


9.3 Relative Importance of Resistances 


The relative importance of the three resistances included in this 
analysis can be determined by neglecting terms in Eq. (8.27) before 
performing the integration. This procedure of dropping off terms was 


carried out and the results illustrated in Figs. 57 and 58. Curve 1 in 
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each graph represents the heat transfered through drops greater than 

diameter, D, taking into account the three resistances: curvature, inter- 

facial, and drop conduction. Curve 2 predicts the heat transfer neglecting 

the curvature resistance, and curve 3 the heat transfer neglecting the 

interfacial resistance. Curve 4 represents the heat transfer for the 

most simplified case of considering only the drop conduction resistance. 
The magnitude of each of the three resistances can be computed from 


data of Figs. 57 and 58 by the equation 


AT 


O/A : Co 1) 


Ko 


Using the heat flux predicted by curve 1, the total resistance (drop 
conduction, interface, and curvature) can be calculated. When the heat 
Flux from curve 2 is used, the resistance found from Eq. (9.1) is for 
drop conduction plus interface. Similarly, curve 3 gives the drop 
conduction plus curvature resistance, and curve 4 gives the drop con- 
duction resistance alone. By subtracting the resistance found from 
curve 4 from that of curves 2 and 3, the interface and curvature re- 
sistances can be computed respectively. 

The magnitudes of the resistances so calculated are displayed as 
bar graphs A and D of Fig. 59. The total resistance for low pressure 
condensation is about 3 times greater than for atmospheric conditions. 
The relative magnitudes of the individual resistances (i.e., comparing 
the drop conduction resistance for low and atmospheric pressures) are 
deceptive. Because the drop distribution is different for the two pres- 
sures, the values of the resistances for a single drop are "weighted". 
The drop conduction resistance averaged over the surface is over twice 
as large for low pressure as compared to atmospheric pressure conditions 


mainly because there are more large "dead" drops at the lower pressure. 
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The true effect of property variation (as opposed to drop distribu- 
tion plus property variation) can be realized by comparing bar graphs A 
with B and C with D. The low pressure distribution, L-3, was kept 
constant and the pressure varied between the two extremes Ce of 88 °F 
and 212 °F) and the results plotted as graphs A and B. The drop conduction 
term decreases by about 10 percent in going from low to high pressure be- 
cause of the drop in the thermal conductivity of water. The interfacial 
resistance is nearly 7 times greater for the low pressure conditions. 
This is purely due to property changes and not due to alteration of the 
value of the condensation coefficient (kept at 1 for both cases). The 
magnitude of the constriction resistance is about the same for the two 
pressures and results in less than 3 percent of the total resistance in 
either case. 

A similar calculation was performed holding the atmospheric distri-~ 
bution, A-3, constant and varying the properties. Graphs C and D of 
Fig. 59 show similar trends as mentioned for the low pressure distribution. 

The effect of the drop distribution on the relative magnitudes of 
the resistances can be obtained by comparing the bar graphs A with C 
and B with D. For graphs A and C, the thermal properties at T.., = 88 °F 
were held constant and the distribution varied. The drop conduction re~ 
sistance drops by a factor of one half in going from the low to the at- 
mospheric distribution. As stated before, this is because there is more 
"dead" surface area due to large drops for the lower pressure. The other 
resistances are little affected by distribution changes. 

The conclusion which can be drawn from Fig. 59 is that pressure 
affects dropwise condensation by varying both the drop distribution 


and drop growth rates. In general, when considering the magnitude of 
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the resistances averaged over the surface, both drop distribution and 
growth rate must be taken into account. 

Prememet. 59, it is evident that the concueeten resistance repre- 
sents the major portion of the total resistance for dropwise condensation. 
Therefore, for the same drop distribution fluids with higher thermal 
conductivities should exhibit greater heat transfer. 

The heat transfer for dropwise condensation of mercury was calcu- 
lated using the drop distribution data measured for water. The results 
are presented in Fig. 60. As expected, the total resistance for mercury 
is less than that for water. The conduction resistance is approximately 
10 times less for mercury compared to water. However, the interfacial 
resistance is much larger for mercury. These trends are purely due to 
the property differences between the two materials. 

It should be noted that due to the large variation in wetting 
properties for mercury and water, the drop distributions for the two 
materials will be greatly different. There will be far fewer nucleation 
sites for mercury err thus a smaller number of active drops. For this 
reason, it is doubtful whether drop condensation for mercury will even 
be an equal to that for water. The purpose of the above analysis was 


merely to show the effect of large differences of fluid properties. 


9.4 The Effect of Pressure on Dropwise Condensation Heat Transfer 


The preceding section indicates that both the drop growth and the 
drop distribution vary with pressure. How well one can predict the 
variation of drop condensation heat transfer with pressure by considering 
only the drop growth rate effects is illustrated in Fig. 61. The dashed 
line represents the experimental heat transfer data. The lower solid 


line is the predicted curve varying the properties but keeping the drop 
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distribution (measured at a saturation temperature of 88 °F) constant 
over the complete pressure range. The upper curve varies the themno- 
dynamic properties in Eq. (8.27) while holding the atmospheric pressure 
distribution constant. 

Equation (8.27) predicts the heat transfer correctly only at a Tat 
of 88 °F and 212 °F where the distributions were measured. In between 
these values, the measured and predicted trends agree but the absolute 
values do not. This again confirms the statement that it is necessary 
to know the effect of pressure on both drop growth rate and distribution 


in order to correctly predict the variation of drop condensation heat 


transfer with pressure. 


9.5 The Effect of AT on Dropwise Condensation Heat Transfer 


The effect of AT on drop growth rate has been discussed in section 
9.1 of this chapter. The observation was made that the rapid decrease 
in growth rate for low AT's will result in a decreasing heat transfer 
coefficient. The dashed curves of Figs. 62 and 63 represent the predicted 
H versus AT according to Eq. (8.27)." For these calculations, the drop 
distributions measured at a AT of .5 °F for each pressure were assumed 
to be constant over the range of AT's from 0 to 4 °F. 

Refering to Fig. 62, one can conclude that since the measured and 
predicted curves deviate, the drop distribution varies with AT. The 
predicted curve shows a decreasing coefficient for AT's below .5 °F. 
However, the drop in the predicted coefficient is not as rapid as the 
measured curve. This indicates that nucleation is probably becoming 
important in this low AT range. Thus one can conclude that AT changes, 
ae solid curve represents the curve chosen by least squares technique 


through the experimental heat transfer data (see Section 6.1 for 
details). 
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especially for very low AT's, alter both the drop growth rate and drop 
distribution. 

The AT effect on drop distribution is not as evident for atmospheric 
pressure condensation (see Fig. 63) since the predicted and measured 
curves are nearly the same shape. This indicates that nucleation is 
Mee as much of a limitation for atmospheric as compared to low pressure 


condensation. 


9.6 "Active" Surface Area 

Information concerning the amount of "active’’ surface area is useful 
when discussing dropwise condensation. Table 15 summarizes the results 
obtained by examining Figs. 37, 38, 55 and 56. For low pressure conden- 
Sation, drops smaller than 150 microns occupying about 35 percent of the 
surface area transfer 90 percent of the total heat. At atmospheric 
pressure, drops of diameter less than 40 microns covering 23 percent of 
the surface transfer 90 percent of the heat. In each case about 10 
percent of the surface is bare. 

Fifty percent of the heat is transfered through drops smaller than 
10 microns for low pressure and 4 microns for atmospheric pressure conden- 
sation. In both cases half the heat is transfered through 5 percent of 
the surface area. 

The data on "active" area is useful when estimating the magnitude 
of the constriction resistance developed by Mikic [28]. Sample calcula- 


tions and a discussion of constriction effects are supplied in Appendix G. 


9.7 Minimum Value of the Condensation Coefficient 
For all calculations performed up to this point, it has been assumed 
that the condensation coefficient is equal to one. In order to check 


the validity of this assumption, the coefficient was varied from the 


Tas 
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lowest value reported (.006) up to one. The results of these calcula- 
tions ‘are summarized in Table 16. 

For low pressure condensation, the analysis shows that the coef- 
ficient must be greater than .35. Using the maximum drop distribution, 
L-6, and & = .35, the predicted heat flux is only 98 percent of the 
measured value. Since the true drop distribution must be lower than 
distribution L-6, one can conclude that the condensation coefficient 
must be greater than .35. 

A further discussion on the correct value of condensation coef- 


ficient is supplied in Appendix E. 


9.8 Conclusions 

Since some guess-work was necessary to obtain the drop distribution 
over the complete range of drop sizes, it cannot be concluded at this 
time that the proposed model successfully predicts dropwise condensation 
heat transfer. An analytical prediction or further experimental measure- 
ment of the drop distribution of the small drops is necessary before the 
model can be verified quantitatively. However, the scheme of dividing 
the problem into drop distribution and drop growth effects results a 
a model which explains all the parametric trends observed to date. 

The validity of the information drawn from the heat transfer analysis 
feerning the relative magnitude of the resistances, pressure and AT 
effects, and the minimum value of the condensation coefficient also de- 
pends on the "correctness" of the chosen drop distribution. The author 
feels reasonably confident that the chosen drop distribution approximates 
the true distribution close enough to insure that the discussions in 


the above areas are valid. 


-98- 


X. CONCLUSIONS 


mUet Heat Transfer Measurements 

1. For the three condensing surfaces tested (Teflon, mirror smooth, 
and rough), the heat transfer coefficient rose sharply with AT in the 
range from .1 °F to 1 °F and then remained relatively constant for 
ms to 3 °F. 

2. A critical AT is necessary to initiate the dropwise condensa- 
tion process. This AT rit is .075 °F for mirror smooth surfaces and 
.03 °F for 500 grit roughened surfaces. 

3. For a AT of 3 °F, the heat transfer coefficient decreases 
linearly with the saturation temperature of the vapor according to the 
relation: h = 220 Or eat - 32). For a mirror smooth chemically promoted 
vertical surface, the coefficient decreases from 37,000 Btu/hr ft °F 


hor T pee eee eto 13,000 Btu/hr £t-°r fore. Of OG fs 
sat sat 


mOS2 Visual Investigation 

1. The drop distribution spectrum can be divided into four regions 
according to the governing drop growth mechanism. They have been de- 
fined as the nucleation, condensation-coalescence, coalescence, and 
departure regions. Within each of these regions, the number of drops 
of a particular diameter varies with diameter according to the form: 
N=enpD-. 

2. The drop distribution varies with the system pressure and the 
surface to vapor temperature difference. Excess vapor velocity, sur- 
face roughness and inclination can also alter the distribution. 

3. For a AT of .5 °F during low pressure condensation Cl wae = 88 °F), 


the nucleation site density was measured to be about 200 million per 
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square centimeter. At the same AT for atmospheric pressure conditions, 
Bae density is increased by a factor of about three. Below a AT of 


.5 °F, the nucleation site density drops sharply with AT. 


10.3 Prediction of Heat Transfer 

1. Heat transfer for dropwise condensation can be modeled by con- 
sidering the growth rate of individual drops and the drop distribution 
to be the two important factors affecting the process. The heat transfer 
through a single drop can be calculated from the drop growth rate. The 
average heat flux through the surface can be found by integrating the 
heat transfer through a single drop over the complete range of drop 
sizes. The drop distribution is needed to perform this calculation. 

2. The growth rate for a single drop can be predicted by assuming that 
the resistances affecting the heat transfer act in series. The curvature, 
interfacial mass transfer, and drop conduction resistances have been in- 
cluded in this analysis. The resistances due to constriction and pro- 
moter effects can be added in series with the average resistance due 
to all the drops on the surface. The non-condensable resistance has 
been ignored throughout this analysis. 

3. For low pressure condensation ae = 88 °F), drops smaller 
than 150 microns in diameter occupying 35 percent of the surface area 
transfer 90 percent of the total heat. At atmospheric pressure, drops 
of diameter less than 40 microns covering 23 percent of the surface 
transfer 90 percent of the heat. In each case, about 10 percent of 
the surface is bare. Fifty percent of the heat is transfered through 
5 percent of the surface area for both low and atmospheric pressure 
condensation. 


4. The resistance due to conduction through the drops accounts for 
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70 percent of the total resistance for a vertically oriented, mirror 
smooth, chemically promoted, copper surface during low pressure lies = 
88 °F) condensation. For atmospheric pressure, the drop conduction 
resistance is nearly 85 percent of the total. 

2. The interfacial mass transfer resistance is nearly 8 times 


greater at low pressure (T. = 88 °F) as compared to atmospheric 


at 
pressure. This increase is due solely to property changes. 
6. The condensation coefficient must be greater than .35. A 


value of one was used for all heat transfer calculations throughout 


iis thesis. 


10.4 The limiting Heat Transfer Mechanisms 


1. When determining why a parametric change alters the heat transfer 
for dropwise condensation, one must consider how it affects the drop 
growth rate of individual drops and also how it affects the drop distri- 
bution. Frequently a change in operating conditions has an effect on 
both growth rate and drop distribution. 

2. The decrease in heat transfer with decreasing pressure is due 
to a slowing down of the drop growth rate and an alteration of the drop 
distribution. The changes in thermodynamic properties with decreasing 
pressure result in the lower growth rate. The less desirable drop 
distribution at lower pressures is caused by poorer drop nucleation 
and drop departure. 

3. The rapid decrease in heat transfer coefficient for AT's below 
.5 °F is due to a slowing down of the drop and also to nucleation 
limitations. 

4. The promoter limits drop condensation by offering an added 


conduction resistance in series with the condensation resistance. For 





-101- 


Teflon surfaces this is a serious limitation. 
2. Low thermal conductivity of the substrate causes a non-uniform 
heat flux over the condenser surface and limits the heat transfer by 


introducing a constriction resistance. 


10.5 Recommendations for Further Study 


The determination of the drop distribution over the complete 
range of drop sizes is the most important piece of information needed 
to prove the validity of the heat transfer model proposed in this thesis. 
At very low AT's where the number and size of the active drops are 
greatly decreased, the entire drop distribution could be measured using 
standard microscopic techniques. However, for AT's in the normal 
operating range (above .5 °F), special experimental techniques will be 
required to obtain this data. The analytical prediction of the drop 


distribution in the small size range is equally difficult. Computer 


simulation of the drop condensation process should prove to be a nsec. 
approach. . 





1 


ie. 


2. 


iS 


14. 


i. 


IG. 


eo 2 


REFERENCES 


Schmidt, E., W. Schurig, W. Sellschopp, 'Versuche iiber die 
Kondensation von Wasserdampf in Film-und Tropfenform", Tech. 


Mech. Thermo-Dynam., 1, p 53 (1930). 


Jacob, M., "Heat Transfer in Evaporation and Condensation II," 


MNeenanical Engineers, 58, p 729 (1936). 
Eucken, A., Naturwiss., 25, p 209 (1937). 


Emmons, H., ‘The Mechanism of Drop Condensation", Transactions 
A.I.Ch.E., 35, p 109 (1939). 


Fatica, N., D. L. Katz, "Dropwise Condensation", Chem. Eng. Progress, 
45, p 661 (1949). 


Baer, E., J. M. McKelvey, "Heat Transfer in Dropwise Condensation", 


Delaware Science Symposium (1958). 


Welch, J. F., J. W. Westwater, "Microscopic Study of Dropwise 
Condensation", Int. Dev. in Heat Trans., p 302 (1961). 


Kast, W., "Heat Transfer with Drop Condensation", Chenie~Ingenieur- 
Technik, 31, p 163 (1963), (Translated from German - NSTIC 
Translation No. 1564). 


Sugawara, S., I. Michiyoshi, "Dropwise Condensation", Memoirs of 


Faculty of Engineering, Kyoto Univ., 18, p 84 (1956). 


Silver, R. S., "An Approach to a General Theory of Surface: 
Bendensers , Proc. Inst. Mech. Engrs., 178, p 339 (1963). 


Sugawara, S., K. Katsuta, "Fundamental Study on Dropwise 
Condensation", Third Int. Heat Trans. Conf., 2, p 354 (1966). 


Tammann, Von. G., W. Boehme, "Die Zuhl der Wassertropfchen he der 
Kondensation auf verschiedenen festen Stoffen'', Annelen der Physik, 
peep 22 (1935). 


McCormick, J. L., E. Baer, 'Dropwise Condensation on Horizontal 
Surfaces", Devel. in Mech. (Proc. of 8th Midwestern Mech. Conf. 
1963) (1965). 


McCormick, J. L., J. W. Westwater, "Nucleation Sites for Dropwise 
Condensation", Chem. Eng. Science, 20, p 1021 (1965). 


McCormick, J. L., J. W. Westwater, "Drop Dynamics and Heat Transfer 
During Dropwise Condensation of Water Vapor on a Horizontal Surface", 


Chem. Eng. Prog. Sym. Series, 62, p 120 (1966). 


Peterson, A. C., J. W. Westwater, 'Dropwise Condensation of Ethylene 


Glycol", Chem. Eng. Prog. Sym. Series, 62, p 135 (19'G6)% 





7. 


1S. 


Ie 


ZO . 


Zak, 


Za. . 


ZS 


24. 


25. 


ZO. 


Zi « 


28. 


Zo. 


SU. 


Bl. 


52. 


=1103= 


Baer, E., J. A. Koutsky, A. G. Walton, "The Nucleation of Water 
Vapor on Surface Films Formed by Adsorption on Metal Surfaces", 
plumiacemeuerence, 3, p 280 (1965). 


Koutsky, J. A., A. G. Walton, E. Baer, "Heterogeneous Nucleation 
of Water Vapor on High and Low Energy Surfaces", Surface Science, 
Seep too (1965). 


Erb, R. A., "Heterogeneous Nucleation on Single-Crystal Silver 
and Gold Substrates in Cyclic Condensation of Water Vapor", 
Ph.D. Thesis, Temple Univ. (1965). 


Erb, R. A., "Wettability of Metals under Continuous Condensing 


Conditions", J. of Physical Chemistry, 69, p 1306 (1965). 


Umur, A., P. Griffith, "Mechanism of Dropwise Condensation", 
ASME, Paper No. 64-WA/HT-3 (1964). 


McCormick, J. L., E. Baer, "On the Mechanism of Heat Transfer in 
Dropwise Condensation", Journ. of Colloid Science, 18, p 208 (1963). 


Le Fevre, E. J., J. W. Rose, "A Theory of Heat Transfer by Drop- 
wise Condensation", Proc. Third Int. Heat Trans. Conference, 2, 
p 362 (1966). 


Rose, J. W., "On the Mechanism of Dropwise Condensation", Int. J. 


Heat Mass Transfer, 10, p 755 (1967). 


Gose, E. E., A. N. Mucciardi, E. Baer, ''Model for Dropwise Conden- 
sation on Randomly Distributed Sites", Int. J. Heat Mass Transfer, 
10, p 15 (1967). 


Hurst, C. J., "Transient Droplet Growth During Dropwise Condensation", 
Ph.D. Thesis, Penn. State University (1966). 


Glicksman, L. R., "An Analytical Study of Dropwise Condensation", 
To be published in Int. J. Heat Mass Transfer. 


Mikic, B. B., "On Mechanism of Dropwise Condensation", To be 
published in Int. J. Heat Mass Transfer. 


Sugawara, S., I. Michiyoshi, "Dropwise Condensation", Memoirs 


of Faculty of Engineering, Kyoto Univ., 18, p 84 (1956). 


Hampson, H., N. Ozisik, "An Investigation into the Condensation 


of Steam", Institution\of Mech. Engrs., 13, p 282 (1952). 


Dolloff, J. B., N. H. Metzger, "Dropwise Condensation of Steam 


at Elevated Pressures", ONR Tech Rpt #2, Contract Nonr 3357(02) 
(May 15, 1968). 


Rose, J. W., "Dropwise Condensation - The Effect of Surface 
Inclination", To be published in Int. J. Heat Mass Transfer. 





3. 


34. 


BD . 


B10. 


7. 


38. 


0 


40. 


41. 


42. 


43. 


44. 


45. 


46. 


oe 


48. 


=a 


cemtvemy ob. f., Chem. Eng. Science, 10, p 1 (1959) and 17, p 55 
(1962). 


McCormick, J. L., E. Baer, "On the Mechanism of Heat Transfer in 
Dropwise Condensation", Journ. of Colloid Science, 18, p 208 (1963). 


MacDougall, G., C. Ockrent, Proc. Royal Society, 180A, p 151 (1942). 


Le Fevre, E. J., J. W. Rose, ‘An Experimental Study of Heat Transfer 
by Dropwise Condensation", Int. J. Heat Mass Transfer, 8, p 1117 
(1965). ; 


Krause, J., "Dropwise Condensation - Quantitative Research", 
Heat Div. Paper No. 28, National Eng. Lab., East Kilbride (1950). 


Wenzel, H., "Versuch uber Tropfenkondensation", Allg. Warmetech, 
pee 53 (1957). 


Tanner, D. W., C. J. Potter, D. Pope, D. West, ‘Heat Transfer in 
Dropwise Condensation. Part I The Effects of Heat Flux, Steam 
Velocity, and Non-Condensable Gas Concentration", Int. J. Heat 
Mass Transfer, 8, p 419 (1965). 


Nagle, W. M., G. S. Bays, L. M. Blenderman, and T. B. Drew, "Heat- 
Transfer Coefficients during Dropwise Condensation of Steam", 


Trans. Amer. Inst. Chem. Engrs., 31, p 593 (1935). 


Gnam, E., "Tropfenkondensation von Wasserdampf", Forschungsh. Ver. 
Misch. Img., 382, p 17 (1937). 


Fitzpatrick, J. P., S. Baum, W. H. McAdams, "Dropwise Condensation 


of Steam on Vertical Tubes", Trans, Amer. Inst. Chem. Engrs., 35, 
p 97 (1939). 


Shea, F. L., N. W. Krase, "Dropwise and Filmwise Condensation of 


Steam", Trans. Amer. Inst. Chem. Engrs., 36, p 463 (1940). 


Kirschbaum, E., G. Winckelsesser, A. K. Wetjen, "Neues uber den 
Warmeanstausch", Chem.-Ing.-Tech., 23, p 361 (1951). 


Costas, L. P., "Heat Transfer by Dropwise Condensation of Steam", 
Ph.D. Thesis, Univ. of Maryland (1959); 

Citakoglu, E., J. W. Rose, ''Dropwise Condensation - Some Factors 
Influencing the Validity of Heat Transfer Measurements", Int. J. 
Heat Mass Transfer, 10, p 1267). 

Ma, J. T. S., Ph.D. Thesis, Iowa State College (1959). 

Brown, A. R., M. A. Thomas, "Filmwise and Dropwise Condensation at 


Low Pressures", Third Int. Heat Trans. Conf., 2, p 300 (1966). 





49. 


=). 


pi. 


2. 


53. 


54. 


Ds 


DO . 


Dy 


Do. 


Se 


60. 


ol. 


-105- 


Tanner, D. W., D. Pope, C. J. Potter, D. West, "Heat Transfer in 
Dropwise Condensation at Low Steam Pressures in the Absence and 


Seeoenes Ot Non-Condensable Gas", Int. J. Heat Mass TranstCiye al. 


p 181 (1968). 


Furman, T., H. Hampson, "Experimental Investigation into the 
Effect of Cross Flow with Condensation of Steam and Steam—Gas 


Mixtures on a Vertical Tube", Proc. Instn. Mech. Engrs., 173, 
mmkay (1959). 


Hampson, H., "The Condensation of Steam on a Tube with Filmwise 
or Dropwise Condensation in the Presence of a Non-Condensable 
Gas", Int. Devel. in Heat Transfer, 2, p 310 (1961). 


O'Bara, J. T., E. S. Killian, L. H. S. Roblee, "Dropwise Conden- 
sation of Steam at Atmospheric and Above Atmospheric Pressures," 


Chem. Eng. Science, 22, p 1305 (1967). 


iemner, Do W., C. J. Potter, D. Pope, D. West, "Heat Transfer in 
Dropwise Condensation. Part II Surface Chemistry", Int. J. Heat 
Mass Transfer, 8, p (1965). 


Griffith, P., Man Suk Lee, "The Effect of Surface Finish and 
Material Thermal Properties on Dropwise Condensation", Int. J. 
heat Mass Transfer, 10, p 697 (1967). 


mepper, L., E. Baer, J. Colloid Science, 10, p 225 (1955). 


Babco, R. P., A. L. Gosman, "Promotion of Dropwise Condensation 
of Several Pure Organic Vapors", Am. Soc. Mech. Eng., Paper No. 
57-3-2 (1957). 


Misra, B., C. F. Bonilla, "Heat Transfer in the Condensation of 
Metal Vapors: Mercury and Sodium up to Atmospheric Pressure", 


Chem. Eng. Progress Sym. Series, 52, p 7 (1956). 


Ivanovskii, M. N., V. I. Subbotin, Yu. V. Milovanov, "Heat Transfer 
with Dropwise Condensation of Mercury Vapor", Teploenergetika, 
ip 81 (1967). 


Drew, T. B., W. M. Nagle, W. G. Smith, "The Conditions for Drop- 


wise Condensation of Steam'', Trans. Amer. Inst. Chem. Engrs., 31, 
peel (1935). 


Blackman, L. C. F., M. J. S. Dewar, H. Hampson, "An Investigation 
of Compounds Promoting the Dropwise Condensation of Steam", J. 


Popiemcnen., 7, p 160 (1957). 


Osment, B. D. J., D. W. Tanner, "Promotion for the Dropwise Con- 
densation of Steam", First European Symposium - Fresh Water from 
the Sea, Athens, 1962. 





D2 


tore 


64. 


G5. 


66. 


Oy. 


68. 


G9. 


70. 


7A 


ao 


ior 


74. 


aye 


viol. 


-106- 


Osment, B. D. J., D. Tudor, R. M. M. Speirs, W. Rugman, "Promoters 
of the Dropwise Condensation of Steam", Trans. Instn. Chem. Engrs., 


eve p 152 (1962). 


Peipsomembhe Ge. Il., D. C. P. Birt, C. W. Homour, BOY W. Ashle the 
Promotion of Dropwise Condensation by Montan Wax", J. Appl. 
Siew... L2, p 539 (1962). 


famdenr, D. W., A. Poll, J. Potter, D. Pope, D. West, “The Promotion 
of Dropwise Condensation by Montan Wax. II. The Composition of 
Montan Wax and the Mechanism of Promotion", J. Appl. Chem., 12, 
peod? (E962). | 


Watson, herG. Heed. 0. Brunt, D. C. P. Birt; Drepwise Condensation 
of Steam", Proc. 1961-62 Int. Heat Trans. Conference, p 296 (1963). 


Fox, R. M., "A Review of Literature on the Promotion of Dropwise 
Condensation", US Navy Marine Engineering Lab Report 71 106 (1964). 


Erb, R. A., E. Thelen, "Dropwise Condensation on Hydrophobic Metal 
and Metal Sulfide Surfaces", Amer. Chem. Soc., p 100 (1965). 


Bromley, L. A., J. W. Porter, S. M. Read, “Promotion of Drop-By- 
Drop Condensation of Steam from Seawater on a Vertical Copper Tube", 
Peon n. Journal, 40, p 245 (1968). 


morenz, J. J., B. B. Mikic, ‘The Effect of Thermocapillary Flow on 
Heat Transfer in Dropwise Condensation", To be published in the 
Trans. ASME Journ. Heat Transfer. 


Nabavian, K., L. A. Bromley, "Condensation Coefficient of Water”, 


faem. Eng. Sci., 18, p 651 (1963). 


Lord Kelvin, Proc. Roy. Soc. Edin., 7, p 63 (1870), and Phil. Mag., 
feeep 440 (1871). 


von Helmholtz, R., Wied. Ann., 27, p 508 (1886). 


Hatsopoulos, G. N., J. H. Keenan, Principles of General Thermo- 
dynamics, John Wiley and Sons, Inc., New York (19oa)e 


Shrage, R. W., A Theoretical Study of Interphase Mass Transfer, 
Columbia University Press, New York (1953). 


Adt, R., "Bubble Formation and Growth, A Study of the Boundary 
Conditions at a Liquid-Vapor Interface through Ineversible Thermo- 
dynamics", Sc.D. Thesis, M.I.T. (1966). 


Wilcox, S., "Film Condensation of Liquid Metals - Precision of 
Measurement", Ph.D. Thesis, M.I.T. (1969). 





ia . 
78. 


9. 
80. 
el. 
B2), 
oo. 


84. 


-10/7- 


Petter, Moyea. Weber, Z. Physik. Chem., 119, p 277 (1925). 


Becker, R., W. Doring, Ann. Physiikigu24, pe 719 eC tO ab). 


Hirth, J. P., G. M. Pound, Condensation and Evaporation, Nucleation 
and Growth Kinetics, Pergamon Press, Oxford (1963). 


Adamson, A. W., Physical Chemistry of Surfaces, Interscience 
Publishers, Inc., New York (1960). 


Sigsbee, R. A., G. M. Pound, "Heterogeneous Nucleation from the 


Pao weedvan- Colloid Interface Sci., 1, p 335 (1967). 


Gibbs, J. W., Collected Works, Vol. I, Thermodynamics, Yale 


University Press, New Haven (1948). 


Volmer, M., Kinetik der Phasenbildung, Edwards Brothers, Inc., 
Ann Arbor, Mich. (1945). 


Chakraverty, B. K., G. 


jks 


Pound, Acta Met. , 12. pecotmChJG.). 





-108- 


APPENDIX A HEAT MEASUREMENT TECHNIQUES 


The precautions the author took to ineutemte cunts heat transfer 
measurements will be discussed in this appendix. Chapters 4 and 5 have 
emphasized how the presence of non-condensable gases and errors in 
temperature measurement can drastically affect dropwise condensation 
heat transfer experiments. Table 9 outlines the major steps taken be- 


fore the required accuracy was obtained. 


A.l Errors Due to Non-Condensable Gases 

The original condensing chamber and boiler section of the apparatus 
was poorly designed. As described in Section 5.1, the original system 
was large and contained numerous potential dead air spaces. In addition, 
the vapor output from the boiler was small resulting in a low vapor 
velocity past the condensing surface. 

The preliminary experiments run with the original system resulted 
in very low measured heat transfer coefficients. At a heat flux of 
30,000 Btu/hr re ane heat transfer coefficient was approximately 1,500 
Btu/hr et?°r, Temperature fluctuations as great as 1 °F were observed 
at the surface. The low coefficients and the temperature fluctuations 
are indicative of non-condensable gas effects. 

An attempt was made to reduce non-condensable gas effects in the 
original system by blocking off part of the chamber and increasing the 
output of the boiler. These modifications resulted in a doubling of 
the measured heat transfer coefficient. For a heat flux of 30,000 


2 


Btu/hr ft~, the AT decreased from 20 °F to 10 °F. This error in the 


measured AT of 10 °F can be solely attributed to non-condensable gases. 





-109- 


Rather than continue to modify the original system, a completely 
new condensing chamber and boiler were designed and constructed. This 
feo eeysten, much smaller than the original one, is fully described in 
Section 6.1. Tests made with this system resulted in a measured coef- 
ficient comparable to those by Rose [36,46] and Tanner [39] and their 
colleagues (see Figs. 9 and 10). With this system, the author feels 
confident that the heat transfer measurements were free of non- 


condensable gas effects. 


A.2 Temperature Measurement Techniques 


e e 2 2 
Due to the high heat transfer coefficients (up to 50,000 Btu/hr ft’ °F) 
encountered in dropwise condensation, extreme accuracy in temperature 
measurement is required. This section will mention the precautions taken 


in this experiment to obtain this accuracy. 


A.2.1 Thermocouple Arrangement 


The vapor Sonera” and the temperature at seven locations with- 
in the test section were measured with .005 inch diameter copper- 
constantan thermocouples. The nine thermocouples were made from the same 
spool of highest precision thermocouple wire and installed in the system 
with identical thermocouple circuits (see Fig. A-1). All nine cold 
junctions were placed in individual glass tubes immersed in crushed 
ee. 

Before and after each run, the output of the thermocouples were 
compared at room temperature. All agreed to within .05 °F. A newly 
calibrated Leeds and Northrup K-2 potentiometer was used to measure 
Ene thermocouple output. 


* 
Two thermocouples were used to measure the vapor temperature. 
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hee Installation of Thermocouples 


Significant errors in temperature measurement can result from 
improper thermocouple installation. One of the most important errors 
is due to heat conduction along the thermocouple lead. Even if the 
thermocouple leads are positioned in an isotherm, there can be an error 
due to insufficient immersion depth. 

Figure A-2 illustrates the magnitude of the error that can result 
from insufficient immersion depth of thermocouples. The key to proper 
thermocouple installation is to maximize the heat transfer from the in- 
side of the hole to the thermocouple junction as compared to the conduc- 
tion heat transfer from the junction to the leads outside the hole. 

The ideal situation is to have long, low conductivity thermocouple wire 
and a tight fit between the junction and the hole. 

A temperature drop of 10 °F was experienced when a copper constantan 
thermocouple wire of diameter .01 inch was traversed through the .034 
inch thermocouple holes. The temperature versus distance plot is in- 
dicated in Fig. A-2(a). When a .005 inch diameter wire was used doubled 
over once, the temperature drop was decreased to 5 °F. It was not until 
the effective immersion depth was increased to 1.6 inches with the .005 
inch wire that a flat temperature gradient was measured (see Fig. A-2(c). 

The installation technique used in the experiment is illustrated 
in Figs. A-2(d) and 19. The thermocouple wire was doubled over twice 
to give an effective immersion depth of 1.6 inches. The hole was filled 
with an anti-seizing compound called "silver goop" to enhance heat trans- 
fer from the inside of the hole to the thermocouple junction. The ends 
of the thermocouple holes were plugged with heavy copper wire. The 


resulting error due to conduction along the lead was thought to be less 
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han .05 °F. 
‘To insure that there would be no error due to conduction along 
the thermocouple wires used to measure the vapor temperature, three 
inch long leads were used inside the condensing chamber. The wires in- 
side the glass tubes used for the cold junctions were immersed 12 inches 


me thie ice bath. 


A.2.3. Error Due to Location of Thermocouple Junction 


Wilcox [76] has carried out a detailed analysis to determine the 
temperature measurement error due to uncertainty of the precise location 
of the thermocouple measuring junctions within the hole. Figure A-3 
shows how uncertainty in temperature measurement within each hole can 
result in a sizeable uncertainty in the determination of the wall 
temperature. 

Wilcox has shown that the error in wall measurement can be de- 
creased by using a large number of widely spaced, small thermocouple 
holes and a high conductivity bulk material. Since seven thermocouples 
(with spacing between the two most distant thermocouples of 1.8 inches) 
were embedded in the copper test section used in this experiment, high 
accuracy in the measurement of the wall temperature was obtained. 

The details of the analysis of Wilcox will not be repeated here. 
Figure A-3 shows that a Gaussian distribution of possible temperature 
measurements about the correct one at the center of the hole, iT: was 
chosen. The resulting temperature distribution at the wall will be 
Gaussian with the magnitude of the standard deviation depending on the 


number, spacing, and size of the thermocouple holes and the conductivity 


of the bulk material. 
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The error in wall measurement for the test section used in this 
experiment was calculated from the analysis of Wilcox. The resulting 
error in wall measurement divided by the heat flux (CT ~ T)/Q/A) was 
meO X os °F/Btu/hr ae This analysis showed that the error in wall 
measurement was always less than 10% of the measured temperature 


difference (T -T ). 
S W 


A.2.4 Conclusions 

It is not possible to estimate exactly how efficient the final 
design of the apparatus was in ridding itself of non-condensable gas 
effects. The author feels certain that the high velocities past the 
condensing surface were sufficient to remove a build-up of gases. Thus 
it is thought that the results reported in this thesis are for "gas free" 
steam; 

A typical plot of the measured linear temperature gradients are 
plotted in Fig. A-4. It is believed that the total error in temperature 


measurements was always less than 10 percent of the measured vapor to 


surface temperature difference. 
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APPENDIX B DISCUSSION OF DROP DEPARTURE MECHANISMS 


fel Departure Criterion for a Static Droplet 


The subject of drop departure has been discussed in papers by 
Fatica and Katz [5] and Sugawara and Michiyoshi [9]. In these two 
papers, models have been proposed which predict the maximum size drop 
which can exist in static equilibrium on an inclined surface. In both 
of these models, no vibrational motion of the drop due to coalescence 
is considered. 

Figure B-l illustrates the forces acting on a droplet. It is 
evident from this figure that the values of the contact angles (advancing 
and receding), the surface tension, and the weight of the drop determine 
its departure size. 

Just before a drop begins to depart, components of the forces in 
the plane of the surface can be equated to give a departure criterion. 


The resulting force balance is 


S cos 9 = W sin 8 + S cos @ 
r a 


Bare 
where ( ) 


W = ogv 


The volume of a drop depends on its shape (contact angles) as well 
as its diameter. Fatica and Katz and Sugawara and Michiyoshi have 


shown that the volume of a drop can be computed from the relation 
72 Dee Onn. (B.2) 
ar 
Substituting (B.2) in (B.1), the maximum drop size can be found to be 


p= i= eee (B. 3) 


max og sing 
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Here £(0) is a functional relationship including both the advancing and 


receding contact angles. According to Fatica and Katz, 





0 
1 + — 
9 pa 2 cos 0. - cos 8 
£(0)° = 5 ; = (B.4) 
tan — Eon 84/2 
2 ———-—-———— £(0) + £(6) 
tenn 0/2 
and to Sugawara and Michiyoshi, 
16 sin"8 0 - 8. 
E(o) = 9 790 agin 20° 
Mm m 
(Boo) 
where 6 acid 
ea 
m 2 


Fatica and Katz summarized their results for the case of a water ~ 
droplet on a vertical surface in graphical form (reproduced here as Fig. 
B-2). From this figure one sees that the maximum drop size decreases 
as the advancing and receding contact angles are increased. This 
figure may be used for other fluids and different temperatures by 


scaling the ratio of the two property terms, - 


B.2 Departure Criterion for a Dynamic Droplet 


When two drops coalesce, they form a single drop located at the 
center of gravity of the two original drops. Thus, coalescence causes 
the drops to move toward each other. When a large "dead" drop is sitting 
on a surface, it is constantly being bumped into by the small "active" 
drops about its periphery. Each time a coalescence occurs, the large 
drop moves slightly. The random nature of the small drops coalescing 
into the big drop will set the big drop into some sort of a random 


vibratory motion. 
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As the condensing rate is increased, the rate at which active drops 
grow and bump into the large drops increases. Thus, it can be assumed 
that the large drops will vibrate faster at higher heat fluxes. 

The motion imparted to the big drops by the coalescing action makes 
them less stable than static drops. Thus, the dynamic drops will depart 
from the surface prematurely. The lowering of the departure size due to 
the dynamic effects has been observed by Brown and Thomas [48] and this 


author. 


B.3 Vapor Velocity Effect on Departure Size 


This author and several other experimenters [39,52] have noticed 
that high vapor velocities past the condensing surface increase the 
heat transfer coefficient. The explanation for this increase in per- 
formance is the early removal of departing drops. 

A high vapor velocity across the condenser surface introduces a 
shear force on the large drops which tends to sweep them from the surface. 
The appropriate component of this shear force can be added to the force 
balance of Eq. (B.1) and the departure size predicted. Most experimenters 
have blown the vapor across the condensing surface perpendicular to the 
downward track of departing drops. The velocity effect would be more 


dramatic if the vapor flow was directed in the downward direction. 


Peeeeriftect of Vibration of the Surface on Departure Size 


The preliminary results of the surface vibration study mentioned 
in Section 6.7 demonstrated an increase in heat transfer performance 
with increased frequency. Mechanically vibrating the surface introduces 
an acceleration force on the drops which tends to shake them from the 


surface. The resulting early departure is similar to that caused by 
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the dynamic effect of the coalescing action. 

The direction of the surface vibration with respect to the departure 
direction as well as its amplitude and frequency will determine the 
effectiveness of this heat transfer enhancing technique. The magnitude 
of the decrease in departure size can be predicted by introducing the 


acceleration force into Eq. (B.1). 


B.5 Conclusions 

This analysis of departure mechanisms is by no means complete. It 
would be very beneficial to be able to predict the drop distribution for 
the large departure size drops on the surface. How changes in the distri- 
bution of large drops affect dropwise condensation heat transfer is 


another area which needs further investigation. 
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APPENDIX C DISCUSSION OF CURVATURE RESISTANCE 


C.1l° Kelvin - Helmholtz Equation 


The equilibrium between two phases of a pure substance across a 
curved interface is different than that for a flat interface. Lord 
Kelvin [71] and von Helmholtz [72] have shown that the curvature of the 
surface separating the two phases as well as the temperatures determine 


the equilibrium pressure. Thus, 


Bre =) fCl or : (Esi5) 


The relationship stated in Eq. (C.1) can be derived by refering to 
fae. C-l. 

Assume that the system illustrated in Fig. C-l is at a uniform 
temperature despite the variation in pressure caused by hydrostatic 
forces. For a liquid which does not wet the capillary tube, the liquid 
vapor interface will be a curved surface at a level of height L below 
the flat surface. This system simulates drop-like curvature. 


The pressure of the liquid below the curved interface is 
Deere Peumear = (C. 2) 
For dynamic equilibrium across the curved interface 
mr’ ae + 27rd = mr’ Pro : (C23) 
For pressures well below the critical pressure, 
PP’ ee 


gc S 


since the hydrostatic pressure eee is small. Thus, 
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Peo 7 Ps / Pée 7 Pec i (C.4) 


Substituting (C.4) in (C.3) and rearranging terms, 


Dee. = ae : (C.3) 
fembining (C.2) with (C.5), 
_ 20 
or 
ZG 
= —-=—— Gia 
co (C.6) 


For an element of vapor, dy, in the capillary tube, 


dp = ce gdy . (C.7) 


Assuming the perfect gas law is valid and since the temperature of the 


system is uniform, then 


P 
oa = om : CGes) 
& gs 
Be@estituting (C.8) in (C.7), 
) Se = oO gdy . (G22) 
S p gs 


Integrating Eq. (C.9) from y = 0 to y = L gives 


Pp 
p & —= = p_ gL. (C.10) 
Ps 
Again assuming that the perfect gas law holds, 
Pp = 
—= = RT. (C.11) 
0 S 
gs 


By combining Eqs. (C.10) and (C.11), the following relationship is found: 





Sl die 


Pp 
oa = 25, . ( Gaab2)) 
Equation (C.12) is known as the Kelvin-Helmholtz Equation. Since 
r > 0 for a drop, then this equation states that the pressure of the 
vapor above a curved interface is greater than the pressure above a 
flat interface. The analysis is equally valid for the case of a liquid 
which wets the capillary tube and rises to a concave liquid surface. 
This situation simulates bubble-like curvature. 
The Kelvin-Helmholtz equation is displayed graphically in Fig. C-2. 
The dash lines give the pressure of the vapor in equilibrium with drops 
of radius, r. The dot and dash lines give the pressure of the liquid 


in equilibrium with bubbles of radius, me 


Peoeecemperature Difference for Drop-Like Curvature 


The equilibrium condition for the curved interface is compared 
to the flat interface in Fig. C-3. The two systems are at the same 
temperature, To? but not the same pressure. For system 1, the system 


eeu? Pa. » as read from the pressure gauges, is equal to the satura- 
il 


tion pressure, Ps» corresponding to the uniform temperature, T° However, 


for system 2, the system pressure, ae » 1s equal to the vapor pressure, 
Z 


eee in equilibrium with the liquid drop at the temperature, TS 


The Kelvin-Helmholtz equation states that 


P 

tn Bo a (c.13) 

Ps Pr O 

ea De 
Since ——— 1, then 

S P P - p 

ee x ae AS e (C214) 
Ps Pg 


SUbsmrmutiime Cc. 14) in (C.13), 
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20 P. 


Pers Oe ; (Can5) 
gc Ss rp, RT, 


The saturation temperature corresponding to the pressure of system 


zy ee can be found from the Clausius-Clapyron equation 


dT | Veg*o 
ie 
‘ i ae 
If it is assumed that eo is constant from p_ to p then, 
a S gc 
ae | ee 
O SG Weare. <0 
Ss gc fg 


20 P, VY, Ty 
ore - Ty = rp, RI, a 2 (C.18) 


fg 
Pe ‘fe 
For a perfect gas, RT = i: ThuseEq. (Gi lS)setmpieties to 
: 20 Ty 
ae - To =: rH ; (Cory) 


re is the saturation temperature corresponding to the pressure measured 


for system 2. Ty is the saturation temperature of system 1. Thus, 


ZO 
aa 
e = = a ae : (E220) 
t H 
sat. sat, © fe OP, 
Here ee is the saturation temperature above a curved interface and 
2 
Ce is the saturation temperature above a flat interface. Equation 
1 


(C.20) states that the vapor above a drop must be supersaturated compared 
to the vapor above a plane interface at the same temperature in order 


for the drop to be in thermodynamic and mechanical equilibrium. 
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The two systems illustrated in Fig. C-3 are at the same temperature, 
Ee The systems can be modified so that the pressure instead of the 
temperature is kept constant. If the pressure is the same at both the 
curved and the flat interface, then the saturation temperature corre- 
sponding to the curved interface will be less than the saturation tempera- 


ture at the flat interface. The resulting temperature difference is 


given by Eq. (C.20). 


Setting AT. = sat, - "sat, ’ then the minimum size drop which can 
exist on a surface is 
Z0 T. 
1 2 ae : (G.218) 


min ge Pr. AT. 


A simple expression for the temperature drop due to curvature can be 


found by combining Eqs. (C.20) and (C.21): 


AT = AT : (GAZ 





This temperature difference is illustrated in Fig. 44. 


C.3 Temperature Difference for Bubble-Like Curvatures 
For the sake of comparison, the temperature drop due to bubble- 
like curvature will be derived. Refering to Fig. C-4, one sees that the 

pressure of system 2 is now the pressure of the liquid. (For the case 
of drop-like curvature, the system pressure was the pressure of the 
vapor phase.) 
The pressure of the vapor inside the bubble can be found from the 
Kelvin-Helmholtz equation, 
20 P, 


po o-pFeooeS SCO. (G2 23) 
gc S XP, RT. 


The pressure of the liquid near the interface is 


e 


- —_ «& Se ee 
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Combining (C.23) and (C.24), the liquid pressure (and thus the system 


pressure) is 


2 
2 Ps 20 


2 oe kts ) Seze) 


Since p,/P-RT ol, thenseq-. (C.23) Simplifies to 


20 


Be — Pee = om , (C, 26) 


The Clausius-Clapyron relation can be applied to (C.26) in a 


similar way as in Eq. (C.18). The resulting temperature difference is 


20 “FE Ty 
T -T = = —&— ., (cay 
oO ic r H 
fg 
For pressures not near the critical point ee ~ Me = Lie and thus 
ime (C.27) becomes 
20 TS 
TI -T = =a : (C28) 
O fc r H 
fg "g 


Te is the saturation temperature corresponding to the pressure 


measured for system 2. The saturation temperature for system l is TS" 


Thus, Eq. (C.28) may be rewritten as 


20 T 
Sat, 
ic -T = a : (CG5 28), 
Sat, sat, re ee oe 
iT is the saturation temperature above a curved interface and T 
sat, ea 


the saturation temperature above a flat interface. os TSePe tempera — 
ture of the vapor within the bubble. 
Equation (C.29) states that the vapor inside a bubble must be 


superheated as compared to the saturation temperature of the liquid. 








DO eee a eee ae 


=e 


Since the bubble AT involves the reciprocal of the vapor density (as 
compared to the reciprocal of the liquid density for the droplet AT), 
Ene temperature drop predicted by Eq. (C.29) is large compared to that 
pmedicted in Eq. (C.20). | 


The minimum bubble size will be 


20 Te 
tr. = —_— 
min ee De AT. 


where (C.30) 


eae Eifect of system Pressure on Minimum Size Drop and Bubble 


The minimum size drop which can exist in equilibrium is 





40 T 
D . = ae ® (C235) 
nin) ae Pe AT. 
For a bubble, the minimum size is 
46 T. 
Dae = SS ae : (C.32)) 
min, ee a AT. 
At the same operating conditions, 
Die 
min 0 V 
Se oS, (Cc. 33) 
dD. QO Vv 
min, f g 


This relation shows why the minimum drop size is smaller by more than 
three orders of magnitude than the minimum bubble size. 

The effect of varying the pressure of the minimum drop and bubble 
size is displayed in Fig. C-5. It is evident from this figure that 
pressure has little effect on the minimum size of drops but a large 


effect on minimum bubble sizes. The reason for this is that the liquid 
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is nearly independent of pressure while the vapor density is 





y pressure dependent. 
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APPENDIX D DISCUSSION OF NUCLEATION 


Heat transfer engineers encounter the nucleation phenomenon in 
nucleate boiling and dropwise condensation. For the case of boiling, 
the nucleation criterion has been well developed and the predicted 
trends have been experimentally verified. However, droplet nucleation 
has not been as well covered. Because of the great difference in the 
size and Eber of nucleation sites, droplet and bubble nucleation are 
not entirely analogous. 

For both nucleate boiling and dropwise condensation, the nucleus 
has been observed to form at imperfections such a5 cavities or scratches 
on the surface. For boiling, the size and shape of the cavity which 
will produce a stable vapor nucleus for a given system pressure and 
surface to liquid temperature difference has been predicted and verified 
by experimental observations. If the same analysis is used for drop 
condensation, the critical size cavity can be predicted. At atmospheric 
pressure, the cavity size for droplet formations is nearly 2000 times 
smalier than for bubble formation at the same AT. As the pressure is 
lowered to 34 mm Hg, the ratio of critical bubble size to droplet size 
increases to 40,000. 

The purpose of this discussion is to show that one can expect sur- 
face chemistry to play a vastly different role in the nucleation process 
for dropwise condensation as compared to nucleate boiling. The arguments 
concerning the size, shape, and number of surface imperfections which 
serve as nucleation sites for nucleate boiling are not likely to be 


appropriate for dropwise condensation. 


Because of the disparity between the nucleus size for boiling and 
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condensation, the standard approach used to predict bubble nucleation 
will not be used in this discussion of droplet nucleation. Instead, 
the author will revert back to the more classical approach of comparing 
heterogeneous nucleation of a droplet on a substrate to homogeneous nu- 


cleation of a droplet from a uniform temperature vapor. 


D.1 Homogeneous Nucleation 


feveracé of homogeneous vapor-liquid nucleation can be found in 
the works of Volmer and Weber [77], Becker and Doring [78], Hirth and 
Pound [79], and many others. Much of the material presented in this 
appendix was taken from the reviews of nucleation by Adamson [80] and 
Sigsbee and Pound [81]. 

Gibbs [82] determined that there is a certain free energy of forma- 
tion involved in the transformation of vapor to a liquid droplet. The 
reversible work required to form a droplet from the vapor phase involves 


a volume and surface contribution and can be expressed as 


Ne = 4) 3nr AG. + ane? (D.1) 
where 
eee ey (D.2) 
V be - Ps 


A plot of AG versus radius is given in Fig. D-l. By maximizing 
Eq. (D.1) with respect to r, the critical size nucleus, r*, and the 
Gibbs free energy of formation of a nucleus of critical size, AG*, can 


be found to be 


r* = - _20 (D.3) 


and 





== 





~ (D.4) 


By substituting (D.2) into (D.3) and rearranging terms, the following 


equation is obtained: 


p_ 20v . 
Qn : - Tar (1 25) 
Ss 


Equation (D.5) is simply the Kelvin-Helmholz equation for the effect of 
curvature on vapor pressure (see Appendix C). AG* can be considered an 
"energy barrier'' which must be overcome in order for nucleation to take 
place. As AG* is lowered, nucleation becomes easier. Figure D-2 shows 
that increasing the supersaturation ratio, S = p/P,» lowers the critical 


free energy of formation and thus causes more nucleation to take place. 


wee Weterogeneous Nucleation on a Plane Substrate 


The free energy of formation involved in the transformation of 
vapor to liquid in the presence of a solid surface is lower than for 
the case of no third phase present. The equation for homogeneous 
nucleation can be modified by taking into account the contact angle 
the liquid nucleus makes with the plane substrate. The resulting critical 


radius and Gibbs free energy of formation as derived by Volmer [83] are, 


_ _ 20 

Bee CDG) 

Vv 
and 
16 O 

Ac,* a Tt AG >(8) (D.7) 

where 
eB 
Mey (2 - 3 cose + cos 8) (D.8) 


The result of this analysis is that the radius of curvature of the 





-128- 


critical nucleus is not affected by the presence of the solid interface. 
It is a function only of the surface tension and the Gibbs free energy 
difference, AG. However, the “energy barrier" is lowered. The contact 
angle function, $(8), plotted in Fig. D-3, varies from 1 for 6 = 180° 
to 0 for 6 = 0°. When the liquid completely wets the surface (0 = 0°), 
there is no energy barrier to nucleation, whereas for the other extreme 
(6 = 180°) the energy barrier is the same as for homogeneous nucleation. 
The effect of contact angle on nucleation is illustrated in Fig. 
D-4. At the same supersaturation ratio, lowering the contact angle de- 


creases the energy barrier to nucleation. 


D.3 Heterogeneous Nucleation on a Substrate at a Ledge 
The Gibbs free energy for the formation of a nucleus at a ledge is 
lower as compared to that on a plane surface. Chakraverty and Pound [84] 


determined that the critical Gibbs free energy at a 90° ledge is 


3 
r 
AG, = Bn g 5 
AG 
V 





o, (8) (D.9) 


where 
9 9 9 he 2 
o, (8) = al (sin 6 - cos 8) + (2/1T)cos 8(sin'8 - cos 86) 


4 (2/t™) (cos 8 Bee apie 8) - cos 8 Bice (D.10) 
r*sin 0 ; ; ; 1/2 

- (2/1r*) sin “[r* cosO/(r*° - Y) } | 
r*cos 86 


Y is a variable of integration. 


x 
When @ < 45° and 50° < 8 < 105°, AG, will be less than AG,*. Thus, 


providing the contact angle is within these ranges, nuclei will form 





ee 


first at ledges and other surface faults. 

The effect of a substrate on nucleation depends on its surface 
Smerey. One Can consider the substrate as contributing its surface 
energy to help overcome the energy barrier to the formation of a stable 
nucleus. The higher the surface energy, the less "push'"’ is needed in 
the form of a AT to initiate nucleation. Surface flaws such as scratches 
and pits have a higher "energy" than the average for the entire surface. 
This is why the nucleation rate is greatestat imperfections on the 


surface. 


D.4 Nucleation Equation in Terms of AT 


Heat transfer engineers prefer to work with a temperature difference 
as opposed to a supersaturation ratio p/p: By using the Clausius- 
Clapyron equation, the term 2n p/P, can be written in terms of AT. 


The Clausius-Clapyron equation can be written as © 


Vv iL 
qr _ _fg s . (D.11) 
dp H 
fg 
Ve T 
If it is assumed that a is constant from Pe to p, then 
fg 
T-T Ve T 
Ss. aaa . (D.12) 
P Ee fe 
2 ae 
et 5 < 1, then 
Ss Pp - Pp 
gnP ~ ——— : Goal) 
Ps Ps 
Combining (D.12) and (D.13), 
(T - T )H. 
Qn Ps ee ; (D.14) 
Ps Ps “gs fg 


If the vapor is considered a perfect gas, then 
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pee . (D.15) 


CDEio) 


where AT =T - To: 

The nucleation equations outlined in the previous sections can be 
rewritten ' terms of a temperature difference by using Eq. (D.16). 
Performing this substitution results in the following equations for 


homogeneous nucleation: 


~H 
0 ee (D.17) 
Vv V i 
£ Ss 
20 Ve T. 
Se rare 
a Wa (D.18) 
ie 
and : . 
Ack = x no ¢ ap (D.19) 
Pe M Es 


For heterogeneous nucleation on a flat substrate the critical 


Gibbs free energy of formation is 
* 
AG ‘ia AG* $(86) 
where (D.20) 


3 
2 - 3 cos 0 + cos @ 
6(8) = (2 = 3 cos 2H cos 
In the case of heterogeneous nucleation on a ledge then the contact 
angle terms o, (8) » as specified in Eq. (D.10),should be substituted in 


Eq. (D.20). 


D.5 Nucleation Rate 


Volmer, Becker, and Doring [78 and 83] have determined what the 
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nucleation rate should be for a given nucleus of critical size. Details 
of this rather extensive analysis can be found in the above references. 


The nucleation rate for heterogeneous nucleation is 


_-AG*}/ RI 


eee (D.21) 


In this expression, J is the number of nuclei per cubic centimeter 
forming per second. For first approximations, A can be considered a 
constant independent of pressure and temperature. 

The term AG*¢/RT isthe quantity which must be evaluated to determine 
the rate of nucleation of a liquid droplet from a supersaturated vapor. 
The appropriate contact angle function (Eq. D.8 or Eq. D.10) should be 
used depending on how one models the surface. 

No attempt will be made in this survey to quantitatively predict 
nucleation rates as parameters are varied. The author does, however, 
want to specify the trends which can be expected when the system pressure 


(at a constant AT) and the AT (at a constant pressure) are varied. 


D.6 Effects of Pressure and AT on the Nucleation Rate 
In order to determine the parametric effects on nucleation, the 
term AG*)/RT must be examined. From Eq. (D.19) this term may be 


1 
written as 


Goal 
AG*o 16 S 
7) °°”6CU GUD (D. 22) 
RO H AT 
i@ fg 


where ¢ is given by Eq. (D.8) or (D.10). 
Decreasing AG*}/RT has the effect of increasing the nucleation 
rate. Equation (D.22) shows that as the AT is increased, AG*$/RT 


decreases and the nucleation rate increases. This is certainly the 
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expected trend, as AT can be considered the driving force which over- 
comes the energy barrier to nucleation. 

The pressure effect on nucleation can be realized by examining 
Fig. D.5. Keeping the AT constant, the term AG*/RT decreases with 
increasing pressure. This predicts superior nucleation at higher 
pressures. | 

The contact angle function, 9, also decreases with increasing 
pressure. Figure D.6 shows that for a given liquid-substrate combina- 
tion, the contact angle decreases with temperature. Since it has been 
shown in Fig. D-3 that the contact angle function, $(6), decreases with 
contact angle, then it can be concluded that the contact angle function 
decreases with pressure at a given AT. 

Since the data of contact angle versus temperature is limited and 
since the precise contact angle function for an irregular surface is 
unknown, the relation between AG*¢/RT and T will not be drawn in Fig. 
D-5. However, it can be concluded that a curve of AG*O/RT versus T 
would have a sharper decreasing slope than the one drawn for AG*/RT. 


Thus, for heterogeneous nucleation at a constant AT, the nucleation 


rate can be expected to increase with pressure. 
{ 


D.7 Conclusions 

A typical condensing surface has nucleation sites of varying degrees 
of "potency". At a constant pressure, when the surface to vapor tempera- 
ture difference is increased, more and more of these nucleation sites 
become active. In a similar way, at constant AT, increasing the pressure 
activates an increased number of sites. How rapid the number of active 
nucleation sites changes with AT and pressure depends on the characteris- 


tics (i.e. surface energy) of the surface. 
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APPENDIX E DISCUSSION OF INTERFACIAL MASS TRANSFER RESISTANCE 


Past experiments have observed that there is a finite temperature 
difference at the liquid-vapor interface during the condensation and 
evaporation processes. The magnitude of this temperature difference 
can be computed by examining the mass transfer phenomenon at the inter- 
face. The most complete coverage of this subject ean ne found in the 
book by Schrage [74]. Recent works by Nabavian and Bromley [70], Silver 
[10], Umur and Griffith [21], Le Fevre and Rose [23], and wens | have 
interpreted Schrage's equations. 

The above references derive the expression for the interfacial 
heat transfer coefficient, One by considering the mass flow equations 
at the interface. The mass flux toward and away from a condensing 


surface are respectively: 


™ a— ye 
W = Py (M/27RT_) (E.1) 
and 
Fa 2 
am abe (M/2mRT ) (E.2) 
A net mass flux, w, may exist toward or away from the surface, 


Here a is the "condensation coefficient" or "evaporation coefficient." 
The condensation coefficient is defined as the fraction of vapor mole- 
cules which enter the liquid phase after striking the liquid surface. 


Substituting (E.1) and (E.2) in (E.3), one gets 





* BD 

No attempt will be made at this time to present a complete literature 
survey of this topic. Doctoral theses at M.I1.T. By Adi 75] and 
Wilcox [76] cover the interface problem in great detail. 
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mie 
w = ap. - P,,) (M/27RT _) : (E.4) 
Equation (E.4) can be written in the form of temperature differences 
instead of pressure differencesby utilizing the integrated form of the 
Clausius~Clapeyron equation: 


H 
= — E saa 
P Pp ae = (T. T,) . (E.5) 


V Ss es 


The interfacial heat transfer coefficient is defined as 


W H 
he ie : (E.6) 
> V 


Combining (E.4), (E.5) and (E.6) gives 





7 2. Hy 
h, = “OnE Tv ve : (E.7) 


For net condensation, the vapor has a finite flow toward the 
surface. This would cause more vapor molecules to stick to the surface 
than would occur for an equilibrium system. 

To account for this effect, Schrage determined that the condensa-~ 


tion coefficient used in Eq. (E.7) should be 


ee aoe 
i — 2 ee. O e (E.8) 
Thus, | 
lieD H 
ee _fg | (E.9) 


i Ze 271RT 


The value to use for a has been a topic of controversy for several 
decades. The condensation (or evaporation) coefficient as measured in 
over a dozen experiments has ranged from .006 to 1.0, depending on the 
pressure. It has been shown conclusively in the recent research by 


Wilcox [76] that the low values for a at high pressures can be attributed 





> 


to experimental errors. The accepted value used today for the condensa- 
tion coefficient is 1, and is believed to be independent of pressure. 
Assuming @ to be 1, the interfacial heat transfer coefficient can 
be evaluated from equation (E.9). Values for h, as a function of the 
interface temperatures have been calculated and plotted in Fig. E-l. 
Figure E-1 shows that the interfacial heat transfer coefficient 
increases sharply with temperature. It should be emphasized that this 
increase in h, is due solely to property changes. The value for the 
condensation coefficient has been kept constant at 1 throughout this 


Pemperature range. 
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APPENDIX F DISCUSSION OF DROP CONDUCTION RESISTANCE 


- The conduction resistance for a hemispherical drop of radincwer, 
can be approximated by comparing it to the resistance for a right- 
circular cylinder of radius, r, and height, r (see Fig. F-1). Assuming 
that the top and bottom surface temperatures, T and To» are uniform 


for the cylinder, then its conduction resistance can be defined as 





@f.1) 


Since the average conduction path is shorter for the case of the drop, 


then its resistance assuming uniform surface temperatures is 





(a2) 


where C < 1. 

Fatica and Katz [5] analyzed the drop conde ction problem for an 
isolated spherical segment assuming the vapor-liquid interface tempera- 
ture, T,> and the liquid-solid interface temperature, T.> are constant. 
For any spherical segment they found the heat transfer coefficient 


through the drop to be 


k 
Dae = £(0) D (F.3) 


where £(8) is a drop shape factor. They presented the value of the 
shape factor as a function of contact angle, 8, in graphical form (re- 
produced here as Fig. F-2). 

miewcenstant, C, of equation (F.2) is related to the shape factor 


of Fatica and Katz by the relation 


2 
CEG pi 





= oie 


mrom Fig. F—-2 and equation (F.4), C can be found to be .241. 

Mikic [28] analyzed the exact solution of the conduction problem 
performed by Umur and Griffith [21] and found that if C = .25, the 
simple equation (F.2) aioe ely approximated the exact solution. 

McCormick and Baer [13,34] and Rose [24] accepted the value of 
.241 for C in their analyses of the drop conduction problem. This 
author used the value of .25. Since the two approximate solutions are 
so close, the convenience of using .25 instead of .241 was the basis 


of its selection. 





: =ioc= 
APPENDIX G = DISCUSSION OF CONSTRICTION RESISTANCE 


Until the work of Mikic [28], none of the proposed theories for 
dropwise condensation considered the effect of the thermal properties 
of the condenser surface. Three investigators, Hampson and Ozisik [30], 
Tanner et al. [53], and Griffith and Lee [54] have shown that lowering 
the thermal conductivity of the substrate lowers the measured heat 
transfer coefficient for dropwise condensation. In the latter two 
references, the decrease in coefficient going from copper to stainless 
steel was over five fold. Certainly the effect of thermal conductivity 
of the condenser surface must be considered when talking about dropwise 
condensation. 

The model proposed by Mikic states that the cause of the decrease 
of the measured coefficient for lower conductivity substrates is the 
non-uniformity in surface temperature. Because of the variations in 
size and spacing of drops on the surface, the local surface to vapor 
resistance varies over the condenser surface. Consequently, the surface 
temperature and thus the heat flux is also non-uniform over the surface. 

The effect of non-uniform heat flux is illustrated in Fig. G-l. 

At any point on the surface, the temperature driving force is not the 
measured (T. - T but the local temperature difference (T. ~ T om” 
Mikic calls the resistance due to non-uniformities in heat flux a 
"constriction" resistance. If the drop condensation rate is to be 
considered as that due to the temperature difference (T. ~ TDs then 
one must introduce an additional resistance. If Raye is the total 


resistance of the system, R. the resistance due to the drop condensation 


process, and R_ the average "constriction" resistance over the surface, 
cm 


then 
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Pocal ~ ce i oe (G.1) 
or 
ee 7 a ne 7 ao Bc - i 
ee ee a 
Q/A Q/A Q/A (G. 2) 
* 


where ae is now the "true" mean surface temperature. 


The paper of Mikic proceeds to derive an expression for the "con- 


Beriction' 


resistance. Without going into further details concerning 
the derivation, the most usable expression for this resistance will be 


stated to be 
r 


el Y d 
Rom “3 = «G.3) 
(1 - y) 


where ry is the radius of a typical inactive drop, y is the surface 
area covered by inactive drops, and k is the liquid conductivity. 

To evaluate the constriction resistance from equation (G.3), one 
must use some judgment. If one considers inactive drops to be those 
transfering less than 10 percent of the average heat, and chooses the 
typical inactive drop to be approximately on half the maximum drop 
size, then using the results of Figs. 37, 38, 55, and 56, the constriction 
resistance can be computed. 

For low pressure condensation, drops with diameter greater than 
150 microns transfer less than 10 percent of the total heat. According 
to Fig. 37, 56 percent of the surface is inactive. One half of the 
maximum drop size for low pressure condensation is a drop having a 


radius of 750 microns. For atmospheric pressure condensation, Figs. 38 


and 56 show that drops greater than 40 microns in diameter are inactive, 


% e 
Tn is actually not the temperature at the base of drops. There is 


an add Oona A nrouen ne promo a 


this analysis the promoter resistance has been neglected. 
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the inactive area is 68 percent, and the typical inactive drop is one 
of 625 microns. 

The results of substituting the above values into equation (G.3) 
are summarized in Table G-l. 

Table G-1 indicates that even for copper with a very high thermal 
conductivity, the "constriction" resistance is significant for atmos- 
pheric pressure condensation. Mikic showed that for stainless steel 
the constriction resistance would be as much as 80 percent of the total 
resistance. These results show that the effect of thermal conductivity 
of the surface,as exemplified in the "constriction" resistance, cannot 


be neglected. 
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APPENDIX H DISCUSSION OF TEFLON RESISTANCE 


. For chemically promoted surfaces the thermal resistance due to 
the promoter layer is usually so small that it can be neglected. However, 
for Teflon this is not the case. Because the thickness of the physically 


adhering Teflon is at least an order of magnitude greater than for the 
chemical promoters, there is a temperature drop through the Teflon 
which is appreciable. 

If the rate of condensation is to be considered as that due to 
the temperature difference (T. - TD» then one must introduce an ad- 
ditional resistance. The resistance due to the conduction loss across 
the Teflon can be assumed to be in series with the resistance due to 


x 
the dropwise condensation process (see Fig. H-l). Accordingly, 


Frotal ~ Re a 
a (H.1) 
Eos of sete os 1 
S w__s i W 
Q/A Q/A Q/A 
where R is the total resistance, R, is the resistance due to drop- 
total te 


wise condensation, and Bs is the Teflon resistance. The teflon re- 


sistance is given by 


g 
R = aes (H.2) 


where & is the thickness and Ke is the thermal conductivity of the 


Metlon. Substituting (H.2) in (H.1), 


R. + : (H. 3) 


Boal alate 


Ee 





* ° 
The constriction resistance, discussed in Appendix G, has been 
neglected in this analysis. 
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Talking in terms of heat transfer coefficients as opposed to re- 
sistances, equation (H.3) becomes 


1 
Wea ea 


k 
me ay 
area (H.4) 

tC 

Now assuming that the heat transfer coefficient for the drop 
condensation process, hes is the same for Teflon as for chemical pro- 


moters, then the overall coefficient of the Teflon systems is dependent 


x 
on the thickness and conductivity of the Teflon layer. 


The procedure for applying the Teflon on copper described in 
Chapter 5 was re from Du The thickness 


of the resulting layer,after following these procedures, was predicted 
to be 6 x oie inches from past experience of Du Pont engineers. Using 
this thickness and a conductivity of .1 Btu/hr ft °F, the heat transfer 
Mee ificient of the Teflon layer can be computed to be about 20,000 


°F. Assuming that 37,000 Btu/hr ft? °R is the coefficient 


Beu/hr ft 
for atmospheric pressure drop condensation, the overall coefficient will 
be 13,000 Btu/hr £t7°R, This is very close to the measured heat transfer 
coefficient for the Teflon surface of 12,500 Btu/hr Et or. 

Figure H-2 represents the predicted values for the heat transfer 
coefficient for atmospheric pressure dropwise condensation on varying 
thicknesses of Teflon. In all cases the coefficient for the drop con- 
densation itself is 37,000 Btu/hr ft °F. 


Figure H-2 shows that manufacturers must produce commercially 


* | ; 
This will be true provided that the nucleation site density and con- 


tact angle are comparable in both cases. Experimental observation by 


hig ed that the 


chemically promoted surfaces were very similar. 
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available Teflon coatings of thicknesses significantly less than a 
thousandth of an inch in order to make the total resistance (Teflon 


plus drop condensation) less than the resistance for film condensation. 
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APPENDIX I LISTING OF COMPUTER PROGRAMS 


‘The following 7 programs were used extensively during this thesis. 
The programs were written in Fortran IV for use on the IBM 1130 


‘digital computer. 





AGE 


y 4 
13 
14 
15 


16 
17 
18 
19 


20 


21 
22 
23 
24 


26 
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E GRAHAM HEAT TRANSFER DATA REDUCTION PROGRAM=PART I 


THIS PROGRAM WRITTEN TO PERFORM THE FOLLOWING JOBS 
BONVERT MILLIVOLT TC OUTPUT TO DEGREE F 

MAKE LEAST SQUARES FIT STRAIGHT LINE OF TEMP GRADIENT 
COMPUTE HEAT FLUX» SURFACE TEMPsHEAT TRANS COEFFe 


DATA INPUT 


REAL MVLT(10930) 

DIMENSION DIST(10) s TMP (10930) 9 Y¥(30) sPT(30) sSLOP(30) 
DIMENSION HTFX(50) sHTC( 50) sVTMP( 50) »STMP(50) sDTMP(50) 
DIMENSION DEV(10930) 

IR=2 

IW=3 

READ(2991) NORUN 

FORMAT(I2) 

READ(2*92) NOPT 

FORMAT (I2) 

COND=215 

NO=NOPT=1 

READ(IRs9) (DOIST(I) »IT=19NO) 

NN=NOPT 

NM=NORUN 

READ(IR99) C(MVLT(I9J) 9 T=l9NN) sJ=1l5NM) 

FORMAT (7F5e3) 


CONVERSION OF MILLIVOLT TC OUTPUT TO DEGREE F USING FOURTH 
ORDER CURVE FIT OF CU=CONST TC TABLES 

DO 15 J=1sNORUN 

DO 15 I=lsNOPT | 

CALL CCMVT (MVLT(I9J) sTMP( I 9J) sDUMMY) 

CONTINUE 


LEAST SQUARES FIT OF STRAIGHT LINE TEMP GRADIENT 
DO 20 J=1l»sNORUN 

DO 18 I=1»sNOPT 

Y(I) =TMP(I oJ) 

CALL SJHLS(NOPT=1»sDIST*#YsPOINsSLEP) 

PT(J)=POIN 

SLOP (J) =SLEP 

DO 20 I=l»sNOPT 

DEV( IJ) =TMP(I9J) = (PT(J)+SLOP(U)#DIST(T)) 
CONTINUE 


COMPUTATION OF HEAT FLUX» SURFACE TEMPsHEAT TRANS COEFF 
DO 27 J=lsNORUN 

HTFX(J)==12*COND*SLOP (J) 

VTMP (J) =TMP(NOPT 9J) 

STMP (J) =PT(J) 

OTMP(J)= VTMP(J)=STMP(J) 

HTC(J) =HTFX( 4) /OTMP (J) 





Ww Rh) — © oO 
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3 GRAHAM HEAT TRANSFER DATA REDUCTION PROGRAM=PART I 
CONTINUE 


DATA OUTPUT 

DO 32 J=lsNM 

WRITE(1W933) J» HTFX(J) sOTMP (JU) sHTC( JU) 

DO 31 T=lsNN 

WRITE( IW934) IsDIST(I) sMVLT(CI 9J) sTMP(I9U) sDEVII oJ) 

WRITE(IWs35) SLOP(JU) » PT(JU) 

FORMAT (1H1940X'RUN NO'9I3 // LOX'HEAT FLUX=' sF9GelslOX'TEMP DIF='» 


lF4e2slOX*HEAT TRANS COEFF=!sF9el // 9X'POSITION NO'912X'DISTANCE'» 
2 LIX*MILLIVOLT' 913X'TEMP(F)' sl1X'DEVIATION' ) 


FORMAT(I2092F200e3 sF200e2sE2003) 

FORMAT (//1OX'TEMP = 'sF8e2s!' X DIST + '5F8e2) 
CALL EXIT 

ENO 





AGE 


00 


ee - 
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2 GRAHAM HEAT TRANSFER DATA REDUCTION PROGRAM=PART 2 


THIS PROGRAM WRITTEN TO PERFORM THE FOLLOWING JOBS 
FIT STRAIGHT LINE THROUGH DELTA T = HEAT FLUX DATA 
DETERMINE RELATION BETWEEN HEAT TRANSFER COEFFe VS DELTA T 


DATA INPUT 


DIMENSION T(200)s Q(200)sTT(40) 9H(40) 
READ (291) M 

FORMAT(13) 

DO 100 J=loM 

READ(291)N 

READ(292)(T(I) sIT=19N) 
READ(292)(Q(I)sIT=1l9N) 

FORMAT(8F10¢3) 


LEAST SQUARES FIT STRAIGHT LINE THRU DELTA T = HEAT FLUX DATAc 


CALL SJHLS(N sQ9T 9A2B) 
WRITE (393)A9B 
FORMAT(ILH1' INTERCEPT! »9F1l00591l0Xs'SLOPE',E1L5e5//S//) 


DETERMINATION OF H VS DELTA T RELATION 


TT(2)=201 

TT(3)=205 

TT(4)=0207 

TT(5)=0l 

DO ll K=295 
H(K)=(lemA/TT(K))/B 

CONTINUE 

DO 10 K=6939 
TT(K)=TT(K=-1)4001 
H(K)=(LemA/TT(K))/8 

CONTINUE 

WRITE(394) 

FORMAT(*' DELTA T'sl0Xs'HEAT TRANSFER COEF.!) 
WRITE( 395) (TT(I) sH( 1) 91=2931) 
FORMAT(5X9F1l00e2920XsF1l0e02) 
CONT INUE 

CALL EXIT 


END 
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PGE Z GRAHAM COMPUTATION OF DROP GROWTH RATE DISTRIBUTIONS 


C THIS PROGRAM WRITTEN TO PERFORM THE FOLLOWING JOBS 
COMPUTES RELATIONSHIP BETWEEN DROP GROWTH RATE AND DROP DIAMETER 
C COMPUTES RELATIONSHIP BETWEE DIAMETER AND TIME 


mo 


C DATA INPUT 


DIMENSION Y(1001) 
TSAT=6720 
HFG=970.3 
DENL=59428086 
CON=2388 
STEN=2e0904032 
VV=26.78 
QONAM=35000¢6 
S1G=240 
OT=1le 
M=18¢ 
R=1544, 
Q=4e17E8 
O=( (QM) /(6e283#R*¥TSAT) ) #% 05 
B=(HFG**2)/(TSAT#VV) 
HE=778 e*SIG*O*B 
H= 1e8E6/HE 
G= 1252/CON 
E=2e9E~2 
Set 3292 oe XB e*TSATESTEN/S (DENL#®HFG#DT ) 
DMIN=2e%*F 
F=0e20 
| H=020 
WRITE(3926) DTsDMIN 
6 FORMAT(1H1 // 15X92E15e5) 
WRITE(3910) TSATsHFGsVVsDENL sCONsSTEN 
Q FORMAT(//10X92F10 els 2F10e4sF1l0e3 E1104) 
WRITE(3925) SIGsHE»sQONAM 
5 FORMAT(/// 15X93E15¢5) 
WRITE(3919) 
9 FORMAT ( /// 6&X* DIAM! slOX'RATE1]' plLOX'RATE2Z' »10X'RATES3' 
2 LOX*RATES!) 
READ(293) D 
FORMAT(F10e2) 


( COMPUTATION OF DROP GROWTH RATE VSe DIAMETER 


RATE] = E*DT/(DENL#HFG) *((1@F/D) /(G*D+H) ) 
RATE2 = EXDT/(DENL*HFG) /(G*D+H) 

| RATE3 = ExXDT/(DENL*HFG) #((1-F/0)/(G*D)) 
RATES = EXDT/(DENL*HFG) /(G*D) 


| WRITE(3933) DsRATE1] sRATE2 sRATE3 sRATE4 9 DD 
33 | FORMAT(F100294E15059F 1503) 
TF (0e1=-0)2913913 
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3 GRAHAM COMPUTATION OF DROP GROWTH RATE DISTRIBUTIONS 


CONTINUE 

WRITE(3918) 

FORMATCLHI/// 6X'DIAM' s10X#RATEL'510X#RATE2' 510X'RATES' 5 
LOX RATES ' 93X' DIAM SQUARED!) 


COMPUTATION OF DIAMETER VSe TIME 


00 14 I0=1,10 

D=ID*.l 

OD=D#*2 

GG= G¥(e5* (DF) #24 2HF* (DMF) +FHH2HALOG(D=F) ) 
HH = H*(D=-F+FeALOG(D=F) ) 
GI=G#F*#*2#(265+ALOG(F) ) 

HISH*F#(1+ALOG(F)) 

TIMEL= DENL*HFG/ (E*OT) *( GG+HH=-GI =H] ) 
TIMES3=DENL#HFG/(E*OT)*(GG=GI) 

TIME2 = DENL#¥HFG/(CE# OT) #(G/2% (DH 2—4HF HH) +H (DR 2HF) ) 
TIMES = DENL¥HFG/(E* DT) *¥G/2* (DEH 2—4HF #2 ) 
WRITE(32933) De TIMELsTIME2sTIME3sTIME4s00 
CONTINUE 

DO 11 JOD=192192 

D=ID0 

OD=D0**2 

GG= G¥(e5* (DF) HH 2424 F*( DMF) +F RH ZHALOG(DF) ) 
HH = H¥(D=-F+F#ALOG(D=F) ) 
GI=G¥FHH2% (2 e5+AL0OG(F) } 

HI=H*F*(1+AL0G(F) ) 

TIMEL= DENL#HFG/ (E*#DT) *( GG+HH=GI HI ) 
TIME3Z=DENL#*¥HFG/(E*DT)*(GG=GI) 

TIME2 = DENL#HFG/(E*#DT) #(G/2% (D¥R2—4KF REZ) +H*(D—2*F) ) 
TIMES = DENL*¥HFG/ (EX DTI XG/2* (DF 2-44 F ERD ) 
WRITE( 3933) DsTIMELsTIME2sTIME3,TIME4s0D 
CONTINUE 

DO 12 IA=205200910 

D=IA 

OD=0*#2 

GG= G¥(o5%(D=F) #¥242HF %( DMF) +FRHZ*ALOG(D=F)) 
HH = H*(D=-F+FHALOG(D=F) ) 
GI=GH#Fe*2%(2,.5+ALOG(F) ) 

HI=H*F*(1+ALOG(F) ) 

TIMEL= DENL#HFG/(E*DT) *(GG+HH=GI—H] ) 
TIME3Z3©=DENL*HFG/(E*0T)*(GG=GI) 

TIME2 = DENL#HFG/( E* OT) *(G/2* (DER 2—GRFRE2 ) +HE(D—2*F ) ) 
TIMES = DENL#HFG/(E#DT) *¥G/2*# (DEH 2—4HFR*2 ) 
WRITE(3933) De TIMELsTIME2+TIME3sTIME4s0D 
CONT INUE 

CALL EXIT 

END 
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2 GRAHAM INTEGRATION OF HEAT FLUX DISTRIBUTION 


THIS PROGRAM WRITTEN TO PERFORM THE FOLLOWING JOBS 
INTEGRATE THE PREDICTED HEAT FLUX DISTRIBUTION 
COMPARE PREDICTED WITH MEASURED HEAT FLUX 


DATA INPUT 


OIMENSION Y(1001) 

TSAT=67200 

HFG=970.3 

DENL=59e8086 

CON=2388 

STEN=2004032 

VV=26.78 

QONAM=350906 

SIG=2e0 

OT=le 

M=18. 

R=15446 

Q=4e17E8 

O=((Q#M)/(6e283*R*¥TSAT) ) #¥05 

B=(HFG**2)/(TSAT#VV) 

HE=778e*SIG#O*B 

H= 1.8E6/HE 

G= 1252/CON 

E=2e9E~2 

Fao5 #392 o%#Be¥TSATHESTEN/ (DENL*HFGEDT ) 

DMIN=20*F 

F=02e0 

H=020 

WRITE(3926) OTsDMIN 

FORMAT(1H1 // 15X92E15¢e5) 

WRITE(3910) TSAT®HFGoaVV»DENL s CON sSTEN 

FORMAT(//10X92FlOels2FlOe4sFlOe3 sElLLe4) 

WRITE(3925) SIGsHE »QONAM 

FORMAT(/// 15X%93E15e5) 

WRITE(3919) 

FORMAT(/// 1L1IX* DIAM! 9 7X'HEATFLUX ! 9 1LOX'RATIO'» 
OX'FACTOR! 512X'EXP! ) 


INTEGRATION BY SIMPSON'S RULE 


SUM=0 6 

ANT=0 

XU=25006 

ANT=0.0 
READ(293) DeJ 
FORMAT(F100e1912) 
IF(J=1)89894% 


4 READ(2922)29A 


22 FORMAT(F10049E1504) 
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3 GRAHAM INTEGRATION OF HEAT FLUX DISTRIBUTION 


XH=(XU=0)/10. 

DO 102 IT=1lsl1l 

X= D+XH*(1=1) 

YC I) SE*ARXHH( 22-2) *OTH# (Lem F/X) /( GHX+H) 
ANT=ANT#Y(1) 

DO 103 1=2910 
ANT=ANT#Y(1) 

DO 104 1=291092 
ANT=ANT#Y( I) 4Y(1) 
ANT=ANT*XH/3 6 
SUM=SUM+ANT 
RATIO=SUM/QONAM 
XU=D 

IF (J=-1)88588544 


DATA OUTPUT 


FORMAT (F15029E15059F15049E15049F 1503) 
WRITE( 39333) DeSUMsRATIOsAsZ 

IF (201"D)299913 

WRITE(39333) DsSUMsRATIO 
TF(001°0)299913 

CONTINUE 

CALL EXIT 

END 
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2 GRAHAM DETERMINATION OF DROP SIZE DISTRIBUTION 


THIS PROGRAM WRITTEN TO PERFORM THE FOLLOWING JOBS 
COMPUTES COMMON LOG OF NO OF DROPS AND SIZE OF DROPS 
FITS STRAIGHT LINE THRU LOG NOe VS LOG DIAM 


DATA INPUT 


DIMENSION 0( 100) 200(100) »AN(100) sANN(100) 
READ( 291) N 

FORMAT(13) 

READ(292)(D(I)sI=loN) 

FORMAT(10F8e0) 

READ(293)(AN( I) sIT=loeN) 

FORMAT (8E10e¢3) 


COMPUTE LOG NOe AND LOG DIAMe 


DO 4 I= 159N 

DD(I)=ALOGT(D(1)) 

ANN(T)=ALOGT(AN(TI)) 

CONTINUE 

FITS LEAST SQUARES STRAIGHT LINE THRU LOG NOe VS LOG DIAMe 
FORM OF EQUATION IS N=A#D##2 


CALL SJHLS(NsDDsANNs POINT» SLOPE) 
PINT=1l0e**POINT 

WRITE(397) 

FORMAT(1H1ls INTERCEPT! 910Xs'SLOPE!') 
WRITE(395)PINTsSLOPE 
FORMAT(E1545910XsE15-45) 

CALL EXIT 

END 
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AGE 2 GRAHAM CONVERSION OF THERMOCOUPLE OUTPUT TO DEGREES F 


SUBROUTINE CCMVT (TMVoTFoTR) 
THIS PROGRAM WRITTEN TO PERFORM THE FOLLOWING JOB 
THIS SUBROUTINE CONVERTS THE OUTPUT OF A COPPER CONSTANTIN 
THERMOCOUPLE» IN MILLIVOLTSs TO TEMPERATURE IN DEGREES 
FAHRENHEIT» OR DEGREES RANKINEe 
IF(TMV45e379) 19293 
1 TF==3100+( TMV4+50379)/00009 
GO TO 50 
2 TF=-3106 
Ge 10 50 
3 TFCTMV4+52284) 49596 
G TF==3000+(5e2844+TMV) /020095 
GO TO 50 
5 TF=-3006 
GO TO 50 
6 IFC TMV44ell1l) 79899 
T TFE= 061336559E4— ,LLI7L49BE4#ETMV= »4090264E3*%TMV¥#2—= 
1 eFO0764B86E2Z*TMV*¥*E3Z—@ e34611LIBELXTMVE#G 
GO TO 50 : 
| 8 TF=—2006 
GO To 50 
9 ITF(TMV422.559) 10911912 
10 TF= 05024062E2+ s6708685EZ* TMV e654206E1* TMV¥¥2 
Ll + el3A66154E1L*XTMVHE3+ 346861 1LEMLETMVEHG 
GO TO 50 
| ll TF=-100. 
GO TO 50 
~~ 12 TFCTIMV400e67) 13914915 
~—6©13 TF= 63235749E2+ e4761L4OLE2Z*#TMV= e81BO0 7X TMV¥2 
1 + 6 2912321*TMV*#3= 066431592E-2* TMV ERG 
. GO TO 50 
14 TF=0.20 
GO TO 50 
1 aE5 IF(TMV=1.2517) 16917918 
—16 TF= 0 3200279E2% e4679051E2*TMV= 061382542E1* IMVR*2 
Ll + e9235143 E-1LETMVERZ= o2356772EMLETIMV ERG 
Be TO 50 : 
17 TF=1006 
GO TO 50 
18 IF(TMV=-34967) 19920921 
19 TF= 063562773E2+ e4092755E2*TMV+ 01955488E]1 
Ll *TMV#*#2— og 742077B*#TMVEHS4 oe TODZS66EM1LX TMV EES 
GO TO 50 
20 TF=2006 
GO TO 50 
21 TF(TMV=62647) 22923924 
22 TF= 04580235E2+ e3610Z307TE2Z*TMV+ 10644181 *TMV¥*2 
Ll = ep 29T4O3Z3* TMVER3+ slLGOSLZBIDEMLETMV ERG 
GO TO 50 
23 TF23006 


[rx x OX Fw 
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PGE 3 GRAHAM CONVERSION OF THERMOCOUPLE OUTPUT TO DEGREES F 


GO TO 50 
(24 TF(TMV=92525) 25926927 
25 TF=— eL375061E3+ e1323431EB*TMV= 1745894 E2HTMVHR2Z+ 
aq elLGOB9SBS4ELBHIMVEHS meh 33329 LEMLETMV HG 
GO TO 50 
26 TF=4006¢ 
GO TO 50 
27 TFC TMV"126575) 28929930 
28 TF= e69992Z1LE2+ 03403 852E2*TMV+ o2749543%TMVE*2 
lL = e2620201LEM1L#TMV#ES+ 2642380 36E 3H TMV ERG 
GO TO 50 
29 TF=5006¢ 
GO TO 50 
30 IFC IMV"15e773) 315932933 
31 TF= e31lOTOBIEZ—= e3299984E2Z*TMV+ oe TZ06085E1L*TMVEX2 
1 = e3568681*¥TMVEXS+ ec6394187TER2ZHTMV ERG 
GO TO 50 
32 TF=6006 
GO TO 50 
33 IF(TMV"19.1) 34935536 
34 TF= e2632777E3=— e929301LELRIMV+ e38965S6ELHTMVe82 
L = el619IO59#TMV¥#Z+ 2241077 2EM2*# TMV HRY 
GO TO 50 
35 TF=7006 
GO TO 50 
36 TF=— 6427098E4+ e9391044GE3*TMV= 270486 7EZHTIMV#R2 
Ll + e2430854E1LXTIMVER3Z—= e3147985E—1L#TMVEHG 
~50 TR=TF+459.69 
RETURN 
END 
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2AGE Z SAKHUJA’ THIS SUBROUTINE FITS A FIRST ORDER CURVE 


SUBROUTINE SJHLS (NoXsY¥oAoB) 
DIMENSION X(1).5 Y61)sAA(2), BB(292) 
SET SUMS TO ZERO 
YY=0.0 
XX=020 
XY=040 
X2=040 
E FORM SUMS 
Bor i f=isNn 
XX=XX +X(1) 
YY=YY+ Y(J) 
XYSXVY#X(T)*Y( J) 
X2=X2+ X(1)**2 
1 CONTINUE 
SETUP THE INPUT FOR SIMQ 


«3 


AA(1)=YY 
AA(2)=XY 
BB(1lsl)=N 
BB(1l»2)=XxX 
BB{(251)=xXxX 
BB(292)=xX2 


SOLVE MATRIX 

SIMQ 1S A 18M SCIENeSUBe PACKAGE WHICH SOLVES SIMe LINe EQSe 
CALL SIMQ(BBsAAs2 9K) 

A=AA(1) 

B=AA(2) 

RETURN 

END 


Se ey 
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APPENDIX J DROP COUNTING DATA AND PROCEDURE 


In each photograph, the number of drops of a certain size was 
counted. The average number of drops per photo of this size was ob- 
tained by summing the number counted on each photo and dividing by the 
number of photos. 

The average number of drops per photo was then converted to the 
average number per square centimeter by multiplying by a "magnification 
factor’. The magnification factor equals the reciprocal of the actual 
surface area in the field of view of the photograph in square centi- 
meters. The number so obtained represents one data point on the plots 
of AN versus diameter (Fig. 33). If the same size drop was counted at 
two different magnifications, the average of the two results was ob- 
tained before plotting. 

The best straight line was fitted by least squares technique through 
the data of log AN versus log D in each drop distribution region. The 
resulting form of the drop distribution equation in each region was of 
the form AN = np “1, The equation for N versus D was obtained by 


dividing AN by the band width AD. The data presented is for a band 


PaageAD = .4D. Thus, N = ap 2 Dd”. 
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fe.) Drop Counting Data for Low Pressure (Lae = 88 °F) Condensation 


J.1.1 Magnification = 400X. 150 photos taken. 25 photos counted 


Photo No. Number of Drops in Photo of Diameter 
50u 2511 10u 
iL 0 6 10 
2 0 3 30 
3 1 a 8 
4 0 6 11 
5 0 16 5 
6 iL Z is 
7 Z 2 12 
8 0 0 0 
2 0 2 
10 1 2 
11 2 3 20 
i2 0 1 4 
13 ue 0 12 
14 0 2 20 
15 0 0 0 
16 1 5 La 
7 2 4 21 
18 0 6 38 
9 2 0 15 
20 0 0 0 
21 0 0 0 
De 1 3 10 
23 0 3 ZZ 
24 0 0 0 
25 0 1 1 
Total 14 70 278 


mn 

ov 

nO 
“ 
_ 
a 
_ 


Average 





-158- 


el. 2 Magnification = 200X- 100 photos taken. 25 photos counted 


Photo No. Number of Drops in Photo of Diameter 
50 30u 20u 
it 7 1 11 
2 Zz 4 17 
3 11 8 14 
4 il 4 13 
5 i 5 11 
6 i 7 10 
7 6 9 30 
8 2 11 ne 
9 15 1? 12 
10 8 15 10 
11 6 12, 16 
12 3 6 10 
13 6 
14 1 0 0 
15 4 13 1S 
16 3 16 37 
17 4 ] 8 
18 0 13 18 
19 5 11 L5 
20 4 5 7 
ZA 6 1 13 
»? 2 19 28 
23 ZZ TZ 12 
24 5 10 12 
Z5 L 4 33 
Total 125 241 379 


Average 5 20 LSa2 
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J.1.3 Magnification = 120X . 55 photos taken. 25 photos counted 


Photo No. Number of Drops in Photo of Diameter 
. 200 100u 50u 
1 1 3 8 
Z 0 Ht 14 
3 0 3 16 
4 3 6 10 
5 2 0 ia 
6 1 Zz 11 
7 0 0 0 
8 0 0 15 
9 0 2 14 
10 0 0 0 
11 Zs 3 15 
2 2 3 15 
13 3 5 3 
14 0 i 9 
15 1 1 10 
16 0 2 Hz 
17 1 5 11 
18 0 3 14 
19 Z 6 8 
20 1 3 As 
21 2 4 10 
Ze 1 2 Ie 
23 2 1 11 
24 0 2 18 
Zo it 3 6 
Total 23 67 269 
Average 92 2.68 Oe 
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J.1.4 Magnification = 60X . 40 photos taken. 25 photos counted 
ee ee eee Oem United 


Foto No. Number of Drops in Photo of Diameter 
500u 300u 200 100 

1 3 2 0 12 
2 0 0 1 6 
3 2 0 4 9 
4 0 1 3 10 
5 1 3 2 11 
6 1 il 3 14 
7 0 0 1 6 
8 0 1 il 5 
9 0 0 0 0 
10 0 6 3 LS 
aL 0 1 5 13 
12 0 0 4 ell 
13 0 2 5 6 
14 0 1 2 33 
15 2 0 2 15 
16 0 il 7 20 
17 il 1 Z 3 
18 0 5 10 7 
19 2 Z 2 8 
20 0 0 34 
za 2 0 11 
22 0 4 11 2 
23 | 2 0 1 21 
24 0 3 10 6 
Z2 0 6 3 8 

Total 16 40 81 293 

»64 1.6 ee On), 


Average 
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J.1.5 Magnification = 10X. 24 Photos taken. 24 photos counted 


Photo No. Number of Drops in Photo of Diameter 

2000u 1000u 
i 1 iy 
2 3 15 
3 2 0 
4 it iy 
5 0 9 
6 0 14 
7 0 12 
8 it 14 
9 0 5 
10 2 16 
11 0 3 
eZ 0 
13 é 
14 0 
15 1 12 
16 ir 11 
17 3 6 
18 0 0 
19 0 8 
20 3 / 
Di 0 7 
22 i 16 
23 2 19 
24 0 5 

Total 23 242 


10.1 


\O 
oN 


Average 
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J.1.6 Magnification 


Photo No. 


wo won Do UN FEF WH NY FF 


Pee Ee Fe 
Un & WH F O 


Total 


Average 


Dw 


Sloe 


15 photos taken . 


15 photos counted 


Number of Drops in Photo of Diameter 


2500 
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J.2 Drop Counting Data for Atmospheric Pressure Condensation 
‘J.2.1 Magnification = 200X. 100 photos taken. 50 photos counted 


Photo No. Number of Drops of Diam. Photo No. Number of Drops of Diam. 


50u 30 20u 50u 30u 201 
il 3 5 7 26 1 3 1 
2 4 24 48 oy 4 6 31 
3 7 7 22 28 6 5 10 
4 4 19 29 2 2 13 
5 2 8 13 30 9 1 13 
6 0 det 130 31 0 0 0 
7 14 9 9 32 4 13 15 
8 4 12 33 4 8 31 
9 1 21 34 , 7 8 
10 i 10 24 35 y 8 21 
tel 3 9 25 36 2 11 7 
12 5 15 37 1 2 12 
13 4 20 38 4 13 7 
14 2 12 21 39 4 13 22 
15 i 5 28 40 6 4 15 
16 2 12 23 41 1 3 9 
17 3 4 14 42 10 6 12 
18 2 19 4 43 it 9 7 
19 7 11 13 4 3 2 18 
20 2 10 30 45 3 2 20 
21 2 60 57 46 9 18 12 
22 3 6 14 47 4 2 11 
23 2 7 31 48 5 4 10 
24 3 5 8 49 7 6 19 
25 3 8 15 50 5 2 10 
Total 84 275 640 Total 106 160 550 


Average 3.4 11.0 2520 Average ee 6.4 14.3 
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J.2.2 Magnification = 120X. 50 photos taken. 50 photos counted 


Photo No. Number of Drops of Diam. Photo No. Number of Drops of Diam. 
200. 100. 50 200 1001 50u 

1 3 3 5 26 0 0 52 
2 0 3 2 Jal 0 0 16 
3 1 4 1s, 28 - 0 1 18 
4 1 3 Zo 1 2 5 
2 0 i 30 0 3 25 
6 a 3 18 ai 1 2 8 
7 0 4 14 oz i 6 IE 
8 0 5 24 33 4 3 14 
9 0 6 2 34 1 5 8 
10 0 aa 9 35 0 5 abaL 
11 il 6 9 36 il 13 9 
2 a 0 11 37 0 iz 3 
tS 3 3 9 38 0 Zz 26 
14 0 9 5 3 0 4 3 
iS 2 2 7 40 0 0 11 
16 1 5 8 41 3 dal 
iE / 0 1 2 42 4 it 7 
18 2 2 16 43 i 4 
19 Z 3 5 44 0 17 
20 4 3 14 45 0 1 20 
21 2 iL 46 0 1 18 
ne 1 4 9 47 2 7 10 
23 4 0 2 48 0 Zz 7 
24 2 8 49 0 di: bil 
Zz 0 0 | 50 3 5 9 

Total oil 88 236 Total Ze 98 304 

Average ts 2 S20 9.4 Average ag 340 1735 Il 





—— 
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J.2.3 Magnification = 60X. 40 photos taken. 40 photos counted 


Photo No. Number of Drops of Diam. Photo No. Number of Drops of Diam. 
500u 300u 200u 100pu 500u 300n 2200p 100, 
i} i us 0 i 21 0 0 5 22 
2 0 > i> 11 Ze 0 0 0 0 
3 0 0 0 3 Zs Z 0 3 13 
4 0 0 0 45 24 1 4 2 3 
5 0 0 3 50 Zo 0 0 ii 50 
6 il 4 0 8 26 if 1 0 10 
7 il 3 1 16 zy 4 3 i 15 
8 0 0 0 0 28 0 Z 0 42 
2 0 Z Ul 6 Zo 0 Z 5 1 
10 2 0 3 1S: 30 3 Zz 0 1/7 
ee 1 16 15 31 0 0 0 0 
2 EZ i 3 7 Sy 0 0 1 45 
3 4 1 1 13 a3 0 Ai 0 10 
14 Z 2 0 za 34 1 1 us 9 
15 0 iy 9 8 35 0 0 3 4 
16 2 0 1 29 36 0 0 4 10 
7 Z 2 1 Tt oy 4 1 1 7 
18 4 1 1 7 38 0 0 0 0 
19 i 6 2 13 39 0 0 0 0 
20 0 6 9 8 40 2 Z 5 9 
Total Zc 48 72 Z95 18 19 42 Zig 


Average 1.1 Zn. 0 14.7 a 1.0 2.1 eESsc 
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J.2.4 Magnification = 10X. 24 photos taken. 24 photos counted 


Photo No. Number of Drops of Diameter 
2000 1000 

1 i 3 

2 0 7 

3 0 5 

4 1 9 

3 0 0 

6 i 0 

7 i 4 

8 1 3. 
9 0 iz 
10 0 16 
ml 0 

2 0 
13 0 Z 
14 1 2 
15 0 9 
16 1 7 
17 0 3 
18 0 8 
19 0 6 
20 0 14 
21 Z 7 
22 ae! 2 
23 0 4 
24 0 0 
Total 10 149 


Gnz 


= 
NO 


Average 
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J.2.5 Magnification = 5X. 16 photos taken. 15 photos counted 


Photo No. Number of Drops of Diameter 
2500 
i 3 
2 1 
> 0 
4, 3 
5 Z 
6 2 
7 il 
8 5 
9 1 
10 2 
ey Lh 
12 2 
is L 
14 1 
es 4 
Total 2] 


Average io 
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DIST. 
NUMBER 


L-5 


LA4 


1 T3 


L-2 


Lel 


L-0 


DIAMETER 


RANGE 


500 < D < 3000 


30 < D < 500 
10 < D < 30 
1<D< 10 
D<1 
1<D< 10 
D<l 
L << 10 
D<l 
1<D< 10 
D<il 
1<D< 10 
D<1 
1<D< 10 
D<l 
1<D< 10 
D<i 
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AN 


5.73x10°D 
Balyx LO De 
8.72x10°D 
8,72x10°D~ 
8,72x10°D~ 
8.72x10°D_ 
8.72x10° 


2.00x10° 


le 


91 


83 


1e47 


4&7 


le%&7 


-2.00x10°D~ 


1.50x10°D~ 


1.50x10° 
9,00x10°D™ 
9,00x10°* 
6.00x10°D_ 
6.00x10" 
2.98x10"* 
2.98x10" 


4 7 


& 8 


30 


TABLE 12 DROP DISTRIBUTIONS FOR LOW PRESSURE CONDENSATION 


N = AN/AD 


me 
1.43x10’D , 


7.92x10°D~ 


Ze 
2.18x10°D : 


2.18x10°D~ 


2e 
2.18x10°D~ . 


2.18x10°D~ 


2.18x10°D ; 


5.00x10°D_ 


1 
5.00x10°D" 


3.75x10°D~ 


3.75x10°D 
2.20x10°D 
2.20x10°D 
1.50x10°D_ 


1.50x10°D~ 
7.36x10°D~ 
7.36x10°D~ 
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TABLE 13 DROP DISTRIBUTIONS FOR ATOMOSPHERIC PRESSURE CONDENSATION 


DIST. 
NUMBER 


A-6 


 AS5 


ARS 


A=-2 


A=-1 


A-0 


DIAMETER 


RANGE 


500 < D < 2500 


10 < D < 500 
ims D < 10 
D<1l 
iD < 10 
D<1l 
1 <D< 10 
D<l 
1 <D< 10 
D<1 
r- p) < 10 
D<1 
1 <D < 10 
D<1 
1 < D < 10 
D<1l 


AN 


2.07x10°D~ 
2.11x10°D~ 
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2.11x10°D 


2.11x10°D~ 


-2.11x10° . 
1.00x10°D~ 
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280 
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69 
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5.30x10°D 
5.30x10°D~ 


Ze 
2.50x10°D~ 
1¢9 


2.50x10°D~ 


Ze 
1.75x10°D~ 
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209 10 
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